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Abstract 


Changes  in  the  stiffness  and  strength  of  notched  quasi- isotropic 
graphite-epoxy  laminates  were  recorded  and  related  to  the  fatigue  damage. 
Two  different  laminates  were  considered  and  the  effects  of  stacking  sequence 
were  compared.  Nondestructive  testing  techniques  such  as  X-radiography , 
moire  technique,  acoustic  emission,  deply  technique,  and  stiffness  change 
were  performed  to  observe  damage  development.  Results  show  that  the  mechan¬ 
ical  response  and  the  fatigue  damage  depend  strongly  on  the  stacking  sequence 
of  laminates.  In  general,  residual  strength  increased  remarkably  for  both 
laminates  due  to  stress  redistributions  while  the  continuous  stiffness 
change  curve  is  typical  for  each  laminate  and  reflects  damage 
characteristics.  Buckling  effects  as  well  as  matrix  cracking  and  delami¬ 
nations  contribute  to  stiffness  changes. 
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1.  INTRODUCTION 


The  study  of  hole  related  fatigue  damage  in  materials  is  motivated  by  Mi- 
need  to  include  holes  into  structures.  For  example,  in  an  airplane,  oven  * 
critical  plate  or  shell  may  have  to  provide  access  from  one  side  to  anoth-v 
for  electric  lines,  fuel  pipes,  or  other  supplementary  equipment. 

As  of  to-date,  the  knowledge  of  the  fatigue  response  of  notched  composite 
laminates  is  much  less  complete  than  that  of  unnatche!  laminates.  The  rrn  ■ 
reason  is  that  fatigue  damage,  residual  strength,  and  stiffness  predict]' 
models  (1,2]  for  the  unnotched  situation  are  based  on  uniformities  of  dam a • 
that,  are  reflected  by  terms  such  as  ’characteristic  damage  state  (CDS)',  ; 
instance.  These  uniformities  of  damage  do  not  apply,  however,  for  notch 
specimens.  An  idea  of  the  complexity  of  the  fatigue  behaviour  of  notched 
cimem;  nay  be  given  by  studying  an  X-ray  photo  of  damage  around  the  hole  of 
notched  specimen  after  h • gh  fatigue  loading. 

The  objectives  of  the  present  work  are  to  doterriae  "h(i  nature  of  dr: 
induced  in  graphite  epoxy  composite  laminates  will*,  c'cte-  holes  by  rye  • 
tensile  loading  and  to  establish  the  influence  ot  such  damage  on  \ : 
strength,  stiffness  and  life  of  laminates. 

The  complexity  of  damage  in  composite  materials  regu  i  re~  ditfer^at: 
ni'pn>s  for  complementary  insight.  Sendecky  j  (3|  *'r'i\v  \  ♦  ha»-  t  he  X-ray  t'*  . 
nine**  makes  damage,  namely  matrix  cracks  and  del  j  i«  -  ■.  >  -  •».«.  .  isible  h-r 
rest  r  j t  ed  t.o  observing  the  x-y  distribution  of  -lamag*'  •  n.  1  v  n-,.  d*n1y  i  • 
nique  complements  the  X-rey  photography  because  oi  >M  "v  *o  idem- • 

sp-'ijfic  planes  of  de  lamination .  X-ray  photog’- adv  ;  *  h«»  de^iy  techn’u 
are  also  opposite  in  the  sense  that  the  former  is  a  n  ••  .i  •  ire  vest  (V; 
method  while  the  latter  inevitably  destroys  the  test  s;-'  or>  n  . 

Tfie  advantage  of  the.  ultrasonic  Oscan  technique,  <  -■  w  r  -  X-  riv  phot'-, 

laphy,  is  that  its  sensitivity  is  not  dependent  -n  -■  ■!,u:m  medium,  t 
enabling  it  to  detect  de  1  aminat  ion  -  I  ike  flaws  -- a  i  w-erd  ’  \-.v  undamaged  mo.  •• 
rial.  Damage  events  during  the  loading  process  -.an  i-  e.  ;  ,.-j  hy  .r  er.v.  1 
emission  recording.  It  yields,  in  comb  ma:  i  %:»  with  ,  •'  n:i  i  t  nneouc 

recorded  load-strain  curve,  information  on  the  sounep-f  md  he  iuterac.'  ’  •  > 
effects  of  subsequent  damage  events. 

Sima*  the  strain  field  in  notched  specimens  mu-un  I  !  e:  .  ;  j;  difficult 

decide  what  strain  measurement  technique  giv'-s  post.  valui!* 

information.  However,  the  moire  technique  (4f  provides  insight  into  the  c-  • 
p’eM»  nonhomogeneous  in-plane  strain  field. 

A  recently  established  investigative  method  is  vhrt  as*'  of  stiffness  ch  i  .; 
record ’tig  as  an  indicator  for  damage  state  or  strength  values  [  ^  j  . 

F.T.  b.'irpono.sch  i  and  W.W.  St inchcomb  fb|  pointed  net  ■  i-.av  a  qu.is  i  -  i  sot  roi-  . 
laminate  exhibits  a  different  response  than  itc  "nmt  ;!"•■  ah  1  ami n a • 
For  instance,  the  response  of  a  (  0 , 90  ,  +4  5  , -*o  ]  <;  laminate  cannot,  he  predict*-: 
bv  knowing  the  behaviour  of  the  [0,90]s  and  of  the  |4S,-A:5Js  laminate.  Stat  ¬ 
ic  strength  predicting  models  for  notched  specimen:-  v-*r'-  developed  that  iib 
not.  require  a  laminate  stress  analysis.  Instead,  the  n  etched  and  the  unno*  - 
r.hrd  situation  are  compared.  J.M.  Whitney  and  H.Y.  K;m  17]  publirdiod  :• 
results  of  a  test  series  which  exhibit  that  notched  failure  mechanisms  d i f t 
from  the  unnotched  situation  They  investigated  the  u  fluenm?  of  two  cliff--' 
ept  stacking  sequences,  one  of  them  producing  tensile,  >  he  other  prodn-  : : 
rnnipir-KKivn  interlaminar  normal  stress  tbrot-ghrut  the  ihidness  of 
straight  edge  on  both  notched  and  unnotched  specinw-ns  Die  data  indicate  *h  - 
unn^f  died  tensile  strength  is  reduced  by  the  straight  *-dge  effect  s  an-1 
strongly  dependent.  on  the  stacking  sequence;  while  die  strength  of  notth- 
spec  inv-i:s  is  independent  of  the  straight  edge  efforts  ■  inr°  tensile  tail::-  • 
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is  -nit  fated  by  the  stress  concentrat ion  at  the  notch  before  straight  edge 
it*'  1  an.inat  ion  comes  into  effect. 

Other  investigations  of  the  fatigue  behaviour  of  notched  specimens  were  con¬ 
ducted  by  Sendeckyj  [31,  Stinchcumb  and  co-workers  [2],  and  Whitcomb  [8}.  The 
work  most  comparable  to  the  present  one  is  that  of  Whitcomb,  who  described 
the  damage  state  and  the  mechanical  response  of  four  different  notched  lami- 
nitos  after  tens  ion-tension- fat igue  loading  and  related  delamination  onset 
to  the  results  of  a  FEM  stress  analysis.  The  load  levels  of  two-thirds  of  the 
initial  strength  of  the  laminates,  respectively,  were  low  enough  to  let  the 
specimens  survive  a  lifetime  of  ten  million  cycles.  At  this  point,  the  dam¬ 
aged  specimens  were  investigated  using  destructive  and  nondestructive 
methods.  The  residual  strength  of  all  laminates  was  equal  to  or  higher  than 
the  initial  strength  of  specimens,  respectively.  Stiffness  degradations  were 
found  to  be  in  the  range  from  zero  to  minus  10  percent.  Whitcomb  also  shows 
figures  comparing  initial  stress  distr ibu :ions  with  de laminat ions  at  differ¬ 
ent  locations  through  the  thickness  on  the  hole  that  suggest  that 
interlaminar  normal  stress  and  shear  stress  distributions  govern  the 
locations  of  first  fatigue  de  laminat  ions .  Later  in  the  lifetime,  the.  direc¬ 
tion  of  delamination  growth  can  be  changed  by  an  altered  stress  distribution. 
In  the  present  study,  two  laminates  with  relatively  little  straight  edge  ini¬ 
tiated  damage  were  chosen  in  order  to  isolate  notch  effected  damage  from 
straight  edge  driven  damage  as  much  as  possible. 

The  damage  induced  by  static  and  by  fatigue  loading  was  observed  by  nonde¬ 
structive  inspection  and  other  methods.  Quasi-static  tests  were  performed  to 
compare  the  development  of  damage  under  cyclic  with  that  under  static 
loading.  Experimentally  observed  data  were  interpreted. 


2.  EXPERIMENTAL  PROCEDURES 


2.1  Material  and  Specimens 

The  material  in  this  study  is  T300-5208  graphite  epoxy.  The  Specimens  were 

ten  inches  long,  1.5  inches  wide  and  had  a  0.375  inch  diameter  hole  in  the 

center.  The  two  stacking  sequences  were: 

[  0,  90,+45,-45)s  type  A 

[+45,  90,-45,  Ojs  type  B 

For  comparison,  static  strength  tests  were  also  performed  on  unnotched  speci¬ 
mens  with  the  stacking  sequence: 


(+45,-45,  0,  90 ] s  type  C 


2.2  Fatigue  Tests 

All  fatigue  tests  were  run  in  the  tens  ion- tens  ion  mode  with  a  load  ratio  of 
R~0.1  and  a  frequency  of  10  Hertz. 

Specimens  were  tested  at  several  different  load  levels  evaluated  as  percent¬ 
ages  of  the  mean  tensile  strength  of  five  samples.  The  load  levels  spanned 
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Cho  range  from  70  to  95  percent;  i.e.,  from  long- 1  i  fe- load  levels  to 
short- life- load  levels. 

Several  specimens  were,  run  to  failure  to  obtain  1  i  f  et  ime-over- load  (s-n) 
curves,  others  were  monotonical  ly  loaded  to  failure  after  a  certain  number  of 
cycles  to  obtain  the  typical  strength  versus  load ing-h istory  curves .  Many  of 
these  specimens  were  subjected  to  nondestructive  inspections,  such  as 
X-rad iography  or  moire  tests,  at  selected  numbers  of  cycles. 


3.  RESULTS 

3.1  Initial  Tensile  Strength 

Tensile  strength  tests  were  performed  on  five  specimens  of  each  type  of  lami¬ 
nate,  including  the  unnotched  laminate. 

In  Figure  I,  the  failure  load  of  each  specimen  is  represented  by  the  length 
of  a  vertical  line.  The  average  failure  loads  of  A,  11  and  0  specimens  wr. 
2620,  2240,  and  4680  lbs.,  respectively  as  indicated  by  the  solid  horizontal 
1 ines . 

The  lower,  dashed  horizontal  lines  indicate  the  averaged  load  levels  for 
which  major  discontinuities  in  the  strain- load  curves,  as  well  as  the  first 
strong  acoustic  emissions  during  the  loading  ( rocr-s,  occurred. 


Figure  1 .  Static  Strength 


The  percentages  refer  to  the  average  failure  load  of  each  type  of  laminate. 
The  large  difference  between  the  strength  values  of  notched  and  unnotched 
specimens  is  due  to  the  stress  concentrating  effect  of  the  hole.  Yet,  it  must 
he  noted,  that  the  strength  of  the  Ospecimens  is  also  reduced  by  straight, 
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free  edge  effects.  But  the  strength  values  indicate  that  a  strength  reduction 
factor  for  the  hole  roust  be  much  larger  than  one  for  the  straight  edge. 
Furthermore,  since  edge  effects  in  general  (whether  for  curved  or  straight 
edges)  decrease  rapidly  with  increasing  distance  from  the  edge,  the  effects 
of  the  straight  edge  and  the  curved  edges  of  the  notched  laminates  are  inde¬ 
pendent.  In  other  words,  we  may  assume  tint  for  the  notched  specimens  used  in 
this  study,  the  effects  of  the  straight  edge  do  not  yield  an  additional  con¬ 
tribution  to  the  reduction  of  static  strength  caused  by  the  hole. 

It  is  important  to  note  that  the  highest  measured  residual  strengths,  related 
to  the  notched  cross-section,  of  A-  (6f>  S  ksi)  and  B-  (71.4  ksi)  specimens 
approach  the  tensile  strengths  of  C-specunens,  thus  indicating  that  fatigue 
damage  redistributes  local  stresses  to  the  extent  of  reducing  the  stress  con¬ 
centration  to  the  hole. 


3.2  Initial  Stiffness 

Stiffness  data  were  taken  from  the  linear  part  of  the  load-strain  curves 
recorded  during  the  tensile  strength  tests  by  means  of  a  1.0  in.  extensomet- 
er  centered  with  respect  to  the  hole.  F.ach  vertical  line  in  Figure  2  repres¬ 
ents  the  stiffness  of  a  specimen.  Average  stiffness  values,  based  on  the  area 
of  the  unnotched  section,  are  indicated  by  the  horizontal  lines. 


ABC 


Figure  2.  Static  Stiffness 

The  stiffness  values  of  A-  and  B-specimens,  based  on  the  cross-sectional 
area,  are  lower  than  the  stiffness  values  determined  for  unnotched  C  type 
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specimens.  This  demonstrates  the  large,  local  deformation  field  associated 
with  the  local  stresses  due  to  the  hole. 


3.3  S-N  Data 

A  lifetime  of  more  than  one  million  cycles  under  the  condition  of  a  sinusoid¬ 
al  load  with  a  ratio  omin/omax  =  0.1  and  a  frequency  of  10  Hertz  can  be 
expected  at  an  80  percent  load  level  (2093  kips)  for  A-spocimens  and  a  70 
percent  load  level  (1707  kips)  for  B-specimens.  Load- log  life  curves  for  the 
two  laminates  are  shown  in  Figure  3.  A-specimens ,  cycled  on  a  93  percent  load 
level  (2485  kips),  and  B-specimens,  cycled  on  an  85  percent  load  level  (2093 
kips)  will  survive  hundred  thousand  cycles. 

3? 
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Figure  3.  S-N  Curves 

Two  A-specimens  had  an  average  lifetime  of  1.1  million  cycles  at  the  80  per¬ 
cent  load  level  while  the  B-specimen,  cycled  at  the  70  percent  level,  had  a 
residual  strength  of  137  percent  of  the  virgin  strength  after  one  million 
cycles.  The  residual  strength  data  for  type  B  laminates  suggest  that  residual 
strength  values  of  this  magnitude'  occur  during  the  final  decade  of  logarith¬ 
mic  life. 


3.4  Stiffness  Degradation 

Figure  »*  on  page  6  shows  stiffness  data  for  type  B  laminates  cycled  at  the.  85 
percent  1  eve  1 . 

The  two  stiffness  curves  show  the  effect  of  gage  length  on  the  stiffness  val¬ 
ues.  A  clip  type  gage  extensometer  centered  over  the  hole  was  used  to  measure 
displacement  over  a  1.0  in.  gage  length.  A  i)CDT  extensometer  centered  over 
the  hole  measured  displacement  over  a  3.625  in.  gage  length.  The  stiffness 
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change  values  based  on  the  short  extensometer  data  are  larger.  However,  if 
the  damage  zone  spreads  outside  the  1.0  gage  length,  the  stiffness  change 
values  determined  from  the  short  extensometer  are  inaccurate. 


Figure  4.  Stiffness  and  Residual  Strength  Versus  Life 

The  process  of  damage  development  around  the  hole  of  a  type  B  laminate  is 
indicated  by  the  stiffness  curves  shown  in  Figure  4.  The  limits  of  the 
regions  I,  II,  and  III  are  somewhat  arbitrary,  but  are  chosen  to  reflect 
transitions  in  the  damage  growth  process  as  represented  by  stiffness  change. 
In  region  I,  the  stiffness  decreases  at  a  decreasing  rate  as  small  matrix 
cracks  develop  near  the  hole.  In  region  II,  the  stiffness  degrades  linearly 
with  log  cycles  and  the  matrix  cracks  in  the  zero  and  45  deg.  plies  extend 
away  from  the  hole  and  de laroinat ions  form  in  the  region  damaged  by  matrix 
cracks.  Region  III  begins  as  the  stiffness  degradation  rate  increases.  In 
this  region,  the  sharp  decrease  in  stiffness  is  associated  with  further 
extension  of  the  cracks  in  the  45  deg.  plies,  and  growth  of  delamination 
along  the  45  deg.  matrix  cracks.  Early  in  the  life  of  the  laminate  (region  I) 
the  damage  is  confined  to  the  approximate  zone  of  stress  concentration  around 
the  hole.  In  later  stages  of  life,  the  stresses  around  the  hole  are  redis¬ 
tributed  and  the  stress  concentration  around  the  hole  is  changed  due  to  the 
growth  of  the  damage  zone  away  from  the  hole.  Radiographs  of  damaged  type  A 
and  type  B  specimens  are  shown  Figure  9  on  page  11  and  Figure  10  on  page  11. 
The  data  in  Figure  4  are  presented  in  terms  of  percent  change  in  stiffness 
because  the  exact  values  of  stiffness  are  not  repeatable  in  situations  where 
the  cyclic  test  is  stopped,  the  extensometer  is  removed  from  the  specimen, 
and  the  specimen  is  removed  from  the  grips  for  nondestructive  evaluation  at 
various  cyclic  intervals.  However,  if  the  stiffness  values  are  determined  at 
the  beginning  and  end  of  each  interval,  or  are  recorded  throughout  a  cyclic 
interval  by  the  data  aquisition  system,  the  change  in  stiffness  during  the 
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cyclic  interval  can  be  determined.  7'he  total  change  is  then  computed  as  the 
sum  of  the  stiffness  changes  for  all  intervalls. 


3.5  Continuous  Stiffness  Recording 

The  plots  of  the  results  of  the  peak  detector  program  on  a  logarithmic  life¬ 
time  scale  confirm  the  estimated  stiffness  curves  based  on  quasi-static 
tests.  The  suggestion  of  a  five  percent  stiffness  loss  after  100  cycles  at  a 
90  percent  cycling  load  for  both  laminates  was  taken  ir*o  account  by  multi¬ 
plying  the  first  stiffness  value,  received  by  the  neak  detector  after  200 
cycles,  by  the  factor  0.95. 

Figure  5  shows  the  regions  II  and  III  of  the  stiffness  degradation  due  to 
fatigue  loading  at  the  90  percent  level  of  failure  load  of  each  specimen, 
respectively.  The  time  axis  is  logarithmically  scaled.  In  this  represen¬ 
tation,  region  II  appears  as  an  almost  linear  line  loading  into  a  sharp  kii'-r- 
which  connects  region  II  with  region  III.  Immediately  tee -id  the  knee,  lie- 
degradation  rate  appears  to  be  highly  increased.  At  t.h*»  end  of  region  III, 
the  stiffness  decay  accelerates  again  until  failure. 
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Figure  5.  Stiffness  Degradation  Versus  Logarithmic  Cycles 

Comparing  the  two  curves  shown  in  Figure  5,  we  find  the  B- laminate  showing  a 
lower  stiffness  degradation  before  the  knee,  a  later  and  more  distinct  form¬ 
ing  of  the  knee,  and  a  higher  degradation  rate  after  the  knee  than  the 
A- laminate.  The  stiffness  loss  at  the  end  of  fatigue  life  is  18  percent  for 
the  B-specimcn  and  40  percent  for  the  A-specimon. 

The  data  have  been  replotted  using  a  linear  scale  for  the  life  axis.  Figure  h 
on  page  8.  Throughout  the  lifetime,  the  A-sper.  imen  has  a  lower  stiffness: 
than  the  B -spec imen .  Between  10  K  and  50  K  cycles,  the  degradation  rates  of 
both  carves  are  similar;  after  that  the  A-specimen  looses  stiffness  more  rap- 
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idly  than  the  B-specimen.  Although  the  total  stiffness  loss  during  its  life¬ 
time  is  lower  for  the  B-specimen,  the  stiffness  degradation  rate  shortly 
before  failure  is  very  high.  Both  laminates  exhibit  three  sudden  stiffness 
changes  during  the  fatigue  life:  the  first  one  occurs  during  the  very  first 
cycles  of  loading;  the  second  one  occurs  within  the  first  half  of  the  life¬ 
time  for  the  A-specimen  and  well  before  failure  of  the  B-specimen,  and  the 
third,  sudden  stiffness  lossed  leads  into  the  failure  of  the  specimens.  The 
A-specimen  has  its  second  sudden  stiffness  loss  at  70  K  cycles  and  the 
B-specimen  at  157  K  cycles.  Both  events  coincide  with  the  onset  of  the  knee 
in  the  logarithmic-scale  stiffness  representation  and  are  generally  followed 
by  an  increased  degradation  rate.  It  should  be  mentioned  here  that  locally 
confined  buckling  effects  were  seen  during  cyclic  loading  of  the  A-specimen. 


Figure  6.  Stiffness  Degradation  Versus  Linear  Cycles 

On  both  sides  of  the  specimen,  the  zero-deg  and  the  90-deg  plies,  as  a  sub¬ 
laminate,  delaminated  from  the  inner  {45,-45|s  block  at  locations  governed  by 
compressive  circumferential  stresses  and  buckled  outwards.  It  could  be 
hypothesized  that  the  forming  of  these  delaminations  and  in  company  the 
stiffness  degradation  is  being  rapidly  accelerated  after  reaching  a  critical 
delamination  size  which  facilitates  buckling.  The  sudden  stiffness  change  at 
70  K  for  the  A-specimen  could  then  be  related  to  the  observed  buckling 
effects.  By  eye  inspection,  no  buckling  effects  were  visible  for  the 
B-specimen.  However,  at  the  end  of  the  fatigue  life,  it  seemed  to  be  obvious 
by  eye  observation  that  the  surface  45-deg  ply  contributed  only  little  or 
nothing  to  the  load  transfer,  since  around  the  hole  matrix  cracks  were  formed 
and  some  of  the  surface  material  confined  by  those  cracks  was  almost  com¬ 
pletely  delaminated.  The  matrix  cracks  prevent  the  in-plane- load  transfer 
and  the  delaminations  make  load  transfer  through  interlaminar  shear  stresses 
impossible.  These  effects  were  visible  at  the  end  of  the  lifetime  of  the 
B-specimen,  thus  suggesting  the  possibility  of  a  relationship  between  them 
and  the  sudden  stiffness  change  at  157  K  cycles.  Finally  it  should  be 
recalled  that  the  90  percent  load  levels  for  the  A-  and  the  B-specimens  are 
represented  by  the  absolute  values  of  2354  lbs  and  2195  lbs,  respectively. 


3.6 


Residual  Strength 


A  number  of  the  cyclic  tests  were  halted  at  various  stages  of  the  loading 
history  for  nondestructive  inspection  of  the  damage  zone  and  residual 
strength  measurements.  This  series  of  tests  provided  a  characteristic  and 
repeatable  curve  of  the  change  in  residual  strength  throughout  the  loading 
history . 

During  the  first  few  (ten)  cycles,  the  damage,  primarily  matrix  cracks, 
reduces  the  effect  of  the  stress  concentration  at  the  hole  and  the  residual 
strength  increase  is  on  the  order  of  ten  percent  The  actual  strength 
increase  does  depend  on  the  cyclic  stress  level,  with  the  higher  stress  lev¬ 
els  corresponding  to  slightly  higher  strength  increases  than  the  lower  stress 
leve  1  s . 

Stiffness  change  and  residual  strength  data  for  a  type  R  laminate  after  100 
cycles  at  a  low  cyclic  stress  amplitude  are  also  shown  in  rigure  4.  Both  the 
increase  in  residual  strength  and  the  decrease  in  stiffness  are  less  than 
those  for  the  same  type  specimen  cycled  at  the  85  percent  stress  level.  The 
second  region  of  the  residual  strength  curve,  approximately  ten  percent  of 
the  fatigue  life,  is  characterized  by  a  slow,  but  constant  increase  in  resi¬ 
dual  strength  over  a  logarithmic  cycles  scale.  At  the  end  of  the  second 
region,  the  increase  in  residual  strength  is  115  to  120  percent.  The  limits 
of  regions  I  and  II  for  residual  strength  are  approximately  the  same  as 
observed  for  stiffness  change. 

The  third  region  of  the  residual  strength  curve  is  marked  by  a  further 
increase  in  residual  strength.  Although  the  data  in  this  region  is 
incomplete,  it  appears  that  the  maximum  residual  strength  is  reached  between 
50  and  80  percent  of  the  lifetime.  The  largest  increases  in  residual  strength 
are  26  and  42  percent  for  the  type  A  and  B  laminates,  respectively.  During 
the  final  region  of  the  fatigue  life,  the  residua)  strength  decreases  until 
the  strength  equals  the  level  of  maximum  cyclic  stress,  as  shown  in  Figure  4 
and  F igure  7  . 


Figure  7.  Residual  Strength  of  an  A-Specimen 
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The  last  two  regions  of  the  residual  strength-log  cycles  curve  correspond  to 
the  third  region  of  the  stiffness  reduction  curve  where  the  stiffness  change 
shows  a  sharp  decrease  on  the  log  cycles  plot.  The  logarithmic  scale  used  to 
represent  the  lifetime  may  cause  some  confusion  in  interpreting  these 
results.  The  residual  strength  data  have  been  replotted  in  Figure  8  using  a 
linear  scale  to  represent  cycles. 


Figure  8.  Residual  Strength  Versus  Linear  Cycles 


3.7  Damage  Patterns 

Static  damage  starts  with  matrix  cracks  usually  in  the  surface  plies  at  1570 
lbs  load  for  the  A-laminate  and  '220  lbs  load  for  the  B-laminate.  At  higher 
load  levels,  cracks  appear  in  all  directions.  The  A-laminate  develops  more 
well  defined  zero-deg  cracks  tangent  to  the  hole  while  the  B-laminate  devel¬ 
ops  more  well  defined  45-deg  cracks.  Note  that  these  are  the  directions  of 
the  surface  layers,  respectively.  In  both  laminates,  A5-deg  and  zero-deg 
cracks  already  exist.  Zero-deg  cracks  in  the  A-laminate  grow  to  a  length  of 
more  than  five  hole  diameters.  When  the  tangent  zero-deg  cracks  approach 
their  maximum  length,  cracks  in  95-deg  and  90-deg  direction  are  also  initi¬ 
ated  on  the  straight  edge.  In  the  final  stages  of  lifetime,  a  regular  pattern 
of  matrix  cracks  is  formed  between  the  tangent  zero-deg  cracks  and  the 
straight  edges.  The  uniformity  of  these  cracks  is  similar  to  damage  patterns 
of  unnotched  laminates,  so-called  characteristic  damage  states  (CDS). 

The  delamination  zone  visible  on  an  X-ray  photo.  Figure  9  and  Figure  10,  con¬ 
sists  of  the  superimposed  images  of  the  de laminat ions  on  the  various  laminate 
interfaces.  It  seems  as  if  the  tangent  zero  deg  cracks  in  both  the  A-  and  the 
R-laminates  cause  de 1 amina t ions  on  adjacent  interfaces.  Those  de laminations 
have  the  largest  continuously  connected  area  of  all  delaminations  in  the  1am- 
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inates .  In  general,  while  the  A-laminate  develops  a  much  larger  delamination 
along  its  tangent  zero-deg  cracks,  the  more  confined  delaminations  of  the 
B-laminate  spread  along  the  45-deg  direction  and  reach  the  straight  edge.  It 
is  emphasized  that  Figure  11  or  page  12  gives  the  spatial  distributions  and 
the  characteristic  shapes  as  well  as  the  proper  magnitude  of  the  delamination 
zones.  Noting  that  the  area  of  delaminations  on  all  interfaces  is  higher  for 
the  A-laminate  than  for  the  B-laminace,  we  recall  that  the  A-specimen  experi¬ 
enced  much  higher  stiffness  losses  during  their  fatigue  life  than  the 
B-specimens  . 
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Figure  11.  Delamination  Areas  on  Interfaces 


3.8  Damage  and  Mechanical  Response 

Figure  12  on  page  13  shows  that  the  tangent  zero-deg  cracks  have  a  strong 
redistributing  effect  on  the  strain  field  near  the  hole,  reducing  the  maximum 
strain  concentration  in  x-direction.  Except  at  the  crack  tips,  the  strain 
distribution  along  the  cracks  is  almost  uniform. 

Figure  13  on  page  14  presents  the  average  strain  concentrations  as  a  function 
of  the  length  of  the  tangent  cracks  at  various  fatigue  stages  for  the 
A-specimen.  After  600K  cycles,  which  is  estimated  as  being  50  percent  of  the 
expected  lifetime  the  strain  concentration  is  reduced  to  K=1.46. 

Figure  14  on  page  14  shows  the  decrea;e  of  strain  concentration  versus  the 
logarithmically  scaled  number  of  cycles.  The  rate  of  decrease  is  very  high 
during  the  first  cycles  of  fatigue  loading. 

Recalling  that  the  rate  of  increase  of  residual  strength  of  both  laminates  is 
also  very  high  at  the  beginning  of  fatigue  life,  it  seems  very  likely  that 
the  tangent  zero-deg  cracks  are  a  main  factor  influencing  residual  strength. 


Figure  12.  Undamaged  Versus  Low  Damage  Situation 

Figure  15  on  page  15  snows  the  results  of  two  strength  tests  performed  on  an 
A-specimen  and  a  B-specimen  versus  the  length  of  the  tangential  cracks 
obtained  from  X-ray  photos  taken  after  applying  the  fatigue  load. 

We  note  that  cracks  of  equal  length  in  A-  and  B-specimen  seem  to  influence 
the  strength  increases  of  A-  and  B-specimens  differently.  As  a  result  of  a 
linear  curve  fit  for  the  data  points  we  find  an  increase-in-strength  rate  of 
44  percent  per  inch  for  the  B-  and  23  percent  per  inch  for  the  A-laminate.  A 
residual-strength  versus  s tra  in-concent  rat  ion  curve  (Figure  16  on  page  15) 
was  constructed  from  the  data  used  in  Figure  14  on  page  14  and  Figure  15  on 
page  15  by  eliminating  the  crack  length. 

Figure  17  on  page  16  and  Figure  18  on  page  17  show  the  increase  of  cracks  in 
the  zero-deg  plies  with  the  number  of  cycles  under  fatigue  load.  The  crack 
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Strain  Concentration  Versus  Length  of  Tangent  2ero-Deg  C; 


Figure  14.  Strain  Concentration 
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lengths  are  the  average  values  of  the  length  of  the  two  cracks  visible  on 
X-ray  photos.  For  two  reasons  it  is  difficult,  however,  to  measure  average 
crack  lengths  exactly.  First,  the  length  of  the  dark  line  on  a  X-ray  matches 
the  actual  crack  length  only  in  the  case  of  complete  penetration  with  zinc 
iodide.  Second,  it  is  difficult  to  decide  what  the  length  of  a  crack  is  since 
a  dark  line  on  a  X-ray  is  the  superimposed  image  of  actually  two  cracks  on 
each  side  of  the  laminate  midplane.  If  these  two  actual  cracks  are  shifted  in 
different  directions,  one  might  observe  a  longer  crack  than  actually  exists. 
Since  the  data  for  the  B-laminate  is  more  complete,  Figure  18  provides  also 
information  on  the  influence  of  the  load  parameter  on  the  growth  rate  of 
those  cracks.  Figure  17  takes  the  initial  crack  length  due  to  static  load 
into  account,  so  the  crack  length  starts  with  the  value  indicated. 


Figure  17.  Zero-Deg-Crack  Growth  of  an  A-Specimcn 

The  phenomenon  of  widely  scattered  fracture  surfaces  for  specimens  with  high 
loading  histories  can  be  related  to  the  shift  of  high  strain  concentrations 
away  from  the  hole.  Since  those  strain  concentrations  are  presumably  highest 
at  the  tips  of  the  cracks,  we  conclude  that  catastrophic  failure  can  be  ini¬ 
tiated  right  there.  Many  zero-deg  fibers  broke  at  a  distance  from  the  hole 
of  6  to  7  hole  diameters  or  from  2.25  inches  to  2.2625  inches.  As  for  the 
phenomenon  of  stiffness  degradation,  we  note  that  the  A- laminate  exhibits 
much  larger  de laminat ions  than  the  B-laminate.  These  de 1  am inat ions  and  also 
the  buckling  described  in  the  section  on  continuous  stiffness  recording  con¬ 
tribute  to  higher  stiffness  losses  in  the  A-laminate  than  in  the*  R- laminate. 
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Figure?  18.  Zoro-Deg-Crack  Growth  of  B-Specimons 
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1.  Introduction 


In  future  spacecraft  application  satellite  structures  in  the 
form  of  antennas  or  platforms  will  become  increasingly 
important.  Because  of  their  attractive  structural  properties 
graphite/ epoxy  or  graphite/polymide  composites  are  the  prime 
material  candidates  for  such  applications.  Besides  their 
specific  strength  and  stiffness  properties  they  offer  low 
coefficients  of  thermal  expansion.  So  that  deformations  of  the 
structure  under  temperature  exposure  can  be  held  to  a  minimum. 
However,  on  the  microscopic  level,  the  large  differences  in 
thermal  expansion  coefficients  between  the  resins  and  the  fibers 
lead  to  thermally  induced  stresses  of  substantial  magnitude  when 
the  material  undergoes  thermal  cycling  as  a  consequence  of  sun 
and  shadow  phases  in  an  earth  orbit.  The  corresponding  upper  and 
lower  temperature  limits  depend  on  the  surface  coating  of  the 
structure  and  on  the  attitude  angle  relative  to  the  sun.  They 
may  range  from  -160°C  in  deep  shadow  to  +100°C  and  even  more 
when  exposed  to  the  sun.  Depending  on  orbit  height,  inclination 
and  service  time  3000  to  5ooo  thermal  cycles  may  occur  during 
the  life  of  a  typical  vehicle  and  thereby  introduce  thermal 
fatigue  problems  in  the  material.  Additionally,  e lec t ronmagne t i c 
and  particle  irradiation  may  degrade  the  properties  of  the 
organic  matrix. 

In  order  to  assess  the  effects  of  severe  thermal  cycling  and  of 
electron  radiation  on  the  mechanical  properties  of  various 
carbon  fiber  composites  a  test  program  was  conducted  with 
emphasis  on  a  comparison  of  their  initial  and  residual  material 
prope  rt ies . 

2.  Thermal  Properties  of  Carbon-Fiber  Reinforced  Plastics  (CFRP) 

The  thermal  expansion  coefficient  of  graphite  fibers  in  the 
axial  direction  is  nearly  zero.  It  decreases  as  the  modulus 
increases  to  a  point  where  high  modulus  fibers  exhibit  negative 
axial  thermal  expansion  coefficients  [  1  ].  Epoxy  resins, 
however,  have  thermal  expansion  coefficients  in  the  order  of 
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2  times  that  of  aluminum.  The  thermal  expansion  of  cured 
laminates  depends  both  on  the  material  choice  and  the  stacking 
order.  Table  1  gives  thermal  expansion  coefficients  of  several 
laminates  prepared  from  different  materials.  All  measurements 
were  performed  over  a  60K  temperature  range  (30°C  to  90°C)  using 
a  quartz-tube  dilatometer  having  an  accuracy  of  about  10  m.  The 
specimens  were  preconditioned  at  70°C  for  14  days  in  dry  air  to 
minimize  moisture  effects.  The  values  shown  in  Table  1  are 
average  values  of  4  or  5  specimens.  The  thermal  expansion 
coefficients  of  all  unidirectional  laminates  are  negative  in  the 
fiber  direction.  In  accordance  with  ref.  [  1  ]  the  contraction 
is  greater  in  laminates  with  high  modulus  fibers  than  in  those 
with  high  tension  fibers.  Perpendicular  to  the  fiber  direction 
the  thermal  expansion  coefficients  of  the  laminates  which 
contain  60%  fibers  by  volume  correspond  to  that  of  aluminum. 

The  data  in  Table  1  supplement  the  publicly  available 
information  on  expansion  coefficients  of  CFRP  in  the  above 
mentioned  temperature  range.  Extrapolations  beyond  this  range 
especially  into  very  low  temperature  regimes  encountered  in 
space,  are  hardly  possible  because  thermal  expansions  do  not 
depend  linearly  on  temperature.  Thermally  induced  stresses  are 
initated  in  a  laminate  during  the  curing  process  when  it  reaches 
a  temperature  at  which  the  resin  changes  from  a  rubbery  gel  to  a 
glassy  solid  and  the  cross-linking  reactions  slow  down.  At  this 
temperature  no  stress  exists  between  the  fibers  and  the  matrix. 
It  is  slightly  below  the  final  cure  temperature  of  the  resin  and 
is  called  the  lock-on  temperature. 

When  the  laminate  is  cooled  down  to  ambient  temperature  the 
different  thermal  expansion  coefficients  of  fibers  and  matrix 
generate  tensile  stresses  in  the  matrix  and  compressive  stresses 
in  the  fibers  along  the  fiber  axes  in  each  lamina.  In  addition 
to  longitudinal  stresses  depending  on  the  fiber  volume  fraction 
the  transverse  thermal  expansion  coefficient  of  the  fibers  and 
the  mismatch  of  the  Poisson's  ratios  of  the  fibers  and  matrix, 
radial  and  circumferential  stresses  can  be  generated.  [  2,3  ] 
Beyond  this  in  cross/plied  laminates  shear  stresses  and  tensile 
stresses  transverse  the  fiber  directions  will  he  generated  by 
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the  constraint  of  the  different  fiber  directions.  Because  of  the 
dependence  of  the  stress  state  on  the  difference  between  curing 
temperature  (respectively  lock-on  temperature)  and  ambient 
temperature  it  may  be  concluded  that  resins  requiring  low  curing 
temperatures  should  be  favoured  for  space  applications.  But 
nonlinearity  of  the  thermomechanical  properties  of  the  laminates 
and  /  or  the  admissible  strain  of  the  matrix  without  crack 
formation  may  repress  the  influence  of  the  curing  temperature. 

3 .  Material  Selection 

For  the  evaluation  of  the  thermal  properties  of  CFRP  as  discussed 
in  chapter  2  different  materials  were  selected  including  systems 
with  high  tension  fibers  as  well  as  high  modulus  fibers.  The 
curing  temperatures  are  ranging  from  120°C  to  210°C.  All  of  them 
were  epoxy  systems  except  for  one  polymide  system  identified  as 
PJ.  A  survey  of  these  materials  is  given  in  Table  2.  The  volume 
fraction  of  fibers  and  resin,  respectively,  was  equal  in  all 
Materials . 

4 .  Laminate  Fabrication 

Of  each  material  at  least  one  laminate  was  prepared  with  the 
size  380x380mm  and  numbered  corresponding  to  Table  2.  The  curing 
cycles  used  in  fabrication  process  met  the  instructions  of  the 
resin  producers  in  all  points  with  the  exception  that  a 
pressclave  was  used  instead  of  an  autoclave. 

The  only  polyimide  system,  laminate  P J ,  was  prepared  by  Dornier 
System,  Friedrichshaf en .  The  laminates  164,  168,  169,  175,  246, 
and  247  were  manufactured  from  industry  supplied  unidirectional 
prepregs  in  the  own  laboratory.  The  prepreg  of  laminate  170A  was 
self /prepared  by  impregnating  the  carbon  fibers  and  winding  them 
on  a  drum. 

The  laminates  consisted  of  8  plies  with  the  stacking  sequence 
[2(*45°,  -45°)]S,  only  laminate  247  was  composed  of  4  plies 
because  of  its  double  prepreg  thickness.  The  total  thickness  of 
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all  laminates  was  1mm.  The  mechanical  properties  of  laminates 
with  the  selected  stacking  order  are  highly  matrix  controlled. 
Therefore,  damages  of  the  matrix  as  matrix  cracks  or  delinations 
•aill  substantially  reduce  the  residual  strength  making  the 
discovery  of  such  damages  possible  by  simple  strength  tests. 

Four  specimens  were  cut  from  the  laminates  164,  168,  169.  175 
and  PJ  respectively.  One  was  for  reference  purposes  while  the 
others  had  to  sustain  different  numbers  of  thermal  cycles.  The 
number  of  cycles  readied  by  one  specimen  of  each  material 
respectively  were  1170,  2295,  and  3480. 

5 .  St ress  Calculation  in  T2(*450,  -45°)  \S  Laminates 

The  calculation  of  thermally  induced  stresses  in  a  +  45°,  -45° 
symmetrical  8- layer  laminate  was  conducted  on  the  basis  of  the 
classical  lamination  theory.  The  lock-on  temperature  was  chosen 
to  be  180°C,  while  thermal  expansion  coefficients  of  fibers  and 
matrix  and  the  elastic  properties  corresponded  to  those  of  914 
C-TS-5.  The  laminate  was  considered  to  extend  infinitely  in  both 
directions.  The  calculation  provided  the  compressive  and  tensile 
stresses  in  the  individual  layer  parallel  and  normal  to  the 
fiber  direction  respectively  as  well  as  the  shear  stresses  when 
the  laminate  is  cooled  down  to  -155°C. 

For  each  layer  the  calculation  indicated  a  tensile  stress  in  the 
matrix  perpendicular  to  the  fibers  of  at  least  71  MPa  while  the 
breaking  strength  of  the  resin  is  only  65  MPa.  Severe  cracking 
parallel  to  the  fibers  should  be  the  consequence  when  exposed  to 
very  low  temperatures  encountered  in  space  provided  that  the 
t hermo- e las t ic  properties  of  the  material  is  linear.  The  shear 
stresses  are  of  the  same  magnitude  as  the  tensile  stresses,  i.e. 
71  MPa.  This  is  well  below  the  92  MPa  shear  strength  of  914  C. 
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6.  Thermal _ Cycling  Tests 

The  thermal  cycling  tests  were  conducted  in  a  vacuum  test 
facility  of  the  DFVLR.  Inside  a  vacuum  chamber  of  approximatly 
100cm  diameter  a  rotatable  cage  was  located.  The  specimens  were 
mounted  on  support  provisions  attached  to  the  outer  surfaces  of 
two  opposite  sides  of  the  cage  (Fig.  l).Bv  automatic  rotation  in 
90°  intervals  the  specimens  were  exposed,  alternately,  to 
infrared  heaters  and  LM-.  -  cooled  plates  simulating  solar 
heating  or  energy  dissipation  by  deep  space  effects  (Fig.  2). 

The  rear  sides  of  the  specimens  were  facing  the  walls  of  the 
cage  whose  inside  was  permanently  cooled  by  LN2  .  In  this  fashion 
a  test  specimen  temperature  of  95°C  was  realized  during  hot 
phase,  and  of  -155°C  during  cold  phase.  The  upper  and  lower 
temperature  limits  were  controlled  by  the  exposure  time  of  9  min 
during  the  hot  phase  and  of  41  min  during  the  cold  phase. 

—  4 

!••  ;u It ing  in  a  cycle  time  of  50  min.  A  chamber  pressure  of  10  Pa 
was  maintained  throughout  the  test.  A  typical  thermal  cycle  is 
shown  in  Fig.  3. 

Because  of  the  viscoelastic  behavior  of  the  specimens  matrix 
realistic  thermal  stresses  will  only  be  obtained  when  the 
simulation  of  the  temperature- time  gradients  are  as  realistic  as 
possible.  This  was  accomplished  in  a  nearily  perfect  fashion  by 
the  test  arrangement. 

7 .  El e c t r o n _ Irradiation  Test s 

Space  structures,  in  general,  are  subject  to  electron  and  proton 
irradiation.  The  effects  of  electron  irradiation  on  the 
laminates  168,  246,  247,  and  PJ  (Table  2)  were  investigated  by 

exposure  of  several  specimens  to  electron  beam  irradiation  of 
6  **  ** 

3x10  rads,  at  100°C,  and  in  10  Pa  vacuum.  The  beam  current  was 
held  at  165nA  for  the  required  run  time  of  twenty/seven  (27) 
hours.  The  test  was  conducted  in  a  Van  der  Graaf  accelerator  at 
the  University  of  Illinois,  USA. 
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8.  Evaluation  of  Environmental  Effects 


The  extent  of  damage  caused  in  the  laminates  by  environmental 
conditions  was  investigated  by  means  of 

-  Ultrasonic  inspections  to  detect  de laminat ions  generated  by 
shear  stresses 

-  Radiographs  for  the  detection  of  matrix  cracks  due  to  tensile 
stresses  perpendicular  to  the  direction  of  the  fibers 

-  Electron-beam  micrographs  to  detect  microcracks  and  debonding 
of  fibers  and  matrix 

-  Measurements  of  the  natural  frequency  and  the  logarithmic 
decrement  of  damping  in  transverse  oscillation 

-  Measurements  of  the  spectral  reflection  properties  of  the 
laminate  surface 

-  Comparison  between  initial  and  residual  strength  and 
stiffness  properties. 

All  of  the  test  methodes  were  applied  before  and  after 
increasing  numbers  of  thermal  cycles.  Strength  and  stiffness 
tests  were  established  both  at  room  temperature  and  at  100°C.  In 
order  to  obtain  reproducible  results  the  specimens  were  dried 
prior  to  testing.  With  respect  to  the  irradiated  specimens  the 
tests  were  performed  prior  and  after  completion  of  the  electron 
beam  irradiation. 

8 . 1  Ultrasonic  Inspections 

In  Fig.  A  to  7  amplitude  scans  are  shown  of  the  laminates 
164,  169,  170A,  and  175  before  and  after  1170  thermal  cycles. 

The  laminate  numbers  and  the  corresponding  material  designations 
are  contained  in  Table  2.  Within  the  accuracy  limits  of 
ultrasonic  techniques,  the  laminates  164,  169,  and  175  are  free 
of  defects.  The  amplitude  scans  of  laminate  170A  show  minor 
disturbances  probably  caused  by  pores  progressing  in  the  fiber 
directions  as  it  was  detected  by  scanning  electron  beam 
investigations.  In  contrast  to  the  other  laminates  effects 
of  thermal  cycling,  if  they  occured  cannot  be  found  by 
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comparison  of  the  two  scans  of  Fig.  6  because  prior  to  cycling 
ultrasonic  inspection  was  only  applied  to  the  reference 


specimen . 

Laminate  168  was  prepared  with  artificially  introduced 
delaminations  at  the  center  of  each  specimen  for  investigation  of 
shear  effects  on  the  growth  of  such  a  damage.  As  evidenced  by 
Fig.  8  no  growth  due  to  thermally  induced  stresses  was 
experienced  after  3480  thermal  cycles. 

In  contrast  to  the  epoxy  based  laminates,  laminate  PJ  exhibited 
a  rather  smooth  amplitude  scan  prior  to  cycling  but,  according 
to  Fig.  9,  indicated  pores  and  extend  de laminat ions  after  1170 
thermal  cycles.  The  higher  curing  temperature  of  this  material 
raises  the  thermally  induced  stresses  by  6  to  77,  over  those  of 
the  former  discussed  graphite/  epoxy  laminates  and  might  have 
contributed  to  these  damages. 

8.2  R_a di o g rapjis 

Without  special  provision  radiographs  of  thin  CFRP  laminates 
exhibit  poor  contrast.  In  order  to  raise  the  quality  of  the 
radiographs  high-density  fluids  like  Tetrabromethan  (TBE )  may  be 
used.  After  penetration  of  the  contrasting  agent  into  cracks  and 
pores  in  a  TBE-bath,  those  defects  become  visible  because  the 
fluid  absorbs  X-rays  to  a  higher  degree  than  the  sorrounding 
CFRP. 

At  the  beginning  of  the  program  no  experience  existed  wether 
this  technique  would  affect  the  structural  integrity  of  CFRP. 
Therefore,  X-ray  photographs  of  untreated  epoxy/based  laminates 
were  taken  before  thermal  cycling.  During  the  thermal  cycling 
test  detailed  investigations  showed  that  the  fluid  can  be 
removed  by  vacuum  and  heat  without  detectable  defects  in  ♦  45°, 

-  45°  laminates.  Therefore,  all  radiographs  taken  after  thermal 
cycling  were  cont ras t - 1 if t ed  by  TBE. 

The  photographs  in  Fig.  10  show  no  or  little  cracking  in  the 
TBE-treated  laminates  164,  169,  and  175  after  cycling.  The 
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photograph  of  laminate  170A  in  the  right  of  Fig.  11  indicates 
substantial  damage,  however.  From  the  matching  pattern  of  the 
striations  in  the  photograph  on  the  left  it  appers  that  damage 
existed  already  before  the  thermal  cycling  commenced.  The 
reversal  of  the  contrast  is  a  consequence  of  the  fact  that  for 
this  photograph  no  TBE  was  used.  Subsequent  investigations  with 
an  electron  beam  microscope  demonstrated  that  the  damage  existed 
in  the  form  of  extended  internal  pores  in  the  fiber  direction  of 
the  individual  plies  rather  than  of  microcracks.  Evidently,  the 
majority  of  the  damage  was  caused  by  faulty  fabrication. 

Severe  matrix  cracking  as  predicted  by  lamination  theory  could 
not  be  observed  in  any  of  the  epoxy-based  laminates.  This 
indicates  that  the  use  of  linear  constitutive  laws  leads  to 
unrealistic  stress  states. 

The  two  radiographs  in  Fig.  12  of  a  polyimide  based  laminate  PJ 
were  also  taken  before  and  after  thermal  cycling.  While  the 
laminate  seems  to  be  initially  free  of  any  defects,  a  multitude 
of  thin  dark  lines  are  visible  after  thermal  cycling. 
Subsequently  taken  electron  beam  micrographs  revealed  that  these 
lines  do  not  represent  internal  pores  as  observed  previously, 
but  rather  microcracks  emanating  from  the  surface  of  the 
specimen  (Fig.  18). 

As  was  expected  no  damage  was  found  by  ultrasonic  inspections 
due  to  electron  irradiation. 

8 . 3  Scanning  Electron  Beam  Micrographs 

The  effects  of  thermal  cycling  on  all  of  the  material  systems 
were  investigated  also  by  scanning  electron  beam  microscopy.  In 
several  specimens  prepared  from  laminate  164  no  defects  were 
detected  prior  to  cycling.  After  thermal  cycling  only  a  few 
short  cracks  extending  through  the  thickness  of  one  ply  were 
visible.  Its  discovery  was  possible  only  after  removal  of  the 
dust  particles  covering  the  polished  surfaces  by  ionic  etching. 
As  can  be  seen  in  Fig.  13  the  cracks  were  not  confined  to  the 
matrix  material  but  sometimes  affected  the  fibers  as  well.  This 


observation  was  made  in  919C-TS-5  and  T3T  F178  laminates.  It  must 
be  recognized,  however,  that  the  irregular  distribution  of  such 
defects  and  the  limited  field  of  vision  of  the  microscope  make 
it  difficult  to  establish  reliable  correlations . 


Fig.  19  and  15  show  accumulations  of  voids  along  the  ply 
interface  in  samples  cut  from  laminates  169  and  170A.  No 
significant  differences  could  be  detected  in  the  charakter  of 
defects  observed  before  and  after  thermal  cycling.  However,  in 
laminate  169  the  voids  are  limited  to  small  local  areas  whereas 
in  laminate  170A  the  voids  seem  to  run  together  and  extend  over 
the  whole  laminate.  Sometimes  even  indicating  the  formation  of 
delaminations.  This  observation  is  consistent  with  the 
evaluation  of  the  ultrasonic  records  in  Fig.  6  and  X-ray 
photographs  in  Fig.  11. 

Samples  from  laminate  175  (Fig.  16)  show  similarity  to  those  of 
laminate  169  except  that  more  plies  seem  to  be  affected  by  crack 
formation.  Again  only  a  few  cracks  were  created  by  thermal 
cycling. 


Fig.  17  show  a  cross  section  and  a  surface  view  of  the  T3T  F178 
laminate  PJ .  They  appear  smooth  and  free  of  defects  before 
thermal  cycling,  corroborating  the  previously  obtained 
ultrasonic  C-scans  and  radiographs.  The  photographs  in  Fig.  18 
taken  after  1170  thermal  cycles  exhibit  comparatively  large 
cracks  forming  a  regular  pattern  on  the  surface  (top  photo). 
After  progression  of  the  cracks  through  the  first  ply  a  tendency 
toward  formation  of  small  de laminat ions  along  the  interface  with 
the  adjecent  ply  could  be  observed  (photo  at  bottom). 
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of  the  interface  between  adjecent  layers  which  lead  to  final 
rupture . 

rig.  20  shows  the  same  secimen  area  as  in  Fig.  19  but  taken  at 
a  higher  magnification.  The  photos  exhibit  a  distinct  weakening 
of  the  bonding  strength  between  fibers  and  matrix  due  to  thermal 
cycling.  When  strength  tested  before  thermal  cycling  the 
fracture  surfaces  of  the  individual  plies  were  rough  with  some 
broken  fibers  visible.  After  thermal  cycling  they  showed  smooth 
traces  of  fiber  pull  out. 

In  the  case  of  the  laminate  175  the  failure  mechanism  was 
different  from  that  of  the  laminate  164  before  thermal  cycling. 
Lower  strength  of  the  pitch  fibers  and  a  high  bonding  strength 
between  fibers  and  matrix  tended  to  create  ruptures  transvers 
all  plies  where  fracture  of  fibers  was  a  normal  mode  of  failure. 
After  thermal  cycling  the  failure  mechanism  was  similar  to  that 
of  laminate  164.  As  Fig.  22  indicates  no  significant  change  in 
bonding  between  fibers  and  matrix  took  place.  No  change  of 
appearance  was  defected  after  electron  irradiation. 

8 . 4  Natural  Frequencies  and  Decrements  of  Damping 

Table  3  shows  experimentally  determinded  natural  frequencies  and 
decrements  of  damping  of  the  candidate  materials  after  different 
numbers  of  thermal  cycles.  In  all  cases  the  natural  frequency 
decreases  with  increasing  numbers  of  thermal  cycles.  However,  as 
shown  in  Fig.  23  and  24,  the  major  change  took  place  soon  after 
the  beginning  of  cycling,  probably  during  the  first  few  cycles. 
This  may  correspond  to  the  development  of  initial  cracks  which 
ought  to  be  expected  during  the  same  time  frame.  Because  of  the 
stress  relaxation  further  cracking  and  crack  extension  will  be 
slower  with  increasing  numbers  of  thermal  cycles. 

The  decrease  of  the  natural  frequency  of  the  poly imide - bases 
laminate  was  three  or  more  times  larger  than  tnat  of  the  other 
materials  as  is  apparent  in  Fig.  24.  This  may  be  a  consequence 
of  the  severe  cracking  observed  in  these  samples.  As  indicated 
by  Fig.  23  and  24  the  decrement  of  damping  of  laminates  164, 
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169,  and  170A  increases  with  increasing  numbers  of  thermal 
cycles  indicating  growing  internal  friction.  The  continous 
increase  of  friction  is  not  only  caused  by  matrix  cracking  since 
the  formation  of  cracks,  as  previously  described,  tends  to 
diminish  because  of  stress  relaxation.  It  is  probable  that  the 
molecular  structure  of  the  resins  is  affected  by  fatigue,  or  the 
interface  properties  between  fibers  and  matrix. 

For  the  contradictory  behavior  of  laminates  175  and  Pj  no 
physically  reasonable  explantion  can  be  offered  at  this  time. 

With  respect  to  electron  irradiation,  only  minor  effects  on 
frequencies  and  damping  characteristics  were  observed.  However, 
these  departures  from  the  non- irradiated  specimens,  shown  in 
Table  4,  are  so  small  that  they  may  be  due  to  the  limited 
accuracy  of  the  test  equipment.  As  only  one  sample  of  each 
laminate  was  available  a  confirmation  of  these  results  by 
statistical  means  was  not  possible  so  their  physical 
significance  is  questinable. 

Table  4:  Change  of  Natural  Frequencies  and  Damping  Decrements 
due  to  Electron  Irradiation 


Laminate  No. 

Material 

Frequency 

f  % 

Damping 

% 

168 

914  C  -  TS  - 

5 

-- 

-6,6 

246 

HYE  -  2034 

D 

*2,4 

-6,7 

247 

T6T  262  -  12  F 

550 

♦1,0 

-3,7 

PJ 

T3T  F  178 

-1,9 

_ ^5 

8 . 5  Relative  Spectral  Reflection 

Both  the  thermal  cycling  and  the  electron  irradiation  caused 
brownish  discolorations  of  the  specimen  surfaces.  Therefore, 
measurements  were  made  of  the  relative  spectral  reflection  of 
environmentally  tested  and  of  virginal  laminates.  Fig.  25  shows 
typical  results  of  these  investigations. 
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All  the  thermally  cycled  laminates  experienced  a  decrease  of 
relative  spectral  reflection.  The  formation  of  surface  cracks, 
as  observed  by  electron  beam  microscopy,  can  only  be  a  partial 
explanation  for  this  phenomenon  because  it  also  occured  in  two 
laminates  which  were  exposed  to  electron  irradiation  without 
experiencing  crack  formation  (Fig.  26).  On  the  other  hand, 
simulated  sun  irradiation  for  150  hours  with  an  intensity 


comparable  to  that  of  a  low  earth  orbit  caused  a  substantial 
increase  of  spectral  reflection  (Fig.  25),  which  may  be 
traceable  to  an  increasing  interlinking  of  the  molecular  chains 
of  the  matrix.  This  is  a  normal  event  due  to  energy  transmission 
to  the  laminate.  During  longer  exposure  times  the  high  energy 
photons  of  the  ultra  violet  region  of  sun's  light  spectrum 
damage  the  molecular  chains  consequently  leading  to  a  blunt 
laminate  surface.  Very  similar  effects  must  be  expected  as  a 
result  of  dense  fluxes  of  high  energy  electrons. 

Furthermore  the  fatigue  stresses  caused  by  thermal  cycling  might 
also  have  affected  the  molecular  structure  of  the  matrix  as  it 
was  supposed  in  chapter  8.4.  This  is  strongly  supported  by 
indications  that  the  decrease  of  spectral  reflection 
corresponded  continously  to  the  increase  of  the  number  of 
thermal  cycles  while  the  formation  of  cracks  which  are 
observable  by  micrography  slowed  down  as  discussed  earlier. 

8 . 6  Residual  Strength 

Tables  5  and  6  present  mean  values  and  their  standard  deviations 
of  the  static  strength  of  test  specimens  obtained  after  various 
numbers  of  thermal  cycles.  For  the  determination  of  the  orginal 
strength  of  the  laminates  six  specimens  were  tested,  whereas 
only  3  specimens  were  available  for  the  residual  strength  tests 
after  1170,  2295,  and  3480  cycles  respectively. 
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Fig.  27  to  30  show  the  ratio  of  residual  to  initial  strength 
against  the  number  of  thermal  cycles.  After  1170  thermal  cycles 
all  materials  exhibited  a  reduction  of  residual  strength 
regardless  of  their  curing  temperatures,  consistent  with  the 
predictable  formation  of  the  microcracks  due  to  discrepant 
thermal  strains.  In  the  high  temperature  curing  materials  the 
reduction  ceased  after  1170  thermal  cycles.  The  strength  of  the 
low  temperature  curing  materials  decreased  further  after  2295 
and  3480  cycles.  Considering  the  higher  intensity  of  prestress 
in  the  high  temperature  curing  laminates,  it  is  surprising  that 
the  reduction  of  tension  and  compression  strength  is  only  3%  and 
10%  in  laminates  164  and  175.  While  the  materials  with  lower 
curing  temperatures  deteriorated  more  both  at  23°C  and  100°C 
testing  temperatures.  Notable  is  the  drastic  strength  reduction 
of  the  polyimide-based  specimens  after  only  1170  thermal  cycles, 
requiring  additional  testing  of  this  material. 

A  comparison  of  initial  and  residual  tensile  strength  of  some 
materials  subjected  to  electron  irradiation  is  given  in  Fig.  31. 
Again  the  degradation  of  the  polyimide-based  laminate  was  the 
most  severe,  although  no  kind  of  defect  was  detectable  by  the 
applied  techniques  including  measurements  of  spectral 
reflection.  On  the  other  hand,  the  epoxy-based  laminates 
especially  the  <?14C  show  small  degrations  of  their  tensile 
strength.  Paradoxically,  the  T6T  262-12-F550  laminate  247  which 
had  been  phased  out  of  the  thermal  cycling  test  because  of  poor 
fabrication  shows  an  increase  of  tensile  strength  after  electron 
irradiation. 

8.7  Stress/Strain  Relations 

Fig.  32  to  39  show  sress-strain  plots  of  all  materials  tested 
after  different  numbers  of  thermal  cycles  and  at  two  different 
temperatures.  With  increasing  number  of  thermal  cycles  the 
stiffness  of  the  specimens  decreased  in  most  cases  especially 
during  testing  at  100°C  indicating  a  higher  toughness 
especially  for  the  high  curing  temperature  materials  due  to 
the  small  strength  reduction,  in  contrast  to  the  formation  of 
cracks  which  took  place  at  the  beginning  of  thermal  cycling  and 
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than  ceased  the  reduction  of  stiffness  progressed  continually 
with  increasing  number  of  thermal  cycles.  Changes  of  the 
molecular  structure  as  already  supposed  former  more  likely  may 
be  a  satisfactory  explanation  for  this  phenomenonthan  matrix 
cracking . 

9  .  Conclusion 

All  epoxy-based  materials  exhibited  only  small  effects  due  to 
thermal  cycling  and  electron  irradiation.  The  very  small 
strength  reduction  especially  of  the  laminates  with  high  curing 
temperatures  does  not  restrict  the  qualification  for  space 
applications  of  these  materials  at  all.  Only  the  continous 
decrease  of  the  stiffness  with  increasing  number  of  thermal 
cycles  may  be  critical  because  of  increasing  deformations  as 
far  as  the  application  for  large  structures  with  very  long 
service  times  is  considered.  But  this  shall  only  be  valid  for 
such  highly  matrix  controlled  laminates  as  the  tested 
[2(+45°  -45°)]S  laminates  were. 

The  only  polyimide  based  laminate  in  the  test  showing  severe 
damage  after  thermal  cycling  should  be  further  investigated  to 
verify  the  observations  made  until  now. 
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X-ray  Photographs  of  Three  Different  Materials  After 
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Fig.  20  Fracture  Areas  of  914C  -  Specimens  Before  and 
After  Thermal  Cycling 
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Typical  Fracture  Areas  of  HYE-2034  D  -  Specimens  Before 
and  After  Thermal  Cycling 


Fig.  23  Change  of  Natural  Frequencies  and  Damping  Characteristics 
in  Percentages  of  Original  Values  as  a  Result  of  Thermal 
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Fig.  26  Relative  Spectral  Degree  of  Reflection 
of  CFRP-Laminate  Surfaces 
a:  Prior  to  Environmental  Testing 
b:  After  Electron  Beam  Irradiation 
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FOREWORD 


This  document  is  a  compilation  of  the  more  voluminous  report 

"Notched  Strength  Predictions  of  Graphite-Epoxy  Laminates" 

(Part  I  +  II  +  6  Appendices,  DFVLR  IB  131-83/21) 

on  research  into  the  applicability  of  fracture  models  for  pre¬ 
dicting  the  notched  strength  of  graphite/epoxy  laminates. 

The  work  was  conducted  from  June,  1982,  to  September,  1983,  spon¬ 
sored  jointly  by  the  Institute  for  Structural  Mechanics  of  the 
German  Aerospace  Research  Establishment  (DFVLR)  and  in-house 
research  at  Drexel  University  in  Philadelphia.  A  Drexel  Universi¬ 
ty  support  group  comprised  of  graduate  and  undergraduate  students 
and  various  technical  assistants  participated  in  different  phas¬ 
es  of  the  program  fhe  principal  investigator  was  Dr.  Jonathan 
Awerbuch;  Dr.  Bergmann  acted  as  the  monitor  of  the  German 

Aerospace  Research  Establishment  (DFVLR). 
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Presented  in  the  Composite  Materials  Workshop,  Katholieke  Universiteit 
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INTRODUCTION 


Extensive  research  has  been  performed  in  recent  years  on  the  fracture 
behavior  and  toughness  of  composite  systems  containing  artificially  induced 
circular  holes  or  cracks.*  These  studies  focus  on  subjects  such  as 
initiation  of  crack  tip  damage  growth,  critical  crack  tip  damage  zone 
size,  notch  sensitivity,  fracture  toughness,  failure  modes  on  the  micro 
and  macro  scales,  crack  arrest  mechanisms,  etc.,  using  various  theoretical 
and  experimental  techniques.  There  are  several  reasons  for  employing  a 
variety  of  techniques  to  the  study  of  fracture  mechanics  in  composites. 

First,  different  composite  systems  demonstrate  different  failure  modes  and 
damage  mechanisms  and  may  require  correspondingly  different  analytical  tools 
and  experimental  techniques.  Second,  there  is  not  yet  a  consensus  regarding 
the  proper  set  of  criteria  for  failure.  Furthermore,  due  to  the  multiplicity 
of  failure  modes  and  the  corresponding  complexity  of  failure  processes  and 
damage  progression  in  laminated  composites,  a  variety  of  analytical  techniques 
have  been  developed  ranging  from  comprehensive  numerical  methods  to  simplified 
semi-empirical  fracture  models.  This  situation  is  undesirable  in  light  of 
the  importance  of  understanding  the  fracture  mechanisms  of  composites  before 
they  can  be  applied  to  primary  aerospace  structures.  Continuing  research 
is  of  the  utmost  importance  because  the  aerospace  applications  of  composites 
are  designed  with  mechanical  fasteners,  and  they  are  subjected  to  many  types 


*  The  terms  'crack',  'notch'  and  'slit'  are  used  interchangeably  in  the 

literature  of  composites  to  define  artificially  induced  notches.  It  should 
be  noted  that  the  typical  crack  which  forms  in  metals  under  cyclic  loading 
(and  for  which  stress  analyses  have  been  formulated)  does  not  form  in 
composites.  Consequently,  the  artificially  induced  cracks  (slits)  should 
be  termed  'notches’,  having  a  finite  notch  tip  radius.  In  this  report, 
however,  the  term  crack'  has  been  adopted  since  it  is  the  term  being  used 
by  most  researchers.  It  should  not  be  confused  with  the  actual  cracking 
which  appears  at  the  tip  of  these  artificially  induced  'cracks'  (notches) 
under  loading. 
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o £  impact  damage  (such  as  dropped  tools,  runway  9tones,  etc.)  under  various 
•loading  functions  and  environmental  conditions.  A  critical  look  at  the 
current  application  of  fracture  mechanics  In  composites  is  given  in  [1]. 

There  is  a  difference  of  opinion  whether  linear  elastic  fracture 
mechanics  (LEFM)  is  also  applicable  to  composites  containing  discontinuities. 
Initially,  research  has  been  directed  toward  the  applicability  of  LEFM  to 
composites,  however,  results  so  far  indicate  that  LEFM  can  only  be  applied 
in  limited  cases.  Several  fracture  models  have  been  proposed  in  the  past 
decade  extending  LEFM  from  metals  to  composites,  primarily  where  the 
critical  stress  intensity  factor  (i.e.  fracture  toughness)  is  being  compared 
to  that  of  traditional  structural  metals.  Also,  the  introduction  of  these 
fracture  models  was  motivated  by  the  desire  to  introduce  simplified,  easy 
to  operate  predictive  tools  for  the  notch  sensitivity  of  composites.  It 
must  be  emphasized,  however,  that  all  these  semi-empirical  fracture  models 
which  attempt  to  predict  the  notched  strength  of  composites  do  not  address, 
but  rather  by-pass,  the  micro  and  macrofailure  associated  with  the  crack 
extension  process.  The  details  of  the  actual  crack-rip  damage  zone  are 
very  complex  and  vary  among  the  different  laminate  configurations  and 

r 

material  systems.  In  the  great  majority  of  cases,  no  self-similar  crack 
growth  is  observed,  as  would  be  the  case  in  metals. 

As  mentioned  previously  LEFM  can  only  be  applied  directly  to  composites 
in  limited  cases.  Wu  [2]  determined  that  only  when  very  specific  conditions 
are  satisfied  can  the  techniques  of  isotropic  fracture  mechanics  be  directly 
applied  to  anisotropic  plates.  The  stated  conditions  are  [2J: 

1.  The  orientation  of  the  flaw  with  respect  to  the  principal  axis 


of  symmetry  must  be  fixed; 
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2.  The  stress  intensity  factors  defined  for  anisotropic  cases  must 
be  consistent  with  the  isotropic  case  in  stress  distribution  and 
in  crack  displacement  modes;  and 

3.  The  critical  orientation  coincides  with  one  of  the  principal 
directions  of  elastic  symmetry. 

The  experimental  results  obtained  by  Wu  [2]  for  unidirectional  fiber¬ 
glass  reinforced  Scotchply  with  centrally  located  cracks  in  the  direction 
of  the  fibers  (loaded  in  tension,  pure  shear,  and  combined  tension  and  shear) 
indicated  that  for  this  particular  case,  fracture  mechanics  is  adaptable  to 
orthotropic  materials.  The  use  of  unidirectional  materials  in  actual 
structural  applications  is  limited,  however;  and  when  general  multi¬ 
directional  laminates  are  concerned,  the  complexity  of  damage  growth  at  the 
crack  tip,  e.g.  Figure  1,  raises  serious  questions  as  to  the  applicability 
of  classical  fracture  mechanics  to  composites.  On  the  microscopic  level 
the  actual  failure  modes  occurring  within  the  crack  tip  damage  zone  appear 
in  the  form  of  fiber  pull-out,  matrix  micro-cracking,  fiber-matrix  inter¬ 
facial  failure,  matrix  serrations,  and/or  cleavage,  fiber  failure,  etc.  On 
the  macroscopic  level  the  major  failure  modes  are  delamination  and  matrix 
cracking  along  the  fiber  direction  in  the  individual  plies  initiating  at  the 
crack  tip,  and  failure  of  individual  plies.  An  example  of  a  typical  crack 
tip  damage  in  [O2/+45/O2/-45/O/9O]  Gr/Ep  laminate  at  two  different  load 
levels  is  shown  in  the  radiographs  of  Figure  1.  Matrix  cracking  (or  fiber- 
matrix  interfacial  failure)  along  the  0°,  90*,  and  245°  plies  is  clearly 
seen  as  well  as  the  delamination  emanating  from  the  free  edges  at  the  crack 
tips.  It  should  be  noted  that  the  failure  process,  type  of  damage  and  its 
progression  strongly  depends  on  intrinsic  material  parameters  such  as 
laminate  configurations,  material  systems,  etc.,  and  on  loading  functions 


SPECIMEN  NO.  273/1 


NOTCH  LENGTH  -  12.51 


15.50  KN  17.41  KN 


Figure  1.  X-ray  radiographs  of  center-cracked  graphite /epoxy 
(0-/+45/0-/-45/Q/90]  laminate,  25  nn  wide. 


\ 


and  environmental  conditions.  Numerous  researchers  studied  the  type  of  damage 
progression  ahead  of  the  crack  tip  and  provided  insight  regarding  the  different 
material  variables  affecting  the  crack  tip  damage  zone; 

The  presentation  is  a  brief  summary  of  a  comprehensive  review  conducted 
by  the  authors*  on  the  notch  sensitivity  of  composite  laminates  containing 
centrally  located  circular  holes  and  straight  cracks,  and  subjected  to  uniaxial 
loading.  The  complete  review  includes  a  detailed  evaluation  of  severa'l  of  the 
most  commonly  used  fracture  models  for  predicting  the  notched  strength  of 
composites,  and  of  the  experimental  results  published,  including  all  the 
relevant  available  information  regarding  material  system,  laminate  configuration 
and  geometry,  testing  procedure,  notched  strength  data,  material  mechanical 
properties,  etc.  The  review  of  both  the  fracture  models  and  experimental  data 
is  sufficiently  detailed  so  that  the  comprehensive  report  is  self-contained  and 
the  reader  does  not  have  to  refer  to  the  original  publication.  All  the  notched 
strength  data  reviewed  were  compared  with  all  the  fracture  models  reviewed  and 
the  various  material  constants  associated  with  these  fracture  models  were 
determined.  For  this  purpose,  special  computer  programs  were  developed  for 
analyzing  all  the  data  through  a  variety  of  procedures.  The  applicability  of 
the  various  fracture  models  and  their  associated  material  constants  as  well 
as  the  parameters  affecting  the  notch  sensitivity  of  composite  laminates  are 
discussed  in  detail. 

In  this  presentation  the  highlights  of  each  of  the  fracture  models  will 
be  briefly  reviewed,  samples  of  experimental  notched  strength  data  will  be 
presented,  comparison  between  predictions  and  experiments  will  be  shown,  and 
a  summary  of  the  major  material  variables  affecting  notch  sensitivity  will  be 
outlined. 

*  The  complete  reviey  can  be  obtained  from  Dr.  J.  Awerbuch,  Department  of 
Mechanical  Engineering  and  Mechanics,  Drexel  University,  Philadelphia, 
Pennsylvania  19104. 
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PREDICTIONS  OF  NOTCH  SENSITIVITY 


OBJECTIVES 


EVALUATE  ALL  EXISTING  SEMI -EMPIRICAL  FRACTURE  MODELS  TO  PREDICT 
NOTCHED  STRENGTH  OF  COMPOSITES  (ELEVEN  MODELS). 


COLLECT  AND  REVIEW  NOTCHED  STRENGTH  DATA  OF  Gr/Ep,  B/Al  AND  Gr/PI 
LAMINATES  CONTAINING  CIRCULAR  HOLES  AND  STRAIGHT  CRACKS. 


COMPARE  NOTCHED  STRENGTH  DATA  WITH  THE  FRACTURE  MODEL  PREDICTIONS 


CORRELATE  THE  VARIOUS  PARAMETERS  ASSOCIATED  WITH  THE  DIFFERENT 
FRACTURE  MODELS  WITH  NOTCH  SENSITIVITY  OF  COMPOSITE  LAMINATES 
AND  EVALUATE  THEIR  APPLICABILITY  AS  MEASURES  OF  NOTCH  SENSITIVITY 

DETERMINE  EFFECT  OF  MATERIAL  PARAMETERS  (E.G.  STACKING  SEQUENCE, 
CONSTITUENTS,  ETC.)  ON  NOTCH  SENSITIVITY. 
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\  OF  FRACTURE  MODELS 


AUTHORS 

REF, 

ABBRV. 

CRITERION 

HOLE 

SLITS 

M.E.  WADDQUPS 

J.R.  EISENMANN 
B.E.  KAMINSKI 

3 

WEK 

LEFM 

✓ 

/ 

J.M,  WHITNEY 

WN 

POINT  STRESS 

/ 

/ 

R.Y.  NUISMER 

4,5 

AVERAGE  STRESS 

/ 

✓ 

R.F.  KARLAK 

6 

K 

POINT  STRESS 

/ 

- 

R.B.  PIPES 

R.C.  WETHERHOLD 
J.W.  GILLESPIE 

7,8,9 

PWG 

POINT  STRESS 

/ 

/ 

J. W.  MAR 

K. Y.  LIN 

10,11 

ML 

N  >1-0.5 

/ 

/ 

C.C.  POE 

J.A.  SOVA 

12 

PS 

STRAIN 

/ 

- 

WEK  FRACTURE  MODEL 
(CIRCULAR  HOLES) 


Application  of  Linear  Elastic  Fracture  Mechanics  (LEFM) 
Experimental  results  indicated  thati  - - 


1.  Post-fatigue  residual  strength  ‘ 

is  equal  to  or  greater  than  the  ( 

static  strength.  if 

t 

2.  "hole  size  effect"  on  notched 
strength.  Figure  2. 

From  Irwin  [13]: 

«) 

From  Bowie  [14]  solution  for 
symmetrical  cracks  emanating 
from  a  circular  hole 

(2)  Kr  “  o”  /ira  f(a/R) 

Values  of  f(a/R)  can  be  found  in  [15] 

(3)  -  [n/a(l-v2) /E  ]  o|j  f(a/R) 

(4)  SGl  /[s4(1-v2)/E]  -  «;  f  (a/R)  -  constant 

1 

At  failure,  from  Eq.  (2): 

(5)  o"  -  yi^fc/n)]  i 

* 

For  unnotched  speclnen: 


Figure  2.  Stress  concentration 
study  [3). 


\  ■  •  OVKWM 


Figure  3.  Fracture  model  [ 3 1 . 


<s>  ••  • 


(7)  0o/o”  *  f (a/R) 

See  comparisons  in  Figure  5. 


Figure  A.  Parametric  study  (3) 
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STREN8TH  RftTIO  GTREN6TH  RATIO 


U) 


Cfc) 


9 


9 


3 


HOLE  RADIUS  <HH) 


1  _  _ 

1  3 


near*  3.  Caspar t«on  b«tv««n  axpariaanta  and  prttflctfuna  fWCK>f raetura  ao<J«l  tor:  «.  |r»Mrt/«Pxy 
(t43/0|s  laalnata;  b.  beron/alualnia  laainai*. 


WEK  FRACTURE  MODEL 
(STRAIGHT  CRACKS) 


The  concept  of  an  Intense  energy  region  at  the  tip  of  the  crack  has  been 
applied  by  WEK  for  the  case  of  composite  laminates  containing  a  straight 
crack  of  length  2c. 


From  Griffith  [16]; 

(8)  KIC  -  ojrf* 

Similar  to  the  Irvin's  plastic 
tone  correction  (13]: 


(9)  K1(,  -  oNi/rr(c  +  ac> 


ac... crack  tip  damage  zone  size  at 
failure,  i.e.,  (c  +  ac)  is  an 
"effective"  half-crack  length. 

When  c  -  0: 

(10)  K  -  o  4* 

O  C 


The  notch  sensitivity  becomes: 


(11) 


o 

o 


c  +  a 


Data  Analysis: 


(12)  c  -  ac[(oo/o’)2  -  1] 


See  comparisons  in  Figures  8  and  9 
for  Eqs .  (11)  and  (12),  respectively. 


Figure  6.  Schematic  presentation 
of  least  squares  fit 


Figure  7.  Schematic  presentation 
of  notched  strength. 
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WN -FRACTURE  MODEL. 


I.  STRESS  CRITERIA 

Point  Stress  Criterion:  Assumes  that  failure  occurs  when  the  stress 

°y>  over  some  distance,  dQ,  away  from  the 
discontinuity  is  equal  to  or  greater  than 
the  strength  of  the  unnotched  laminate  [4]: 


Figure  10,  Schematic  representation 
of  the  "point-stress" 
criterion  (5). 


Figure  1).  Schema t 1 r  representation 
of  the  "average-stress" 
criterion  [  S  )  . 
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II.  STRESS  DISTRIBUTIONS 

OO 

For  an  infinite  orthotropic  plate  subjected  to  a  uniform  stress, o  ,  applied 
parallel  to  the  y-axis  at  Infinity,  the  normal  stress,  o  along  the  x-axis 
ahead  of  the  notch  can  be  expressed  by: 


•  Circular  Hole  of  radius  R;  an  approximate  distribution  is  given  in  [21]: 
(15)  0y(x,o)  “  J-  {2+(|)2+3(|)4  -  (K^-3)[5(^)6-7(|)8])  x  >  R 


where  is  the  stress  concentration  factor  for  an  infinite  plate,  expressed 
by  [20]. 


(16)  1C,. 


1  + 


22 


^A,  ,A 


11"22  '  A12  + 


A11A22"A12 


2A 


66 


And  A^  are  the  orthotropic  ln-plane  stiffnesses  of  the  laminate  [22-26], 


1 


(17) 


(18) 


In  terms  of  effective  elastic  moduli,  Eq.  (16)  can 
form: 

*T  “  1+l2-I(gU/E22>'5-*12  +gll/g12)’i 

Using  the  standard  reference  system  adapted  for 
composites: 

K~  -  l+(  2  -  [  (E  /E  +E  /G  ]**} 

T  1  y  x7  yx  y  yx 

where  E  ,  E  ,  v  ,  and  G  are  the  effective 
x’  y  yx  yx 

elastic  moduli  of  the  laminate  for  the 
reference  system  shown  in  Figure  12. 


be  written  in  the  following 


Figure  12.  Circular  hole  in 
an  infinite  plate  subjected 
to  uniaxial  tension. 
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•  Straight  Cracks  of  length  2c;  an  approximate  distribution  is  given  in 
[18.19]  as: 

K. 


(19)  0  (x.o)  -  -  ■ 

y  /2x (x-c) 


“  /  c 
0  V  2 (x-c) 


where  the  mode  I  stress  intensity  factor  is  defined  by: 


(20) 


Kr  -  o“  ^ 


The  exact  elasticity  solution  has  been  formulated  in  [20]: 
Ktx 

(21)  o  (x.o) - l -  -  °-X- 

y  J  71  27  2  2 


xc(x  -c  ) 


X  -c 


For  large  cracks  the  approximate  expression  of  Eq.  (19)  yields  predictions 
similar  to  the  exact  solution  of  Eq.  (21).  However,  for  smaller  cracks  the 
exact  solution  should  be  applied.  The.  approximate  solution,  Eq.  (19),  is 
independent  of  crack  length  while  the  exact  solution,  Eq .  (21),  does  depend 
on  crack  length  [4,5].  The  "notch  sire  effect"  on  the  stress  distribution 
is  shown  in  Figures  13  and  14. 


Figure  13.  Normal  stress  distribution  for  Figure  14.  Normal  stress  distribution 


a  circular  hole  in  an  Infinite 
Isotropic  plate  [5). 


for  a  center  crack  in  an 
infinite  anisotropic  plate  [5] 
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III.  NOTCHED  STRENGTH  PREDICTIONS: 


A.  Circular  Holes: 

a.  Point  Stress  Criterion: 

Applying  the  "point-stress"  criterion,  Eq.  (13),  in  conjunction 
with  Ea.  (15): 

(22)  o"/ao  -  2/{2+C*+3£*  -  (K^-3)  (5£*-7£®)  } 
where 

(23)  £x  -  R/(R+d0> 

•  when  R  ■*  «*  (large  holes):  o”/oq  “  1/K^ 

•  when  R  ■*  o  (small  holes):  0fj/°0  “  1.0 

CO 

•  when  lCj.  «  3.0  (isotropic  and  quasi- isotropic  materials): 

(24)  o”/aQ  -  2/(2+  £*  +  3£*} 

•  when  R  -*•  »  (large  holes):  o^/o^  ■  1/3 

•  when  R  -*■  o  (small  holes):  c„/o  «  1.0 

N  o 

See  comparison  in  Figure  15. 

b.  Average  Stress  Criterion: 

Applying  the  "average-stress”  criterion,  Eq.  (14),  in  conjunction 
with  Eq.  (15): 

(25)  o"/oo  -  2(l-£2)/{2-£2-£2  +  (K^-3) (£*-6®) > 

where : 

£2  -  R/(R+aQ) 

when  *>  3.0  (isotropic  and  quasi-isotropic  materials): 

(26)  o“/oo  -  2(l-£2)/(2-£*-S*) 

See  comparison  in  Figure  16. 
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Figure  15.  Stress  reduction  in  quasl-isotropic  glass/epoxy  laminates  due  to  the 
presence  of  a  circular  hole,  point  stress  criterion  [4). 


Figure  16.  Stress  reduction  in  quasi-isotropic  glass/epoxy  laminates  due  to  the 
presence  of  a  circular  hole,  average  stress  criterion  (4). 


•  Since  Kj  ■  3.0  Eqs.  (24)  and  (26)  could  be  used. 

•  The  use  of  constant  values  of  dQ  -  1.02  mm  and  aQ  -  3.81  am  yield  very 


good  agreement. 


B.  Straight  Cracks 

a.  Point  Stress  Criterion 

Applying  the  "polnt-atreaa"  criterion,  Eq.  (13),  In  conjunction 
with  Eq.  (21): 

(27)  c“/oo  -  (1  - 

where 

(28)  t3  -  c/(c  +  do) 

Noting  that  for  an  infinite  plate: 

(29)  kq  - 

The  critical  atress  intensity  factor  becomes: 

(30)  Kq  -  ooV/irc(l  -  C3) 

which  is  a  function  of  crack  length.  (From  Eq.  (19):  K_  »  a  /2ird~  -  const.). 

^  o  o 

b.  Average  Stress  Criterion: 

Applying  the  "average-stress"  criterion,  Eq,  (14),  In  conjunction 
with  Eq.  (21): 

(31)  o”/oo  -  ((1  -  tA)/(l  +  E^,)]*1 

where 

(32)  E4  -  c/(c  +  do) 

The  critical  stress  intensity  factor,  Eq.  (29),  becomes: 

(33)  Kq  -  oq  /wed  -  S4)(l  +  C4) 

which  Is  a  function  of  crack  length  (From  Eq.  (19):  “  °0  ***^72  -  const.) 

From  Eqs.  (31)  and  (32): 

(34)  oJJ/oo  -  /(ao/2)/(c+ao/2) 

Similar  to  WEK-fracture  model,  Eq.  (11),  l.e.  aQ  •  2ac. 
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The  effect  of  hole  or  crack  size  on  the  stress  distribution  ahead  of  the 


hole  edge  or  crack  tip  is  shown  in  Figures  17  and  18,  respectively. 


* 

Figure  17.  Stress  distribution  at  failure,  average  Btress  criterion,  quasl- 
lsotropic  glass/epoxy  laminates  with  circular  holes  [4]. 


3 

U-C),  IN 

Figure  18.  Stress  distribution  at  failure,  average  stress  criterion,  graphite/ 
epoxy  laminates  with  center  cracks  [4]. 

•  Significant  differences  in  the  stress  distributions  between  the  two  holes  and 

crack  sizes,  using  Eqs.  (26)  and  (21,33),  respectively,  (a  -3.81  mm)  are  obtained. 

o 

•  The  results  clearly  demonstrate  the  "hole  (crack)  size  effect"  in  composite 


laminates . 


Comparisons  between  experimental  values  of  (Eq.  (29))  and  predictions 
according  to  the  "point  stress"  (Eq.  (30))  and  "average  stress"  (Eq.  (33)) 
criteria  are  shown  in  Figures  19  and  20,  respectively. 


Figure  19.  Critical  stress  intensity  factor  as  a  function  of  crack  length  for 
graphite/epoxy  laminates,  point  stress  criterion  [4).  Data  taken 
from  (28). 


Figure  20.  Critical  stress  intensity  factor  as  a  function  of  crack  length  for 
graphlte/epoxy  laminates,  average  stress  criterion  [4).  Data  taken 
from  (28J. 

•  The  use  of  constant  values  of  dQ  «  1.02  mm  and  aQ  -  3.81  mm  yield  very  good 
agreement . 


Kq  depends  on  crack  length  and  asymptotically  approaches  a  constant  value 

(dashed  lines  in  Figures  19  and  20),  equal  to  that  obtained  by  using  the 

approximate  solution,  Eq .  (19),  as  given  by:  o  /2xd  and  o  /m  /2  for  point 

o  o  o  o  * 

and  average  stress  criteria,  respectively. 
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IV. 


DATA  ANALYSES 


The  values  of  the  characteristic  distances  d  and  a  are  determined  through 

o  o 

error  analyses  to  best  fit  the  experimental  notched  strength  data. 


Linear  Regression:  Can  be  applied  in  the  following  cases. 


(35) 


(36) 


(37) 


a.  Circular  Holes,  "point-stress"  criterion;  For  laminates  having 

From  Eq .  (25) : 

R  -  d  (.6  [-1  +  (1-24(1  -  o  /(OlV11  -  l)'1 

O  1  O  IN 

b.  Straight  Cracks,  "point-stress"  criterion;  From  Eq .  (27): 

c .  do  ui  -  (o;>o)V*  -  if1 

c.  Straight  Cracks,  "average-stress"  criterion;  From  Eq.  (31): 
c  -  ao  1l[(ao/o“)2  -  1) 


1 


•  Apply  linear  regression  to  the  set  of  data  pairs  R  (or  c)  and  the  parenthetical 
terms  in  Eqs.  (35-37).  Graphically,  dQ  and  aQ  are  represented  by  the  slope  of 
the  line.  e.g.  Figure  21. 


•  Error  Minimization  Techniques:  For  laminates  containing  circular  holes  and 

having  »•  3.0  neither  criteria,  Eqs.  (2Z)  or  (25)  yields  a  closed-form 

expression  for  d  and  a  .  For  such  cases  apply  any  of  the  classical  error 
r  o  o 

minimization  techniques  to  determine  the  best  fit  of  dQ  and  aQ. 

•  Comparison  between  experiments  and  predictions  according  to  the  "point  stress" 
and  "average  stress"  criteria  for  graphlte/epoxy  and  boron/aluminura  laminates 
containing  holes,  Figure  22,  and  straight  cracks,  Figure  23,  show  excellent 
agreement.  Similarly,  experimental  results  of  Kq  (Eq.  (29))  agree  very  well 
with  predictions  according  to  both  criteria,  Figure  24.  The  characteristic 
dimensions,  d^  and  aQ,  seem  to  be  independent  of  hole  radius.  Figure  25, 

and  crack  length,  Figure  26.  20 


T  1  0 


(IIP*  •  aOR 


17«^C. 968 


CHARACTERISTIC  LENGTH  (nil)  CHARACTERISTIC  LENGTH  C  n  rl  > 


1  *68 


» 

*h 


m 


HrtLF  CRACK  LENGTH  (tin) 

fi^ur#  25.  Characcarlatlc  dlnanaiona  (W-fraccura  nodal)  as  •  function  of  half  crack  langch  for  graph: :•< 

apony  (0/90/1451,  lanlnaca  <6  -  "point  straas"  critatlon,  dQ  •  Q.51  an;  -  "avarag* 

atrasa"  critarloo.  aQ  •  2. 26  an)  indicating  that  both  can  ba  eoaaldarad  aa  notarial  conacarta 


V.  DISCUSSION 


•  BOTH  STRESS  FAILURE  CRITERIA  ARE  TWO-PARAMETER  FRACTURE  MODELS,  I.E. 

0  AND  d  OR  a  . 
o  o  o 

•  FOR  ANY  NEW  MATERIAL  SYSTEM  THEY  CAN  BE.  DETERMINED  EXPERIMENTALLY  FROM 
TESTING  UNNOTCHED  AND  NOTCHED  SPECIMENS . 

•  MINIMUM  OF  TWO  TESTS  ARE  REQUIRED,  HOWEVER,  SCATTER  IN  TEST  RESULTS 
SHOULD  BE  TAKEN  INTO  ACCOUNT  IN  DETERMINING  THE  SUFFICIENT  NUMBER  OF 
TESTS  REQUIRED. 

•  THE  CHARACTERISTIC  DISTANCES  d  AND  a  DO  DEPEND  ON  THE  SUBJECT  LAMINATE 

o  o 

•  THEY  MUST  BE  DETERMINED  EXPERIMENTALLY  A  PRIORI  FOR  EACH  MATERIAL 
SYSTEM  AND  LAMINATE  CONFIGURATION  INDEPENDENTLY. 

•  ONCE  dQ  AND  aQ  ARE  KNOWN,  BOTH  WN-FRACTURE  MODELS  CAN  BE  APPLIED  TO 
PREDICT  THE  TREND  FOR  THE  NOTCHED  STRENGTH  AND  THE  CRITICAL  STRESS 
INTENSITY  FACTORS. 

•  WHEN  THESE  CHARACTERISTIC  DIMENSIONS  ARE  PROPERLY  DETERMINED,  AN 
EXCELLENT  AGREEMENT  WITH  EXPERIMENTS  IS  ESTABLISHED  FOR  ALL  LAMINATE 
CONFIGURATIONS  AND  MATERIAL  SYSTEMS,  FIGURES  22-24. 

■  RESULTS  INDICATE  THAT  d  AND  a  ARE  INDEPENDENT  OF  NOTCH  SIZE, 

O  0 

FIGURES  25-26. 
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*OU  BU*«T»  (IN.) 


Hgure  27.  Graphical  point-stress 
characteristic  length 
determination  using  exact 
flnlte-vldth  stress  distri¬ 
bution  and  published  data  [6]. 
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Figure  28.  Plot  of  point-stress 

characteristic  lengths  from 
Figure  27  versus  hole  dia¬ 
meter  shoving  parabolic 
relationship  (61 . 


K-fracture  model  assumes  that  d^  is  related  parabolically  to  the  hole 
radius,  Figure  28,  and  through  curve  fitting  determined  that: 

U 

(40)  d  »  k  k 

o  o 

Rewriting  Eq.  (38): 

(41)  o“/oq  -  2(2  +  (1  +  d0R_1)“2  +  3(1  +  d^'W1 
Substituting  Eq.  (40): 

(42)  o”/oq  -  2(2  +  (1  +  k0R'V2  +  3d  +  k^'VV1 
COMMENTS: 

•  The  K-fracture  model  is  a  two  parameter  model,  namely  o  and  d  (or  k  ). 

o  o  o 

•  Values  of  d^  are  determined  for  individual  test  data  from  Eq.  (38). 

•  Values  of  kQ  are  determined  through  error  analysis  to  best  fit  the 
experimental  data. 

•  Linear  regression  technique  can  be  applied  for  the  case  of  IC^.  ■  3.0. 

•  This  model  could  also  be  extended  to  the  WN  "average-stress"  criterion. 
However,  it  requires  numerical  integration  of  the  stress  distribution. 

■  This  model  could  also  be  extended  to  laminates  containing  cracks,  see 
discussion  on  PWG-fracture  model. 

•  Experimental  results  obtained  by  Karlak  (6]  indicated  a  parabolic 
relationship  between  dQand  hole  radius.  Figure  29.  This  relationship 
was  applied  to  Whitney-Klm  data  ( 3 1 j  and  resulted  in  an  excellent 
agreement.  Figure  30. 

•  Comparison  between  notched  strength  data  and  prediction  shows  excellent 
agreement.  Figure  31,  also  for  the  cases  in  which  the  characteristic 

dimension  is  relatively  Independent  of  hole  radius.  Figure  32. 
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MOU  DUJrffTB  (IN.) 


Figure  29.  Tentative  plot  of  point-9tress  characteristic  lengths  versus  hole 
diameter  [6]. 


a  o.l  0J  0.}  0.4  0.1  0.4  0.7  0.1  O.f  1.0 

►*OU  DUMTta  ON. I 


Figure  30.  Comparison  of  original  Whitney-Kira  notched  strength  ratio  based  on 
constant  point-stress  characteristic  length  (dashed  curves)  and 
modification  proposed  here  (solid  curves)  [6). 


PWG-FRACTURE  MODEL 
(CIRCULAR  HOLE) 


I .  INTRODUCTION 

•  Modified  WN-"point-stress"  criterion  for  laminates  containing  circular 
holes  (ICj  i*  3.0). 

•  it  is  assumed  that  the  characteristic  distance,  d  ,  depends  on  hole 

o 

radius . 

•  An  exponential  relationship  between  dQ  and  hole  radius  is  assumed. 

•  Notched  strength-radius  superposition  principal  Is  proposed  which 
allows  superposition  of  all  laminates  and  materials  onto  a  single 
master  curve. 

•  For  this  purpose  two  radius  shift  parameters  are  defined  which  depend 
on  the  notch  sensitivity  of  the  subject  laminate. 

•  Relative  notch  sensitivity  is  discussed. 


II.  NOTCHED  STRENGTH  PREDICTIONS: 
Rewriting  Eq .  (15): 


(43)  o  /o“  -  (1/2)  (2  +  (R/x)2  +  3(R/x)4  -  (K^  -  3)[5(R/x)6  -  7(R/x)8]} 


where 


(44)  kT  -  1  +  { 2[  (E  /E  J*5  ~  V  +  E  /G 

T  1 '  y  x7  yx  y  yx 


Assump t ion : 

(45)  d  a  Rra 
o 
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(46) 


(47) 


(48) 


UI. 


Introducing  a  notch  sensitivity  factor,  C: 

d  «  (R/R  )B/C  R  ...reference  notch  radius  Introduced  so 

o  o  o 

that  (R/R  )  is  nondloensional. 
o 

Using  the  "point-stress"  criterion,  Eq.  (13),  together  with  Eqs.  (43)  and  (46): 

o>  -  2  (2  +  (X)2  +  3(X)4  -  (l£-3)[5(X)6  -  7(X)8])-1 
No  1 

where 

X  “  [1  +  r""1  R0“m  C-1)-1 

and  R  is  a  reference  notch  radius,  for  algebraic  simplicity  chosen  equal 
o 

to  1.0  ln-  Thus,  PWG-fracture  model  is  a  three-parameter  model  (o  ,  C,  m) . 


EFFECT  OF  PARAMETERS  ON  NOTCH  SENSITIVITY: 


The  higher  the  value  of  C,  the  more  notch  sensitive  the  material,  Figure  33: 

•  C  -*■  0:  notch  insensitivity 

•  C  -  °n/0o  ‘  1/K^ 


for  other  values  of  m  different  notch  sensitivity  curves  will  be  obtained. 


•  The  exponential  parameter  is  bounded  between  0  <  m  <  1,  Figure  34: 

•  m  »  0:  WN  "point-stress"  criterion  is  recovered 

•  m  «  Is  notched  strength  is  independent  of  hole  radius 

•  R  >  1.0:  notch  sensitivity  decreases  with  increasing  radius  (for 

C  •  10  in  2) .  Therefore:  PWG-fracture  model  is  applicable 
for  R  f  1.0  inch.  However:  increasing  C  will  shift  the 
cross-over  point  to  the  right. 
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Figure  33. 

Influence  of  notch 
factor  upon  notched 
(7). 

jensltivity  Figure  34. 

strength 

Influence  of  exponential 
parameter  upon  notched 
strength  [7}. 
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Figure  35. 

Influence  of  stress  concentration 
factor  on  notched  strength  (7). 

l*||g  0  IlStlHl 

Figure  36.  Influence  of  processing  on 
notched  strength  (7j. 

Figure  37.  Influence  of  stacking  sequence 
on  notched  strength  {7]  (data 
taken  from  [31 ) ) . 

•  The  stress  concentration  factor,  affects  notch  sensitivity.  Figure  35: 

•  ranges  in  most  case3  between  2.0  <  <  6.0 

•  For  R  <  1.0  in,  k!J  -  3.0  (isotropic)  is  a  good  approximation  (for 

C  «  10  in  ^  and  a  “  0.5). 

•  COMMENT:  Note  that  the  parameters  C,  m,  and  ICj,  are  coupled.  Consequently, 

the  actual  values  of  these  parameters  for  a  given  set  of  notched 
strength  data  may  not  necessarily  Indicate  the  notch  sensitivity 
of  the  subject  material.  A  relative  notch  sensitivity  parameter 
has  been  Introduced,  see  Section  V.  Excellent  agreement  between 
prediction  and  experiments  has  been  established  by  PWG  (7J,  Figure 
36  and  37.  Similar  data.  Figures  38  and  39,  also  show  very 
good  agreement. 
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STRENGTH  <HPA>  STRENGTH  RATIO 


ri«ut#  M.  Coaparlion  between  uptrlMntl  and  prediction  (PWC-f ractute  aodal  for  er*phite/epo*Y  [0/145/901 
laailnaea  (■  -  0.37,  C  -  0.45,  -  1.43),  1 


Ftfure  39.  Comparison  between  expertoenta  and  prediction  (pWC~fracture  modal)  for  §r aphlt*/epo*y  lO/j-5/aOl 
lfliute  '■  -  0.37,  C  -  0.43,  «M  •  1.437 . 


Using  the  two  radius  shift  parameters  ac  and  a^  a  single  oJJ  -  R  curve  can 
be  constructed  for  all  laminate  configurations  and  material  system  (providing 
they  all  have  the  same  k£) ,  e.g.  Figures  42-43  and  Table  1. 


Figure  42.  Master  curve  for  graphite/  Figure  43.  Master  curve  for  boron/ 

epoxy  [8],  Data  taken  from  aluminum  [8], 

Ref.  [3,31]. 


Data  in  Figures  42  and  43  are  taken  from  Figures  36  and  37,  respectively. 


Table  1.  Material  Superposition  Parameters  [8] 


B/Al  I  [0J6  0.36 
B/Al  II  [0] 6  0.50 
B/Al  III  [0)6  0.50 
T300/934  (*45/0/90]a  0.48 
T300/934  (90/0/±45)s  0.40 


Morganite  11/4617  [0/445]^  0.15 


29.4 

1.28 

8.64 

3.607 

10.0 

1.00 

1.00 

3.607 

8.6 

1.00 

0.74 

3.607 

7.5 

1.04 

0.56 

3.000 

11.5 

1.20 

1.32 

3.000 

40.0 

1.70 

16.00 

2.943 
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c.  m  -  constant,  C  •  constant  (K^  is  variable): 

To  superimpose  all  laminates  of  different  orthotropy  into  a  single 

Act 

master  curve  which  correspond  to  -  3.0 
From  Eq .  (A7)  for  K^.  -  3.0 


(59) 

00  A  2  AA  —  1 

<7n/<7q  -  212  +  A  +  3X  } 

A  A  Act 

X  (i.e.  R  )  for  Kj.  -  3.0 

*  4  *.2 

CO 

(60) 

( 3 X  )  +  (X  )  +  2(1  -  1/S)  -  0 

V°o  *  S 

(61) 

X*2  -  (-1  +  (1  -  24  (S-l)  /S  J*5)  /6 

* 

Solving  Eqs.  (59)  and  (61)  for  R  : 

(62)  R*  -  c1/(m'1)  {[(-1  +  (1  +  24(1-S) /S)*4)/  fe]-*5  -  j}1/ <"»-!) 


Consequently,  the  strength  data  of  any  laminate  configuration  and  material 

*  CO 

system  can  be  superimposed  into  a  common  master  curve  of  lij  <•  3.0  and  having 

*  * 

arbitrarily  preselected  values  of  C  and  m  ,  Figure  44. 


Figure  44.  Master  curve  for  all  materials  (6)  (from  data 
shown  in  Figures  42  and  43). 
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V .  RELATIVE  NOTCH  SENSITIVITY : 

To  quantify  the  degree  of  relative  notch  sensitivity  of  the  subject 
material.  Figure  45. 

As  a  reference  condition,  a  notch  insensitive  material  is  chosen,  i.e. 

m*  -  0.0,  C*  -  1.0  inch-1: 

(ft3>  RnS  _  1O810R*  -  lo«10R  (dB 

Using  R  -  acmR  (Eq-  (51)) 

(64)  log1()R*  -  log1()R  -  log10actn 

:r,5)  Rns  "  log10acm 

Using  Eq.  (50)  with  R  •  0.1  inch  (for  algebraic  simplicity),  m*  »  0.0, 

C*  -  1.0  inch-1,  R  -  1.0  inch: 
o 

(66)  Rns  -  m  +  log1QC 

COMMENT :  R^s  is  defined  for  laminates  having  the  s aae  When  laminates 

of  unequal  orthotropy  are  involved,  the  degree  of  notch  sensitivity 
should  be  determined  from  Eq.  (47). 

VI.  DATA  ANALYSES: 

A  computer  program  routine  scans  through  all  possible 
combinations  of  C  and  m  and  chooses  the  values  which 
result  in  the  smallest  error  between  the  expected  data 
and  the  predictions  according  to  Eq  (47).  A  straight¬ 
forward  linear  regression  program  can  e  employed  for  the 
special  case  of  laminates  containing  circular  holes  and 
having  ■  3.0. 
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Values  of  the  relative  notch  sensitivity,  R^  ,  for  the  six  material  systems 

discussed  previously,  Figures  42  and  43,  are  shown  in  Table  2  {8J.  It  should 

be  recalled  that  R^g  is  defined  for  material  systems  having  the  same  K^,  thus 

a  comparison  of  the  elastic  notch  sensitivity  values  can  be  made  only  among 

the  three  types  of  boron/aluminum  or  among  the  three  graphite/epoxy 

laminates.  Comparison  between  the  R  values  of  the  two  material  systems 

ns 

00 

is  meaningless  since  they  have  different  Kj.  values. 

Table  2.  Relative  Notch  Sensitivity  [8] 


Material 

m 

C 

fin.'1] 

R 

ns 

B/Al  111  (0*]6 

0.50 

8.6 

1.43 

3.607 

B/Al  II  [0’)6 

0.50 

10.0 

1.50 

3.607 

B/Al  I  [0*16 

0.36 

29.4 

1.83 

3.607 

T300/934  [i45/0/90]9 

0.48 

7.5 

1.36 

3.000 

T300/934  [90/0/145] s 

0.40 

11.5 

1.46 

3.000 

Morganlte  11/4617  [0/t45]2a 

0.15 

40.0 

1.75 

2.943 

Notched  strength  data  were  compared  with  prediction  and  an  excellent 
agreement  has  been  established.  Figure  46.  The  master  curve  to  which  all 
notched  strength  data  were  shifted.  Figure  46,  was  selected  from  Ref.  [7], 
i.e.  m*«0.5  and  C*“1.0  inch"*”0.4  mm-*.  Figure  47  shows  the  master  curve 

s  to  h 

defined  by  m*-0.5,  C*”0.4  mm“  ,  !ij»3.0  and  the  shifted  notched  strength  data 
of  the  different  graphite/epoxy  laminates  containing  circular  holes  which 
were  analyzed  in  this  review. 
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STREHGTH  RATIO 


LAt-UP:  l02/t*51# 


-  9  '  .  +  "Orl*lnal"  notched  atran|ch 

|  □  "Shifted*  ootcb«d  atrangth 


Pradlctlon  (tq.  (*0) 


HOLE  RADIUS  (nn) 

Fl|ur«  46.  Motchad  atrangch  varaua  hola  radiua  pradictlona  and  uacir  curva  (r%JOf  raccura  nodal) 
for  traphita^apoay  (Oj/14 3 J %  lanloaca. 


HOLE  RADIUS  CM 


PWG-FRACTURE  MODEL 
(STRAIGHT  CRACKS) 


The  subsequent  formulation  is  analogous  to  that  developed  by  PWG  for  laminates 
containing  circular  holes. 

I.  NOTCHED  STRENGTH  PREDICTIONS : 

From  the  exact  elasticity  solution,  Eq .  (21) 

(67)  "  x(x2  -  c2)  **  x  >  c 
Assumption: 

(68)  d  «  c“ 

o 

Introducing  a  crack  notch  sensitivity  factor 

(69)  d^  -  (c/co)m/K  Cq. .. reference  half-crack  length  introduced 

so  that  (c/c  )  is  nondimensional 
o 

Using  the  "point-stress"  criterion,  Eq.  (27),  together  with  Eq.  (69) 

(70)  o”/oo  -  (1  -  X2)*1 

where 

(71)  Xj  “  [1  +  c”"1  c~m  K-1]'1  (analogous  to  Eq.  (48)) 

and  cq  is  a  reference  half-crack  length,  for  algebraic  simplicity  chosen 

equal  to  1.0  inch. 

Thus,  PWG-fracture  model  is  a  three-parameter  model  (o  ,  C,  m) . 

o 
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Since : 


(72)  d  -  cm/K 


(73)  l°810do  ‘  ■  lo«ioc  “  lo810K 


(see  Figure  43) 


Using  0 


x-c+d 


«  0  together  with  Eq.  (67): 
o 


,2,»a 


N  o 


Solving  for  dQ: 


(75)  do  -  c{  [1  -  (0"/o0)2]'!i-l 


The  values  of  K  and  m  can  be 
obtained  by  a  least  squares  fit 
of  the  data.  These  can  be  used 
to  predict  the  notched  strength. 
From  Eqs.  (72)  and  (74): 


(76)  o”/ao  -  {1  -  (1  +  c"  1  (TV2)*1 


The  values  of 
K  -  19.0  inch-1 
and  m  “  0.4  fit 
the  data  best. 
Figure  49. 


*  log  C 

Figure  48.  Relations  between  character¬ 
istic  dimensions  and  the  notr'r 
size  for  T 300/ 9 34  [90/0/14 5 )s 
[9].  (Data  taken  from  [31))- 


Figure  49. 
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Inappropriateness  of  a  constant  charac¬ 
teristic  dimension  for  the  graphite/ 
epoxy  T300/934  [90/0/145]  material 
system  [9J.  (Data  taken  from  Ref.  [31).) 
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II.  EFFECT  OF  PARAMETERS  ON  NOTCH  SENSITIVITY: 


•  The  higher  the  value  of  K,  the  more  notch 
sensitive  Che  material.  Figure  50. 


Figure  50.  Influence  of  the  constant  K  on  the 
notched  strength:  m  ■  0.2  (9). 


•  The  exponential  parameter  is 


decreases  with  in-  Figure  51.  Influence  of  the  constant  m  on  the 

notched  strength:  K  “  10  19]. # 


increasing  crack 
length  (for  K»10) . 

The  higher  the  value  of  K  the  larger  the  crack  length  for  which  the  cross-over 
occurs.  Figure  51.  44 


Notch  sensitivity  curves  in  linear  scales  are  shown  in  Figures  52  and  53. 


Comparisons  between  experiments  and  predictions  Indicate  very  good 
agreement  for  the  case  in  which  the  characteristic  dimension  dQ  is  Independent 
of  crack  length.  Figures  54  and  55.  Excellent  agreement  is  also  established 
however  when  dQ  does  depend  on  crack  length.  Figures  56  and  57. 
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r Ijur •  34.  Comparison  bicvttn  axpar  Irenes  and  prediction  (FW>f  ractura  modal)  for  (raptilea /epoxy  (0/90/543) 
laiuci  (•  -  -0.1,  K  •  2.30.  k  -  46.*). 


Figure  33.  Characteristic  dimension*  (FWG-f recture  model)  aa  a  function  of  half  crack  length  for  graphite/ 
ap«*y  [ 0 / 90 / *43 )  laminate  («  -  -0.1,  l  •  2.30,  k  -  *6.6). 
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Flfur*  56.  Co«p«rlaoa  bacveen  txpirlMacf  ind  prediction  (PVC~t rocturo  nodal)  tot  boron/ el  uninu*  (0/t4JJ 
lalnttt  (•  “  0.23,  t  •  0.62,  -  26.9). 


» 


HALF  C  F  A  C  K  LENGTH  (MH) 


Ftpire  57.  Characteristic  dlMasion  (FVG-f  rectors  nodal)  •«  •  function  of  holt  crock  lon|th  tor  boron/ 
oluelouo  (0/t43|t  loolnote  (a  •  0.23,  t  •  0.62,  •  26.9). 


COMMENT:  Note  that  the  parameters  K  and  m  are  coupled.  Consequently, 


the  actual  values  of  these  parameters  for  a  given  set  of  notched 
strength  data  may  not  necessarily  indicate  the  notch  sensitivity 
of  the  subject  material.  A  relative  notch  sensitivity  parameter 
has  been  introduced. 


HI.  NOTCHED  STRENGTH  -  CRACK  LENGTH  SUPERPOSITION  METHOD: 


■  To  superimpose  all  notched  strength  data  onto  a  single  master  curve 

*  * 

which  is  defined  by  preselected  values  of  K  and  m  . 

•  For  this  purpose  a  crack  length  shift  parameter,  a  ,  is  defined 

Km 

as  follows: 

A  necessary  and  sufficient  condition  for  superposition  of  a  point  on  the  notched 

*  i 

strength  curve  defined  by  K  and  m  to  a  new  point  on  a  curve  defined  by  K  and  u> 

★ 

is  that:  X^  -  X^,  i.e. 

(77)  (c*)(°  _1)  (co)-<n  (K*)"1  -  (c)  (corm  (K)'1 


solving  for  c  : 

(78)  c*  -  (K*/K)a  (cQ)8  (c)  ” 


Define: 


(79)  a^  -  (K*/K)a  (co)S  (c)  " 

(80)  c*  -  a^c 


*  ^ 

where:  a  -  l/(m  -1) 

*  * 

8  “  (m  -m) / (m  -I) 

a  -  (m-1) / (m*-l) 
in 


c* . . . "new"half-crack  length  on  the 

master  curve  defined  by  K*  and 
* 

m  which  are  selected  arbitrarily. 
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EFFECT  OF  PARAMETERS  ON  CRACK  LENGTH  SHIFT  PARAMETER: 


(Use  c  -  1.0  inch) 
o 


(81) 


(82) 


a.  K  •  cons tan Fi£u£e_58£ 


.  (m-1) / (m  -1) 


b.  m  ■  £ons£ant^  £i^u£e_58^: 
(83)  c*  -  aKc' 


(84)  aR  -  (K*/K)1/(m  _1) 


Combining  Eqs.  (82)  and  (83): 


*  m 

(83)  c  -  aKc 


Analogous  to  the  procedure  shown  for  laminates  containing  circular  holes 


Figure  58.  Superposition  principle  for  composite  material  system  (32). 


Figure  59.  Data  for  graphite-epoxy  T300/934 
system:  Q  ,  [±45/0/90]s,  m-0.4, 
K-13. 5; 0,1 90/0/145 JB,  m-0.4, 
K-19.0  [9]  (data  taken  from 
Ref.  [31]). 


Figure  60.  Data  for  glass-epoxy  Scotchp) y 
1002  material  system:  A, 
[0/±45/90]2b,  m-0.2,  K-22.0 
[9]  (data  taken  from  Ref.  [53]). 


Figure  61.  Slit  notch  master  curve  (-): 

Q  ,  T300/934,  [l45/0/90]8, 
m-0.4,  K-13. 5;  O,  T300/934 
[90/0/145]. ,  m-0.4,  K-19.0; 

4,  Scotchp ly  1002,  [0/l45/90]2s, 
m-0.2,  K-22.0 


PWG  [9]  analyzed  notched  strength  data  of  graphlte/epoxy  and  glass 
epoxy  as  shovn  In  Figures  59  and  60,  respectively.  Using  Eq.  (79),  the 
data  for  all  three  laminates  were  superimposed  into  a  single  master  curve 
with  pre-selected  m*-0.5  and  K*“  10  inch'*,  Figure  61,  and  an  excellent 
agreement  has  been  established. 

Similarly,  the  "original"  notched  strength  data  of  unidirectional 
boron/aluainum  can  be  shifted  to  a  pre-selected  master  curve  (m*-0.5,  Re- 
10  inch'*).  Figure  62.  The  restuls  of  shifting  the  notched  strength  data  of 
all  boron/aluminum  laminates  containing  straight  cracks  are  shown  in 
Figure  63.  50 


IV.  RELATIVE  NOTCH  SENSITIVITY 


•  To  quantify  the  degree  of  notch  sensitivity  of  the  subject  material. 

Figure  58. 

•  As  a  reference  condition,  a  notch  insensitive  material  is  chosen,  i.e. 
m*  -  0.0,  K*  -  1.0  inch-1: 

•  The  relative  notch  sensitivity  is  obtained  by  calculating  the  crack  length 
shift  for  a  reference  half-crack  length  of  0.1  inch. 

(86)  *  aKm  (m*  ”  °-0,  K*  "  1,0  inch_1>  c  "  0.1  inch,  CQ  “  1-0  Inch) 
Substituting  Eq.  (86)  into  Eq .  (79)  yields: 

(87)  Rns  "  log10  ®Km  “  “  +  log10K  (analogous  to  Eq.  (66)) 

Values  of  a^  for  the  three  material  systems  analyzed  in  Figures  59-60 
are  listed  in  Table  3. 

Table  3:  Generalized  Notch  Sensitivity  Factor  [ 9 J . 


Material 

Reference 

m 

K 

aKm 

Graphite/Epoxy 

T300/934 

(t45/0/90]s 

31 

0.4 

13.5 

33.9 

Graphite/Epoxy 

T300/934 

[90/0/iA5]s 

31 

0.4 

19.0 

47.7 

Glass/Epoxy 
Scotchply  1002 
t0/±45/90)2s 

5 

0.2 

22.0 

34.9 

V,  DATA  ANALYSES 

A  straightforward  linear  regression  program  can  be 
employed  from  which  the  least  squares  fit  of  K  and  m 
is  determined  by  using  Eq.  (73).  The  scatter  in 
results  is  then  Illustrated  as  shown  in  Figures  55 
and  57. 
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MODIFIED  PWG-FRACTURE  MODEL 
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40 


SMC-R50 

HMC 

HMC/GR 


Figure  64.  Error  as  a  function  of  notch  Insensitive  crack  half  length  [32J. 


Ci  (mm) 


Figure  65.  Crack  notch  sensitivity  factor  as  a  funccion  of  notch  insensitive 
crack  half  length  [32). 


K  (mm"1) 


Figure  66.  Exponential  parameter  aa  a  function  of  notch  Insensitive  crack 
half  length  [32]. 


Replacing  K,  m,  and  cc  in  Eq  . 

(89)  by  K1, 

V  *nd  cc- 

respectively: 

(92)  oN/oo  -  <1  -  (1  +  cc 

C1)-2^ 

* 

c  >  ct 

analogous 

°N/Oo  '  K0 

* 

c  *  ci 

to 

Eq.  (76) 

•  Is  now  always  positive 

» 

•  dy  Increases  with  increasing  half-crack  length 

•  comparison  between  experiments  and  predictions  is  shown  in  Figures  67  and  68. 


•  The  modification  enables  the  determination  of  the  size  of  the  notch 
Insensitive  region. 


55 


PWG  and  WN  fracture  models  do  not  account  for  a  notch  Insensitive  region 
(o”/oQ  -  1.0  at  c  ■  0.0) 

Modified  PWG-fracture  model  does  account  for  a  notch  insensitive  region 

(<yoo  ”  1,0  at  c  “  c*)»  the  oJj/Oq  axis  is  shifted  to  c*.  The 

*  * 

values  of  Kj  and  are  calculated  after  eliminating  the  notch  Insensitive 
region. 

The  modified  PWG-fracture  model  is  also  applicable  for  the  cases  of 
m  >  0,  Figure  69. 


Figure  69.  Notched  strength  of  HMC/graphite  as  a  function  of  crack  half  length 
for  two  values  of  K  and  m  [32]. 
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•  The  relative  notch  sensitivity  factor  becomes  (see  Eq.  (87)) 

<93>  <si  ■  “i  +  lo*10K*i 

where  the  superposition  method  is  applied  in  a  similar  manner  to  that 

previously  described.  However,  the  crack  length  shift  would  be  applied 

*  *  *  * 
from  the  K^,  curve  to  the  master  curve  of  K  ,  m  ,  rather  than  from 

the  initial  K,  m  curve. 


II.  CIRCULAR  HOLES: 

The  formulation  of  the  modified  PWG-fracture  model  for  laminates  containing 
circular  holes  is  identical  to  that  proposed  for  laminates  containing 
straight  cracks  (Eqs.  (88)  to  (93)). 


(94)  R*  -  R  -  R* 
c  i 


*  *  mi  * 

(95)  dQ  -  (Rc/Ro)  Vct 


R  >  R. 


R  >  R, 


R. ..  .critical  notch  insensitive- 
radius 


Rc- . .critical  corrected  hole 
radius 


* 

d  -  0 
o 


R  <:  R. 


Replacing  C,  m  and  R  in  Eq.  (47)  by  C^,  n^,  and  Rc»  respectively: 

(96)  o"/oo  -  2{2  +  X*  +  -  (t£-3)[5X®  -  7X®]}'1  R  >  R* 


V°o  “  1‘° 


R  <  K 


where 

*  mi  "  1  (R  i"*1  rcV1)"1 

(97)  \±  -  [1  _  (Rc)  1  (V  (Ci)  ] 
and: 

(98)  RnSi  "  mi  +  1°8io('i  (similar  to  that  given  by  Eqs.  (66)  and  (93)) 
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ML-FRACTURE  MODEL 

(CIRCULAR  HOLES  AND  STRAIGHT  CRACKS) 


•  The  LEFM  equation  applied  to  homogeneous  materials  is: 

(99)  a“  -  KIC^ 

in  which  the  exponential  is  the  order  of  the  mathematical  stress 
singularity  at  the  tip  of  the  crack. 


ML~fracture  model  proposed  that  the  fracture  of  composites  Is  governed  by: 


(100) 


°N  ’  HC(2C) 


^...Composite  Fracture  Toughness  [MPa  (mm)  ],  a  property  of  the  subject 
material 


n. . ."The  order  of  singularity  of  a  crack  with  its  tip  at  the  interface  of 
two  different  materials" 

•  The  fracture  model  is  based  on  stress  analysis  for  the  case  of  cracks  at 
the  Interface  of  a  bimaterlal,  i.e.  matrix  and  fiber.  Figure  70. 

■  The  order  of  singularity  is  a  function  of  the  ratio  of  the  shear  moduli  of 
the  matrix  and  filament,  and  the  two  Poisson's  ratios,  and  V2> 

Table  4  (10,11]. 

From  Eq.  (100) : 

(101)  log1(J  (o”/oo)  -  log10(Hc/co)  -  n  log10(2c) 

(Ig  and  n  are  determined  from  a  least  squares  fit  using  linear  regression  and 
results  are  plotted  on  a  log  (oJJ/o^ versus  log(2c/W)  or  log  (2R/W)  format, 
Figure  71, 
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MATRIX 


1 


Figure  70.  Positioning  of  discontinuity 

at  oatrix/fiber  interface  (33]. 


Figure  71.  Schematic  plot  on  log-log  format 
of  ML  -  fracture  model. 


Table  4.  Order  of  Singularity  [11] 


wl/u2 

V1 

V2 

n 

0.01 

0.30 

0.20 

0.250 

0.01 

0.35 

0.20 

0.269 

0.25 

0.30 

0.15 

0.263 

0.25 

0.30 

0.20 

0.262 

0.25 

0.30 

0.25 

0.261 

0.25 

0.35 

0.20 

0.280 

0.50 

0.30 

0.20 

0.280 

0.50 

0.35 

0.20 

0.297 

0.10 

0.30 

0.20 

0.310 

0.158 

0.30 

0.20 

0.339 

0.20 

0.30 

0.20 

0.357 

0.10 

0.33 

0.20 

0.319 

0.158 

0.33 

0.20 

0.347 

0.20 

0.33 

0.20 

0.364 

Comparisons  between  prediction  and  experiments  in  logarithmic  and  linear 
scales  are  shown  in  Figures  72-73  for  graphite/epox;'  laminate  containing 
a  circular  hole  and  in  Figures  74-75  for  boron/aluminum  laminate  containing 
a  straight  crack. 
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STREN8TH  RATIO 


I  -J-l- 


_l_ 


9.1  9.Z 


i.9  z.9  i.»  i» .  e».  ■>». 

HOLE  RADIUS  <rtn> 


Figure  72. 


Coaparlaon  bacvaan  exparlaaaca  and  prediction  (ML-fractura  soda  1)  for  graphlta/apoxy  (0/*<o/90)i 
laaiaaca,  a  hewn  on  a  logarithmic  acala  (a  -  0.173,  “  0.773). 


Figura  73.  Coaparlaon  bacvaan  experiment a  aod  prediction  (ML-fractura  aodal)  for  graphite /epoxy  (O/tAJ/901 
laminate,  above  oa  a  Liaaar  acala  (a  •  0.173,  "  0.V73). 


3* 


HALF  CRACK  LENGTH  tnn) 

FLfura  71.  -Cm* arisen  kacvaon  aaforlaanta  aa4  prediction  (ML-f ractura  aodal)  for  boron/alualnua  (0,/i*ll 
laalnaco.  thorn  m>  a  LlMar  ocala  (a  -  0.2*0,  •  0.171).  2  * 
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PS-FRACTURE  MODEL 
(STRAIGHT  CRACKS) 


I.  INTRODUCTION 


DEVELOPED  FOR  LAMINATES  WHICH  EXHIBIT  SIGNIFICANT  NONLINEAR  STRESS-STRAIN 
CURVES,  e.g.  [t45]s  Gr/Ep,  B/Al,  ETC. 


^  •  DATA  ARE  ANALYZED  IN  TERMS  OF  STRAIN  RATHER  THAN  STRESS. 

•  A  "GENERAL  FRACTURE -TOUGHNESS  PARAMETER",  Qc,  IS  PROPOSED. 

•  Qc  IS  DERIVED  ON  THE  BASIS  OF  FIBER  FAILURE  IN  THE  PRINCIPAL  LOAD¬ 
CARRYING  LAMINAE,  I.E.  IS  INDEPENDENT  OF  LAMINATE  ORIENTATION. 

•  THE  RELATIONSHIP  BETWEEN  Kq  AND  Qq  DEPENDS  ONLY  ON  THE  ELASTIC 
CONSTANTS  OF  THE  LAMINATE. 

•  BASED  UPON  Qc,  THE  CRITICAL  STRESS  AND  STRAIN  INTENSITY  FACTORS  AND 

_  NOTCHED  STRENGTH  OF  INDIVIDUAL  LAMINATES  CAN  BE  CALCULATED. 

•  Qc  IS  PROPORTIONAL  TO  THE  ULTIMATE  TENSILE  STRAIN  OF  THE  FIBERS. 

•  CONSEQUENTLY,  THE  CRITICAL  STRESS  INTENSITY  FACTOR  CAN  BE  PREDICTED 
SOLELY  FROM  THE  ULTIMATE  TENSILE  STRENGTH  OF  THE  UNNOTCHED  LAMINATE. 

•  THE  CRACK  TIP  DAMAGE  ZONE  SIZE  CAN  BE  PREDICTED  BASED  UPON  Qc. 

•  PS-FRACTURE  MODEL  IS  BASED  UPON  A  STRAIN  FAILURE  CRITERION. 
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II.  STRAIN  FAILURE  CRITERION 

Stress  and  Strain  Intensity  Factors: 

(102)  K  -  S  /re  sec (wc/W)  (103)  Ke  -  eq  /rc  sec(rc/W) 

Applying  Crack  Tip  Damage  Zone  Size: 

(104)  K  -  S  /rr(c+p)  sec  (wc/W)  (105)  Ke«  cq  /r(c+p£)  sec(irc/W) 

Critical  Stress  and  Strain  Intensity  Factors: 


(105)  Kq  -  Sc  /r(c+pc)  sec( rc/W) 


(T07) 


oC 


/r(c+p^J  sec  (rc, 


) 


Cr ideal  Crack  Tip  Damage  Zone  Size: 

VFtu)2 

(108)  pc  -  ■  -  -  (109)  PeC 

<U0>.  ^  -V[1-KQe/(l,cF2tu>^',1  <1U>  K«j 


(1VEtu> 


W1 


^eQe'^tu” 


Elastic  Critical  Stress  and  Strain  Intensity  Factors: 


(112)  KQe  -  Sc  /rc  sec(rc/W)  (113)  KgQe  -  cqC  *Vc  sec  (rc/W) 


For  Linear  Elastic  Material 


(114) 


*  Although  the  subject  of  finite  width  correction  (FWC)  factor  is  addressed  in  detnil 
in  the  following  Section,  it  has  been  incorporated  into  the  PS-fracture  model  so 
that  the  formulation  will  be  identical  to  that  developed  by  Poe  and  Sova  [12). 


EoC  "  VEy 


eQ 


Vef 


tu 


F„  /  E 
tu  y 


eC 
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DEFINITIONS: 


K,  K£ . stress  and  strain  Intensity  factors,  respectively. 

. . .critical  stress  and  strain  intensity  factors,  respectively. 
K^e,K^e. .elastic  critical  stress  and  strain  intensity  factors,  respectively, 
p,  Pe- . crack  tip  damage  zone  size  for  stress  and  strain,  respectively. 


PC’  peC 


....critical  crack  tip  damage  zone  size  for  stress  and  strain, 


respectively. 


S,  . . remote  (far-field)  stress  and  strain,  respectively. 

o 

S  .  £oC... remote  (far-field)  stress  and  strain  at  failure,  respectively. 

F  ,  £_  ..unnotched  ultimate  tensile  stress  and  strain,  respectively, 
tu  tu  1 


In  order  to  apply  a  strain  failure  criterion,  the  strains  at  the  notch  tip  are 
derived  based  on  stress  analysis  (IS]  for  an  anisotropic  and  homogeneous  laminate 
containing  straight  cracks.  The  singular  stresses  ahead  of  the  crack  tip  for  a 


specially  orthotropic  laminate  subjected  to  model  I  loading  are,  (for  6  »  0, 


Figure  76) : 


Due  to  symmetry,  -  0  and 
the  strains  are  expressed  by: 


\  ♦  ’  I 


I  I  .  ♦ 

Figure  76.  Laminate  and  principal  laminae 
coordinates  [12]. 


(117)  (6) 


f  1/E 


-v  /E 
yx  y 


-v  /  E 
yx  y 

1/E 


1/G 


0 

0 

xy 


Transforming  the  singular  strains  to  the  principal  directions  of  the  1th 
lamina 


E1 

E 

X 

2 

sin  a 

2 

cos  a 

sln2a 

(118) 

E2 

-  (TJt 

e 

y 

where: 

(119)  [T]  - 

cos^a 

2 

sin  a 

-sin2a 

*5112 

1 

Y 

xy 

-*5Sin2ci 

*28^2(3 

-cos2r 

Substituting  Eq.  (116)  into  Eq.  (118) 


'  1 

E1 

K 

SB  — 

E  *r2wr 

V 

V 

W’1 

(120) 

E2 

C2 

where:  (121) 

C2 

-  (T]1(6) 

Ey 

> 

y 

'12 

1 

C12 

1 

0 

ASSUMPTION:  "A  laminate  falls  whenever  the  fiber  strains  reach  a  critical 
level  in  the  principal  load  carrying  laminae". 

Thus:  ij  *$xr  -  CONSTANT  (AT  FAILURE) 

i.e.  Independent  of  Laminate  Lay-Up 

CONCLUSION:  KqUjJj/E  -  CONSTANT  (at  Failure,  t*  Rq> 

(122)  DEFINITION:  Qc  -  "GENERAL  FRACTURE-TOUCHNESS  PARAMETER" 

Special  Case:  Linear  elastic  material,  KEq  ■  Kp/Ey,  thus: 

(123)  Qc  -  K,q(C1)1 
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The  value  of  (C^)^  has  to  be  determined  individually  for  each  principal 
load-carrying  lamina  using  Eq.  (121),  e.g. 


(124) 

a  -  O’: 

(Cl),  -  1  -  v  /E  /Z 
"l'i  yx  x  y 

(125) 

a  -  45*: 

(C,),  "  >5(1  -  v  </z  /E  )(1  +  /e  /E  ) 
"11  v  yx  x'  y' '  y'  x' 

(126) 

o  «  a: 

(C,),  -  (1  -  v  /E  /E  )(/E  /E  •  sin' 
"1  i  v  yx  x  y  y  x 

III. 

DATA  ANALYSIS 

A.  Th£  ,^Generalj2F£a£ture  Tougtin^sjj  £airameter"j_ 

Qc  can  be  calculated  in  three  different  vays,  i.e.: 

1.  Qc  -  keq^i^i  KtQ  def®rmlned  from  Eq.  (Ill),  efcu  and  are  given 

2.  Qc  »  Kq  (Cp)  Kq  determined  from  Eq.  (110),  Ftu  and  Sc  are  given 

3.  Qc  •  Kp(c^)i/Euy  E . . .ultimate  secant  modulus. 

•  Values  of  Qc  calculated  from  the  first  procedure  are  Intermediate 
between  the  upper  and  lower  values  calculated  from  the  second  and 
third  procedures,  respectively.  Figure  77  and  Table  5. 

•  The  latter  two  procedures  yield  equal  values  for  linear  materials. 


With  increasing  nonlinearity  In  the  stress-strain  curves  the  difference 
between  the  upper  and  lower  bounds  would  be  larger. 


w 


T 


Av«r«9« , 


1 


Figure  77.  Values  of  the  "General  Fracture-Toughness  Parameter"  for  Boron/ 
Aluminum  Laminates  [12]. 


Table  5.  "General  Fracture-Toughness  Parameter",  Qc>  for  Five 
Boron/Aluminum  Laminates  [12]. 


Qc,ifC  calculated  from  - 


Laminate 

Orlnetation 

<1 

51KeQ 

ClKQ/Ey 

ClKq/Euy 

<°16T 

0.8409 

0.01328 

0.01200 

0.01347 

[02/i«5]g 

0.7841 

0.01068 

0.00764 

0.01223 

[145/02], 

0.7806 

0.01074 

0.00771 

0.01233 

[0/i45]g 

0.7375 

0.01250 

0.00770 

0.01477 

f±*5>2a 

0.6771 

0.02156 

0.00350 

0.1001 
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B.  Notched  Strength_Predictions2 


The  notched  strength  can  be  calculated  In  three  ways  all  based 
on  the  assumption  that  Qc  for  a  given  material  system  is  constant. 

1.  Kq  -  and  apply  Eq.  (110)  to  calculate  Sp. 

2.  K-  -  Q-E  /£.  and  apply  Eq.  (110)  to  calculate  S„. 

y  L  uy  i 

3.  K  „  -  Q  /£,  and  apply  Eq.  (Ill)  to  calculated  e  and  use  the 

CQ  U  1  Ovy 

stress-strain  equations, 


The  predictions  of  notched  strength  using  Keq  values  agree  well  with 
experiments,  while  those  using  with  E  and  Ey  yield  the  upper  and 
lower  bounds,  respectively,  Figure  78. 


•  The  difference  between  the  upper  and  lower  bounds  predictions  is  largest 

for  the  [i45],  laminate,  and  generally  Increases  with  the  percentage  of 
2s 

the  i45°  plies.  Figure  78. 

C.  Cr_i tical_Damage  Zone_Sizej 

•  The  critical  crack  tip  damage  tone  sizes,  p^.  and  PeC>  can  calculated 
from  Eqs.  (108)  and  (109),  respectively,  when  the  values  of  and  Keq 
are  back-calculated  from  Eqs.  (122)  and  (123),  respectively  using  the 
average  values  of  Q^. 


•  Notched  strength  data  analyzed  in  [40]  for  a  variety  of  center-crack 
(0i/±45j/90k)8  laminates  (Gr/Ep,  B/Ep,  Gl/Ep,  Cr/Pl,  B/M. ,  etc.) 
indicate  that  Qc<  ^  is  reasonably  constant,  i.e.: 

(127)  Qc/etuf  ;  1.5/nm 
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provided  that: 

•  Laminate  failed  predominantly  by  self-similar  crack  extension 

•  Neither  delamiqation  nor  splitting  in  the  0°  plies  at  the  crack 
tip  occurs. 

•  Crack  tip  damage  zone  size  is  small. 


Using  Eqs.  (109)  and  (123): 


n 


(128) 


eC 


1.5c 


tuf 


Ctu 


i.e.,  critical  damage  zone  size  is  not  constant  but  depends  on  laminate 
orientation  (O^.  and  therefore  on  the  laminate  elastic  constants. 


o 


•  The  critical  fiber  strain  ahead  of  the  crack  tip  in  the  load-carrying 

lamina  is:  "  Qc/ /2irr . 

•  Assuming  chat  at  r  »  dQ  (WN-point  stress  criterion). 

Thus:  /2^T  -  Qc/etuf_ 

i.e.  the  characteristic  distance  in  the  "point  stress"  criterion  of  the 
WN-fracture  model  is  proportional  to  the  square  of  Q^/etuf 

•  Since  dQ  varies  significantly  among  the  different  Gr/Ep  laminates 

analyzed  (see  WN-fracture  model)  having  the  same  (“0.01)  it  should 

be  expected  that  varies  among  the  different  laminates  as  well. 

•  However,  if  the  conditions  listed  above  are  satisfied  it  could  be  expected 
that  Qc  is  independent  of  laminate  configuration. 

•  For  a  given  laminate,  a  constant  value  of  Q^, ,  independent  of  crack  length, 
results  in  an  excellent  agreement  with  experiments,  Figures  79-80. 


i 
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STRENGTH  RftTIO 


3 


HALF  CRACK  LENGTH  (HH) 

Figure  79.  Notched  acrength  veraua  crack  length  predictions  (FS-fractura  nodal)  for  graph Ice /epoxy 

(0/*43)t  lamina ca  (aolld  llna:  ^  of  Cha  subject  laminate;  daahad  llna:  obtalnad  from 

avaraga  of  values  given  In  Tabla  16). 


figure  80. 


mm  cr.ck  l.nith  prMietlm,  (M-frmtur«  aod.l)  for  Ir*phlt./.|>o«T 
Umrnt.  (mild  Urn:  t  .1  tk.  mkl.ee  lisln.c;  d.ikxl  lln.i  0,  obc.lo.d 
from  avaraga  of  valuea  given  la  Table  16).  ' 


D.  Discussion: 


(129) 


(130) 


(131) 


(132) 


In  order  to  fully  apply  the  PS-fracture  model  the  elastic  constants 
(EX-  V  V’  str2n8ths  (Ftu.  SC>'  strains  (e^,  etuf,  e^)  and 
ultimate  secant  modulus  (EUy)  must  be  known. 


•  Most  experimental  data  published  do  not  include  several  of  these 

material  properties,  predominantly  the  values  of  Euy  and  eo(;.  . 

•  Consequently,  only  calculation  procedure  02  (for  calculating  Qc> 
and  calculation  procedure  B 1  (for  calculating  Kq  and  Sq)  can  be 
employed. 


•  For  most  Gr/Ep  laminates,  this  may  be  sufficient  since  their  stress- 
strain  curve  is  essentially  linear  to  failure. 

•  Thus,  from  Eqs.  (110),  (112),  and  (122)  [40): 

Sc  /sec  (irc/M)/Ftu  -  (1  +  ircU^/Q^)2}"''5 
using  Qc  -  1.5  ceuf 

Sc  /sec  (rc/W)/Ftu  -  [1  +  wIc^/l.Sc,.^]2)^ 

For  linear  stress  strain  curves;  F  •  e  .E  : 

tu  tuf  y 

Sc  /sec  (ttc/W)/FCu  “  [1  +  rc(51/1.5)2]  ** 

which  yields  the  prediction  of  the  notched-to-unnotched  strength 
ratio  for  linear  materials. 

For  nonlinear  laminates,  applying  Eqs.  (Ill),  (113),  and  (123): 
toC  /sec  (irc/W)/etu  -  {1  +  l'2)*1 

and  using  the  Ramberg-Osgood  stress-strain  equations  or  any  polynomials 

yields  the  strength  ratio. 
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II.  ORTHOTROPIC  FINITE  WIDTH  CORRECTION  FACTOR 


v:% 


•  There  are  several  analytical  methods  to  determine  FWC,  e.g.: 

•  Boundary  collocation  solution  [45] 

.  Boundary  integral  equation  [46] 

•  Finite  elements  [47] 
etc . 

•  No  closed-form  solutions  are  available,  and  the  FWC  for  orthotropic 
materials  require  numerical  solutions. 

•  Therefore,  the  closed-form  expressions  available  for  isotropic 
materials  are  frequently  used  for  orthotropic  materials. 


•  Validity  of  the  applicability  of  isotropic  FWC  factor  to  orthotropic  materials 
is  discussed  below: 

H  *  Y(orthotropic)  ^(isotropic)  “  KI(orthotropic) /KI (isotropic) 

•  The  larger  the  length-to-width  ratio  (L/W)  is,  the  closer  H  approached  1.0  [47]. 

•  When  L/W  >  3.0,  the  error  for  [O’]  Gr/Ep  is  less  than  2*  [48]. 

•  Local  and  global  compliances  of  center  notched  [O']  B/Al  and  BSiC/Tl 
Indicate  that  Isotropic  FWC  factor  (Eq.  (139))  can  be  used  [49]. 

•  Values  of  H  for  a  variety  of  Gr/Ep  T300/5208  laminates  were  calculated 

in  [50]  for  DEC  and  CCT  specimens  using  the  boundary  integral  formulation. 

(2a/W  -  0.1  to  0.7,  L/W  -  5.0) 


•  H  depends  on  specimen  geometry,  lay-up,  and  material  properties. 
Figures  81  and  82 
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SECTION  REDUCTION  I2./W 


As  Che  "free  edge  distance"  becomes  larger.  U  approaches  unity, 
regardless  of  geometry  and  material. 

H  is  largest  for  [i45)g,  decreases  (symmetrically)  for  increasing  and 
decreasing  ply  angle. 

A  correlation  between  H  and  the  ln-plane  shear  modulus  has  been  suggested. 
Figures  83  and  84  ,  [50]. 

For  all  practical  purposes  2c/W  <  0.5  for  which  H  <  1.12. 

H  is  further  reduced  for  multidirectional  laminates,  e.g.  Table  6. 

H  was  calculated  for  four  Gr/Ep  laminates:  [902/i45]g,  [90/±45]g,  [ 0/ ±A5 ) s » 

and  [i45js,  for  laminates  containing  circular  holes,  center-cracks,  and 
cracks  emanating  from  circular  and  elliptical  cutouts,  all  of  various 
sires  [46]. 

<  For  [±45)2g:  H  s  1.06  for  2c/W  -  0.5 
•  For  all  other  laminates  H  s  1.02  for  2c/W  ■  0.5 


Consequently,  in  this  work  the  Isotropic  FWC  factors  are  used,  i.e.  Eqs. 

(135)  and  (139)  for  laminates  containing  circular  holes  and  straight 

cracks,  respectively.  For  PS-fracture  model,  Eq.  (137)  has  been  employed. 

Table  6.  Anisotropy  factor  in  CCT  specimens  of  mixed 
orthotropic  laminates (a) 


(b) 

[0/±45]  1.002  1.006  1.012  1.019  1.028  1.037  1.045 

[04/±45!g(b)  1.001  1.001  1.000  0.999  0.998  0.996  0.994 

a.  Data  taken  from  Table  5  of  [50] 

b.  HTS  graphlte/apoxy .  Lamina  properties  are:  E^  “  144.8  GPa,  “  11.7  CPa, 
Glt  -  4.5  GPa,  -  0.28 


c.  Center-cracked  tension  specimens 
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Figure  83.  Anisotropy  factor  versus  relative  shear  modulus  for  double-edge- 
cracked  specimens  150]. 
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Figure  84.  Anisotropy  factor  versus  relative  shear  modulus  for  center- 
cracked  specimens  [50]. 
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COMPARISON  between  experiments  and  predictions 


OBJECTIVES 


COLLECT  AND  REVIEW  NOTCHED  STRENGTH  DATA  OF  Gr/Ep,  B/Al,  AND  Gr/PI 
LAMINATES  CONTAINING  DISCONTINUITIES;  i.e.  Circular  Holes,  Center 
Cracks,  and  Impact  Damage*. 


•  COMPARE  NOTCHED  STRENGTH  DATA  WITH  THE  FRACTURE  MODEL  PREDICTIONS. 

•  CORRELATE  THE  VARIOUS  PARAMETERS  ASSOCIATED  WITH  THE  DIFFERENT  FRACTURE 
MODELS  WITH  NOTCH  SENSITIVITY  OF  COMPOSITE  LAMINATES  AND  EVALUATE  THEIR 
APPLICABILITY  AS  MEASURES  OF  NOTCH  SENSITIVITY. 

•  ANALYZE  THE  EFFECT  OF  A  VARIETY  OF  MATERIAL  PARAMETERS  ON  NOTCH 
SENSITIVITY  OF  COMPOSITE  LAMINATES. 


•  THE  COMPARISONS  AND  CORRELATIONS  AND  THE  SUBSEQUENT  CONCLUSIONS  ARE 
BASED  ON  A  REVIEW  OF**. 

•  2800  NOTCHED  STRENGTH  DATA 

•  20  MATERIAL  SYSTEMS 

•  62  LAMINATE- CONFIGURATIONS 

•  ALL  UNIAXIAL  TENSILE  LOADING  IN  AMBIENT  CONDITIONS 

The  9ystem  of  data  collection  and  a  sample  of  the  review  procedure  are  shown  in 
the  following  7  pages. 

Material  nomenclature,  lay-ups,  ei:c.  are  listed  in  the  following  pages. 


DATA  FILING  SYSTEM 


File  Names 

Each  set  of  notched  strength  data  collected  from  the  literature  was 
assigned  a  number  for  the  purpose  of  storing  the  data  in  the  computer  in  the 
fora  of  data  files.  This  number,  or  File  Name,  is  comprised  of  nine  digits 
which  represent  five  types  of  information: 

1.  The  first  two  digits  represent  the  material  system,  as  listed 
in  Table  7. 

2.  The  third  digit  represents  the  test  temperature,  as  listed  in 
Table  8. 

3.  The  next  two  digits  (fourth  and  fifth)  represent  the  laminate 
configuration,  or  lay-up,  as  listed  in  Table  9. 

A.  The  sixth  digit  represents  the  shape  of  the  discontinuity,  as 
listed  in  Table  10. 

5.  The  last  three  digits  represent  the  reference  number  from  which  the 
set  of  data  was  taken,  as  numbered  in  the  list  of  references. 

For  example.  File  Name  110044H.056  represents: 

1.  Digits  11  -  from  Table  7;  graphite/epoxy  T300/5208. 

2.  Digit  0  -  from  Table  8;  21®C  (70“F)  -  room  temperature. 

3.  Digit  44  -  from  Table  9;  (0/90/±45]g  laminate. 

4.  Digit  H  -  from  Table  10;  containing  a  circular  hole. 

5.  Digits  056  -  from  the  reference  list;  Reference  No.  56. 

Summary  nomenclature  of  the  file  name  is  given  below: 

X  X  X  X  X  x.xxx 

-  Reference 

- : - Discontinuity  Shape 

-  Lay-Up 

-  Test  Temperature 

-  Material  System 

The  file  n3me  appears  at  the  top  of  each  Figure. 
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Table  7.  File  Name;  Material  System  Table  8.  File  Name;  Test  Temperature 


Material  System 

File  Name(a) 

Temperature 

File  Name^ 

I. 

Graphite/Epoxy 

‘C 

(*F) 

Not  Available ^ 

10 

21 

(70)  (b) 

0 

T300/5208 

11 

10 

(50) 

A 

AS1/3501-6 

12 

38 

(100) 

B 

T300/934 

13 

66 

(150) 

C 

T300/934C 

14 

93 

(200) 

D 

SP-286T300 

15 

121 

(250) 

E 

AS/3501-5 

16 

149 

(300) 

F 

AS/3501-6 

17 

177 

(350) 

G 

204 

(400) 

H 

II. 

Boron/ Aluminum 

232 

(450) 

I 

/V\ 

260 

(500) 

J 

Not  Available' 

2X 

288 

(550) 

K 

5.6/6061F 

21 

315 

(600) 

L 

8. 0/1100F 

22 

343 

(650) 

M 

5.6/ 2024F 

23 

371 

(700) 

N 

4.0/6061F 

24 

399 

(750) 

0 

4.2/6061F('C) 

25 

427 

(800) 

P 

4.2/2024F(c) 

26 

-46 

(-50) 

Q 

5.7/6061F<c) 

27 

-73 

(-100) 

R 

5 .6/606116 

28 

-101 

(-150) 

S 

5.6/2024T6 

29 

-129 

(-200) 

T 

-157 

(-250) 

U 

III. 

Graphite /Poly imide 

-184 

(-300) 

V 

Cb^ 

-212 

(-350) 

w 

Not  Available 

3X 

• 

Celion  6000/FMR-15 

31 

a. 

Flr3t  two  digits  in 

File  Name. 

a.  Third  digit  in  File 

Name. 

b. 

For  cases  in  which 

constituents 

b .  Room 

temperature. 

are  not  outlined  in 

Publication. 

c.  BS1C/A1 
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Table  9 .  File  Name  ;  Laminate  Lay-Up 


Lay  Up 

File  Name(a) 

Lay  Up 

File  Name 

t°«n 

01 

[902/0]a 

26 

(i!5]s 

02 

[+15/0/-15]s 

27 

[±30]s 

03 

[+30/0/-30] g 

28 

i^j]g 

04 

[+45/0/-45] g 

29 

[i6°]s 

05 

[+60/0/-60] g 

30 

[i?3]s 

06 

[+75/0/-75]g 

31 

\  ±90  j 

n 

07 

[ 90 / 0 / 90 ] g 

32 

;o/t5]3 

08 

[0/±45]2s 

33 

£0/±10]s 

09 

( 90/±45  j  g 

34 

[ 0/  ±15  ]  s 

10 

[+45/0/-45/0]  g 

35 

( 0/±30] g 

11 

[0/90]s 

36 

[0/±45]g 

12 

[0/90 ] 2g 

37 

[  0  /  ±  b  0  ]  g 

13 

[0/90]4s 

38 

14 

[0/±45/90]g 

39 

[0/ :902] s 

15 

[0/±45/90]2s 

40 

[i5/0Ja 

16 

[±45/0/90Js 

41 

[ ±10/0] s 

17 

[ 90/0/±45] s 

42 

[±15/0] g 

13 

[90/±45/0]g 

43 

] ±20/0 ] 3 

19 

[0/90/±45 ) 

44 

[ ±25 /0 ] ^ 

20 

[ ±45/90/0  ]g 

45 

r  ±  30/0 1 

S 

21 

[ ±45/0/90] 2g 

46 

r±35/°]g_ 

22 

[ 90/0/±45 ] 2g 

* 

47 

r±45/o)g 

23 

[+45/0/-45/90]2g 

48 

[ ±60/0] 

24 

( 90 / ±45/0] 2g 

49 

[±75/0] 

s 

25 

[  0/ 90/ ±4  5 ] 2g 

50 

(a) 
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Table  9.  Continued 


Lay  Up  File  Name(a) 

[+A5/90/-45/0] 2g  51 

[02/±45]s  52 

[±45/02)s  53  . 

[0/+45/90/-45]..  34 

ZS 

[90/-45/0/+45]2a  55 

'Q/+45/90/-451a  56 

[+45/0/-45/0]s  57 

[0/+45/-45/0]  58 

s 

[0/90/±452]s  59 

[0/i22)g  60 

[ 90/ ±22 ]  61 

[0/-22/0] s  62 


(a)  Fourth  and  fifth  digits  in  File  Name. 


Table  10.  File  Name:  Notch  Shape 


Notch  Shape 


File  Name^ 


Circular  Center  Hole  H 
Straight  Center  Crack  C 
Impact  Damage  (Central)  X 
Circular  Center  Hole  with  Slits  S 
Double  Edge  Notch ^  D 
Single  Edge  Notch^  E 


(a)  Sixth  digit  in  File  Name 

(b)  Not  included  in  this  report. 
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SAMPLE  LITERATURE  REVIEW 

Reference  [64] 

Author:  I.M.  Daniel 

Title:  The  Behavior  of  Uniaxially  Loaded  Graphite/Epoxy  Plates  with 
Holes . 

Publication:  in  Proceedings  of  the  First  International  Conference  on 

Composite  Materials  (ICCM-I),  Boston,  U.S.A.,  19  78,  pp.  1019-30V:. 

File  Names:  15039H.064,  15052H.064 


Materials:  Graphite/Epoxy  SP-286T300 
Lay-Up:  [0/±45/90]g,  [02/±45] s 

Dimensions:  560  x  127  x  1.05  mm  (gage  length  ■  408  mm) 

Geometry:  Circular  Hole  -  0.0,  0.2,  0.4,  1.6,  4.8,  6.4,  12.7,  and  25.4  mm 
in  diameter 

Tabs:  Glass/Epoxy,  3M  1007  Scotchply,  [0/90)g  laminate,  76  mm  long. 

No.  of  Tests:  2-3  tests  per  crack  per  laminate;  total  of  28  notched  tests. 

Test  Conditions:  Stroke  control,  1  mm/min,  friction  grips,  Riehle  testing 
machine . 


© 


Measurements:  Local  (near  crack  region)  strain,  coating,  and  notched  strength 


Lamina  Properties:  Not  available  (taken  from  Ref.  [58], 


Stress  Concentration  Factor,  K^.: 

[0/±45/90] a :  -  3.00 

[02/i45]g:  -  3.49 


Comments:  1.  Unnotched  strength  for  both  laminates  is  not  available.  For 
[0/145/90]  laminate  the  value  given  in  Ref.  [58]  was  taken 
for  the  data  analysis  of  this  report  since  it  is  of  the  same 
material  system  as  that  studied  in  this  paper  [64].  For 
(0^/145  ]g  laminate  the  value  given  in  Ref.  [66]  of  802  MPa  is 
taken  for  the  data  analysis  of  this  report. 

2.  References  [65,66]  provide  additional  details  on  test  results. 
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Results  and  Conclusions: 

The  author  investigated  the  deformation  and  failure  of  uniaxially 
loaded  graphite/epoxy  plates  with  holes  and  correlated  the  notched 
strength  data  with  the  "point  stress"  and  "average  stress"  riteria  of 
the  WN-fracture  model.  The  main  conclusions  are: 

1.  The  "point  stress"  and  "average  stress”  criteria  satisfactorily 
described  the  strength  reduction. 

2.  There  is  a  critical  hole  diameter  below  which  the  laminates 
become  notch-insensitive. 

3-  In  the  case  of  quasi-isotropic  laminates,  strength  reduction  was 
found  to  be  independent  of  notch  geometry,  i.e.  specimens  with 
holes  and  cracks  of  the  same  size  had  nearly  the  same  strength. 

4.  For  {0/i45/90]  laminate:  a  *■  3.8  mm  d  ■  1  ■ 
so  o 

For  [0,/£45J„  laminate:  a  -  5  mm 
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RESULTS  AND  CONCLUSIONS 


•  VERY  GOOD  AGREEMENT  BETWEEN  ALL  FRACTURE  MODELS  REVIEWED  AND  ALL 
EXPERIMENTAL  NOTCHED  STRENGTH  DATA  HAS  BEEN  ESTABLISHED. 

•  THE  FRACTURE  MODELS  ARE  ALL  SEMI-EMPIRICAL,  i.e.  THEY  CAN  BE  APPLIED 
PROVIDING  THAT  AT  LEAST  2-3  (Depending  Upon  the  Fracture  Model)  NOTCHED 
STRENGTH  DATA  ARE  KNOWN  A  PRIORI. 

•  THE  FRACTURE  MODELS  INCLUDE  CERTAIN  PARAMETERS  WHICH  ARE  ASSUMED  TO  BE 
MATERIAL  CONSTANTS,  i.e.  Independent  of  Specimen  Geometry. 

■  THESE  PARAMETERS  STRONGLY  DEPEND  ON  LAMINATE  CONFIGURATION  AND  MATERIAL 
SYSTEM  AS  WELL  AS  ON  THE  VARIETY  OF  INTRINSIC  AND  EXTRINSIC  VARIABLES. 

•  CONSEQUENTLY,  THESE  PARAMETERS  MUST  BE  DETERMINED  (EXPERIMENTALLY)  FOR 
EACH  LAMINATE  CONFIGURATION  AND  MATERIAL  SYSTEM  INDEPENDENTLY. 

^  •  THE  NUMBER  OF  TESTS  THAT  MUST  BE  CONDUCTED  IN  DETERMINING  THE  FRACTURE 

MODELS  PARAMETERS  DEPENDS  UPON  THE  MODEL  ITSELF  AND  THE  LEVEL  OF  ACCURACY 
REQUIRED. 

•  LN  THE  FRACTURE  MODELS  REVIEWED  THE  ACTUAL  PATTERN  AND  DETAILS  OF  THE  NOTCH 
TIP  DAMAGE  ARE  BY-PASSED  BY  SIMULATING  THE  DAMAGE  AS  SOME  "EFFECTIVE"  NOTCH 
TIP  DAMAGE  ZONE  AND  ASSUMING  IT  TO  GROW  IN  A  SELF-SIMILAR  MANNER. 
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THE  VARIOUS  PARAMETERS  ASSOCIATED  WITH  THE  DIFFERENT  FRACTURE  MODELS 
DEPEND  ON  A  VARIETY  OF  INTRINSIC  AND  EXTRINSIC  VARIABLES  SUCH  AS:  • 


INTRINSIC  VARIABLES 

*  LAMINATE  CONFIGURATION 

*  STACKING  SEQUENCE 

*  CONSTITUENTS  PROPERTIES 

*  Fiber  volume  fraction 

*  FIBER-MATRIX  INTERFACE 

*  FABRICATION  PROCEDURE 


EXTRINSIC  VARIABLES 

•  LOADING  FUNCTION 

•  LOADING  RATE 

•  SPECIMEN  GEOMETRY 

•  SHAPE  OF  DISCONTINUITY 

•  TEST  TEMPERATURE 

•  MOISTURE  CONTENT 


HOWEVER,  A  COMPREHENSIVE  EVALUATION  OF  THE  EFFECT (S)  OF  THESE  VARIABLES  ON 
THE  NOTCH  SENSITIVITY  OF  COMPOSITE  LAMINATES  IS  STILL  LACKING,  e.g.  Effect(s) 
of  matrix  toughness,  fiber  strength,  etc. 


ALSO,  ANY  COMPARISON  OF  NOTCH  SENSITIVITY  AMONG  DIFFERENT  LAMINATES  OBTAINED 
FROM  DIFFERENT  SOURCES  IS  OF  QUESTIONABLE  VALUE,  i.e. 

•  Very  few  publications  report  fiber  volume  fraction,  environmental  test 
conditions,  fabrication  procedures,  constituents  (fibers)  properties,  etc 

VERY  FEW  WORKS  ATTEMPT  TO  CORRELATE  THE  RECORDED  NOTCH  SENSITIVITY  WITH  THE 
OBSERVED  FAILURE  PROCESSES  AND  FAILURE  MODES  PRIOR  TO  CATASTROPHIC  FRACTURE. 


SINCE  THE  APPLICABILITY  OF  THE  FRACTURE  MODELS  REVIEWED  DEPENDS  ON  THE 
FAILURE  PROCESSES,  e.g.  delamination,  splitting,  size  of  damage  zone,  etc., 
ADDITIONAL  ATTENTION  TO  THIS  ISSUE  IS  WARRANTED. 
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THERE  IS  NOW  A  GROWING  DEMAND  FOR  LIGHT-WEIGHT  COMPOSITES  IN  "PRIMARY" 
STRUCTURES.  CONSEQUENTLY,  A  MORE  THOROUGH  UNDERSTANDING  IS  REQUIRED  OF 
THE  WAY  IN  WHICH  PERTINENT  LOADING,  ENVIRONMENTAL  AND  MATERIAL  PARAMETERS 
AFFECT  STRENGTH,  TOUGHNESS,  FRACTURE  BEHAVIOR,  FATIGUE  CHARACTERISTICS, 
IMPACT  RESPONSE,  ETC.,  OF  COMPOSITE  LAMINATES. 


?RACTICALLY  ALL  THE  LAMINATES  REVIEWED  ARE  HIGHLY  NOTCH  SENSITIVE.  FOR 
MANY  CASES,  THE  NOTCHED  STRENGTH  REDUCES  BY  AS  MUCH  AS  50*  FOR  NOTCH-LENGTH- 
TO-WIDTH  . ATIOS  OF  0.2  ♦  0.3. 

MOST  OF  THE  EXPERIMENTAL  STUDIES  ON  NOTCHED  STRENGTH  OF  COMPOSITE  LAMINATES 
ARE  LIMITED  TO  TWO  TYPES  OF  NOTCH  GEOMETRIES:  CIRCULAR  HOLES  AND  STRAIGHT 
CRACKS,  CENTRALLY  LOCATED. 

THE  MAJORITY  OF  THE  NOTCHED  STRENGTH  DATA  ON  GRAPHITE/EPOXY  ARE  FOR 
LAMINATES  CONTAINING  CIRCULAR  HOLES,  WHILE  FOR  BORON /ALUMINUM  MOST  OF 
THE  DATA  ARE  FOR  LAMINATES  CONTAINING  STRAIGHT  CRACKS. 
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EVALUATION  OF  RESULTS*  . 

EFFECT  OF  LAMINATE  CONFIGURATION:  The  experimental  notched  strength  data 
indicate  that  laminate  lay-up  strongly  affects  the  notched  strength  and 
notch  sensitivity  of  graphlte/epoxy  and  boron/aluminum  laminates,  e.g. 

Figures  85-90  and  Figures  91-92,  respectively. 

EFFECT  OF  LAMINATE  STACKING  SEQUENCE:  Most  of  the  notched  strength  data  indicate 
that  laminate  stacking  sequence  does  affect  the  notch  sensitivity  of  graphite/ 
epoxy  laminates,  Figure  93,  most  probably  resulting  from  differences  in  the 
failure  processes.  No  sufficient  data  are  available  to  draw  any  conclusions 
for  boron/aluminum  laminates,  Figure  94. 

EFFECT  OF  NOTCH  TIP  RADIUS:  It  has  been  determined  that  notch  tip  radius  (or 
shape  of  discontinuity) has  little  effect  on  notch  sensitivity  of  graphite/epoxy, 
Figure  95,  and  boron/aluminum.  Figure  96  .  This  characteristic  has  frequently 
been  cited  in  the  literature  [208,216,254-256],  and  it  is  attributed  to  the 
notch  tip  blunting  during  damage  progression. 


© 


EFFECT  OF  IMPACT  DAMAGE:  Both  resin  matrix  and  metal  matrix  composites  are 
highly  sensitive  to  impact  damage,  even  when  it  is  nonvlsual  damage.  The 
lateral  damage  is  the  major  cause  of  strength  degradation.  Results  indicate 
that  induced  lateral  damage  can  be  presented  as  through-the-thickness  straight 
cracks  for  post-impact  residual  strength  predictions.  Figures  97-98. 


EFFECT  OF  TEST  TEMPERATURE:  Very  little  data  are  available  regarding  the  effect 
of  temperature  on  the  notch  sensitivity  of  composite  laminates.  Results  shown 
in  Figures  99  -100  indicate  that  at  elevated  temperature  boron/aluminum  can 
become  more  notch  sensitive,  which  might  be  due  to  strength  degradation  of  the 
fibers . 

*  In  the  following  discussion  only  the  predicted  notch  sensitivity  curves  (WN- 
fracture  model,  "average  stress"  criterion)  are  plotted.  The  actual  experimental 
data  are  excluded  for  the  sake  of  clarity. 
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CRACK  LENGTH  (2C/U) 

Figure  91.  Notched  strength  versus  crack  langth  predictions  (WN-fracture  nodal,  "average  acrata" 
criterion)  for  savaral  5.6  all  boron/aliaiinua  6061F  laminates. 


CRACK  LENGTH  C  2C/U > 

Figure  92.  Notched  atrengch  versus  crack  langth  predictions  (UN-fraeture  modal,  "average-stress”  crltet ion) 
for  aevaral  5.6  ail  boron/aluaini«i  606 IF  laminates. 
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Flgura  93.  Nocchad  atrtngch  ««riui  notch  length  pradlctiona  (WN-fractura  mod  ,  "avaraga  atrata 
crlcarlon)  for  graphlca/apoxy  (0/S45/90)  laminates  containing  circular  holaa  and 
straight  cracks. 
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EFFECT  OF  CONSTITUENTS :  Significant  differences  in  notch  sensitivity  (and 
notched  strength)  were  obtained  for  a  seemingly  identical  material  system, 
Figure  1C1 .  Thdse  differences  are  attributable  either  to  different  fabrication 
procedures  [209]  or  to  differences  in  fiber  strength.  Figure  102.  Improved 
notch  sensitivity  of  unidirectional  boron/aluminum  can  be  achieved  through  a 
proper  choice  of  constituents.  Effect  of  constituents  on  notch  sensitivity. 
Figure  103,  indicates  that  the  combination  of  large  fiber  diameter  and  ductile 
matrix  leads  to  improved  notch  sensitivity  characteristics.  When  two  different 
aluminum  alloys,  e.g.  6061F  and  2024F,  are  chosen  as  the  matrix  material,  very 
little  effect  on  notch  sensitivity  has  been  recorded,  Figure  104.  Effect  of 
heat  treatment  and  type  of  matrix,  on  notch  sensitivity  and  deformation 
characteristics  of  boron/aluminum  are  reported  in  [229,234,238].  Similar 
systematic  studies  on  graphite/epoxy  are  not  conclusive. 


•  EFFECT  OF  LAMINATE  CONFIGURATION  OH  FRACTURE  MODEL  PARAMETERS:  The  effect  of 

lamination  angle  on  the  characteristic  dimension,  aQ,  (WN-average  stress 

criterion)  for  angle-ply  [ i9 ] laminates  is  shown  in  Figure  105,  indicating 

that  the  largest  value  (i.e.  the  least  notch  sensitivity)  occurs  for  [ ±45 1 2S 

laminate.  With  the  addition  of  a  O’  ply  (e.g.  [0/i9]g,  [±9/0]s,  and  [+9/0/-6]g) 

the  dependence  of  a  on  9  Is  reversed,  i.e.  the  value  of  a  is  smallest  for 
o  o 

0  =  45",  Figures  106-108  >  indicating  the  highest  notch  sensitivity  for  that 

angle.  Similar  results  were  obtained  for  a  in  the  WEK-fracture  model.  No  such 

c 

correlation  could  be  established  between  the  PWG-fracture  model  parameters  (i.e. 

ac  and  a^)  and  the  angle  9  of  the  four  laminates.  The  effect  of  the  angle  9  on 

the  relative  notch  sensitivity  parameter,  R  ,  is  shown  in  Figures  109-110, 

ns 

indicating  that  it  is  independent  of  ply  orientation  for  most  cases.  Similarly, 
no  correlation  between  the  ML-fracture  model  parameters,  i.e.  n  and _Hr/oQ  and 
the  angle  9  could  be  established,  Figures  111-114. 
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STRENGTH  RftTIO 


HALF  CRACK  LENGTH  CMfl) 

Figure  101.  Notched  etrengch  venue  creek  length  prediction*  (UN- free  cure  aodel,  '’average  etreea'' 
criterion)  for  boron/alualmm  [ 0 1 a  laalnate*. 


H!ur.  102.  Settle  •er«n|Ch  dl.tubuclon  ,[  U«i,nu,l  5.6  ,(1  41„,t.r  bo„„  (lb.t,  ,;0l, 
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THETft  <±8> 

l05.  CxperlMncal  value.  ol  cHaracttrlatlc  dlMntlon,  a  (W-fractur.  uod.l  "mrip  .cr...  1 
criterion)  > ■  a  function  of  angle  for  graphita/epoxy  lt«lf  lanlnaca.  containing  circular 
holta.  Individual  daca  arc  given  In  Tapia  11  and  Appendix  g. 


THET*  <(J/ie) 

Flgura  106.  Ixparlnaotal  valuaa  of  charac tarlatlc  dinanalon,  a0,  (W-f racturi  nodal  "average  atraaa" 

erlear’on)  aa  a  function  of  angla  for  graphlta/apoxy  laalnataa  containing  elrcula- 

holaa.  Individual  data  ara  glvan  in  Tab la  11  and  Appendix  I. 
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O  Mf.  [56] 
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Figure  107.  Experimental  values  of  characteristic  dimension,  a  .  (WN-f ractura  aodal  **avaraga  scrass" 

criterion)  as  a  function  of  angle  for  graphlte/epoxy  [t5/0)_  laminates  containing  circular 
holes.  Individual  data  are  given  la  Table  11  and  Appendix  I. 


a  *  10  30  40  30  60  79  00  30 

THETA  C+e/Os-d) 

Figure  100.  Experimental  values  of  characteristic  dimension,  ep,  (WH-fracture  aodal.  "average  stress" 

criterion)  as  a  function  of  angle  for  graphlte/eposy  I-H/0/-4)  laminates  containing  circular 
holes.  Individual  data  are  given  in  Table  11  and  Appendix  B. 
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109  Enparlnancal  data  ot  ralaclva  nntci  nnnalrlvlty,  .  (MJG-Itactur.  nodal)  aa  a  function 
o I  anfla  for  iraphlca/apoap  (0/U|,  laalnacaa  concalnlni  circular  holaa.  Individual  data 
ara  (Ivan  In  Tabla  13  and  Appandla  0. 


T  H  £  T  ft  C  ±  9  /  0  ) 

Fi(ur«  110.  trparlaancal  data  of  ralaclva  notch  aanalelvlcy.  .  (FU£-fractufa  nodal)  at  a  function 
of  angla  for  graph  lea /apoiry  (*9/01  laalnacaa  containing  circular  holaa.  Individual  data 
ara  glvon  in  Tabla  13  and  Appaadbi  p, 
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a  l  (J  £0  3d  14  34  64  70  84  30 


THETA  O'*  6) 

Figure  113-  Experimental  valuta  of  fracture  toughness,  Re/a0»  (ML-freeture  eodel)  «•  *  function  of  angle 
foe  graphic*/ epoxy  (O/UJ  laxlnats*  containing  circular  holes.  Individual  data  ara  given 
In  Table  14  and  Appendix  I. 
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rigur*  11*.  Experimental  values  of  fracture  toughness,  Ke/c0,  (ML-frectur#  eodel)  as  a  function  of  angle 
for  graphite /epoxy  (♦d/0/-4|  laminates  containing  circular  holes.  Individual  data  are  given 
In  Table  U  and  Appendix  t. 
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CORRELATION  BETWEEN  NOTCH  SENSITIVITY  AND  FRACTURE  MODEL  PARAMETERS:  The 
different  fracture  models  reviewed  in  this  report  utilize  a  variety  of  para¬ 
meters  associated  with  each  model  and  which  must  be  determined  experimentally. 

In  all  the  models  it  is  assumed  that  these  parameters  are  constant,  i.e. 
independent  of  notch  size,  and  are  to  be  considered  as  material  parameters. 

It  is  of  interest  to  identify  which  of  these  parameters  can  be  correlc.ted  to 
the  notch  sensitivity  of  composite  laminates  containing  circular  holes  and 
straight  cracks.  For  this  purpose,  the  ratio  of  notched  to  unnotched  strength 
(a”/aQ)  at  a  specific  notch  size  has  been  recorded  from  the  data  such  as  those 
listed  in  Table  11,  and  it  has  been  correlated  with  the  various  parameters 

associated  with  the  different  fracture  models.  The  notch  sizes  for  which  o.,/o 

N  o 

was  recorded  are:  2R/W  -  0.25  for  graphite/epoxy  laminates  containing  holes; 
and  2c/W  -  0.30  for  boron/aluminum  laminates  containing  straight  cracks.  The 
correlation  between  the  notched  strength  and  the  various  parameters  is  shown 
in  Figures  115  to  134.  Each  point  in  these  Figures  represents  a  set  of  notched 
strength  data  as  shown  for  example  in  Table  11.  Consequently,  these  Figures 
combined  all  the  sets  of  notched  strength  data  reviewed,  irrespective  of 
material  system,  constituents  properties,  fabrication  procedures,  laminate 
configuration,  stacking  sequence,  specimen  geometry,  loading  procedure,  etc. 
Thus,  significant  scatter  exists  in  the  values  of  all  the  parameters  and  it  is 
emphasized  that  conclusions  regarding  the  presence  or  absence  of  correlation  can 
only  be  qualitative,  i.e.  restricted  to  the  specific  specimen  geometry,  loading 
function,  etc.  which  are  analyzed  in  this  study.  It  should  be  noted  that 
although  the  results  shown  in  Figures  115  to  134  are  for  a  specific  notch  size 
(of  2R/W  •  0.25  and  2c/W  -  0.3),  similar  Figures  and  correlations  are  obtained 
for  all  other  notch  sizes. 


107 


A. 


A 


a.  WEK-fracture  model  (Figures  115-116):  A  correlation  exists  between 

notch  sensitivity  and  a  .  The  smaller  a_  ia,  the  more  notch  sensitive  the 
c  c 

subject  material.  This  is  to  be  expected  from  Eq.  (11).  The  few  exceptions 
are  for  the  data  depicted  from  Refs. [204]  and  [208]  as  Indicated  in  Figure  116. 
The  values  of  a^  for  the  laminates  tested  in  [204]  are  based  on  approximately 
20  notched  strength  data  for  each  laminate.  However,  the  data  were  obtained 
from  three  different  specimen  widths,  19.1  mm,  50.8  mm  and  101.6  mm.  It  is 
possible  therefore  that  lumping  all  notched  strength  data  of  the  different 
laminates  widths  is  the  cause  for  this  discrepancy.  The  data  given  in  [208] 
are  for  4.2  mil  borsic/alumlnum  which  is  a  different  system  from  the  5.6  mil 
boron/aluminum  laminates  tested  in  most  other  sAts  of  data. 

b.  WN-fracture  model  (Figures  117-120):  A  correlation  exists  between 
notch  sensitivity  and  dQ,  Figures  117-118,  and  aQ,  Figures  119-120.  The 
smaller  the  characteristic  distances  are,  the  more  notch  sensitive  the  subject 
material,  as  expected  from  Eqs.  (23,27)  and  (25,31),  respectively. 

c.  K-fracture  model  (no  Figures):  No  correlation  could  be  found  between 
the  notch  sensitivity  and  kQ. 

d.  PWG-Fracture  model  (Figures  121-130):  No  correlation  could  be  found 
between  notch  sensitivity  and  the  various  parameters,  i.e.  exponential  para¬ 
meters,  m.  Figures  121-122,  the  notch  sensitivity  factors  C  and  K,  Figures  123- 
124,  the  notch  shift  parameters,  ac  snd  a^  Figures  125-126.  and  the  notch 
sensitivity  factors,  aB>  Figures  127-128.  The  only  exception  is  in  regard  to 
the  relative  notch  sensitivities  Rn#  and  a^.  Figures  129-130,  which  show  a 
clear  correlation  between  their  values  and  the  notch  sensitivity.  The  larger 
Rns  and  *Km  *^e,  the  more  notch  sensitive  the  subject  material  is.  This 
correlation  is  expected  considering  the  definition  of  these  two  parameters 
given  previously. 
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Figure  121.  Correlation  between  notched  etrengch  ratio  (for  21/W  •  0.23)  aod  exponential  ptriMttr,  a. 
(FVG-frmcture  nodal)  for  graphite/ epoxy  laalnatea  containing  circular  holaa. 
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STRENGTH  RATIO 


NOTCH  S  E  N  3 .  FACTOR  (PUG) 

Figure  123.  Correlation  baeveen  notched  terangch  ratio  (for  2*/W  •  0.23)  and  notch  aanaitlvlty  factor.  C 
(PUG-f ractura  nodal)  for  graphlta/apoxy  lanlnatas  containing  circular  holaa. 


NOTCH  SENS.  FACTOR  CPUS? 

Flguta  124.  Correlation  batwaan  notchad  atrangch  ratio  (for  2c/W  -  0.3)  and  noeeh  aanaitlvlty  facte 
K.  (Pvc-f ractura  nodal)  for  boron/alualnun  lanlnataa  containing  atralght  cracks. 
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Fljura  125.  Correlation  batvaaa  notched  acrangch  redo  (for  2*/W  -  0.23)  end  ahift  piristcir,  a 
(FWC-ftaccur*  nodal)  (or  |Taphlta/apox?  laminata*  containing  circular  holaa.  c 


Fltura  126 •  Corralaeion  between  notched  arrangth  ratio  (for  2e/«  •  0.5)  and  ahift  parameter,  a- 
fraceura  nodal)  for  kom/aluaiina  laalnaeaa  concatnlni  atrnl|ht  crack*.  * 
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The  results  shown  In  Figures  121-128  Indicate  that  no  correlation  could 
be  established  between  the  various  parameter  (m,C,K,am,ac,  and  ax)  and  the  notch  sensi¬ 
tivity  C  and  K  of  composite  laminates  containing  either  circular  holes  or 
straight  cracks.  Although  the  parameteric  analysis  shown  in  Figures  33  and  34 
for  holes  and  50  and  51  for  cracks  indicate  that  there  is  a  direct  correlation 
between  m,  C,  K,  and  notch  sensitivity,  i.e.  larger  values  evidence  higher 
notch  sensitivity,  it  should  be  recalled  that  these  parametric  analyses  were 
performed  where  only  one  of  the  parameters  was  a  variable  while  the  other 
parameter  was  held  constant.  As  previously  discussed,  each  laminate  is 
characterized  by  a  different  pair  of  m  ard  C  or  m  and  K  (for  holes  and  cracks, 
respectively).  For  all  practical  purposes,  no  correlation  can  be  made  between 
their  specific  values  for  a  given  composite  laminate  and  the  notch  sensitivity 
of  chat  laminate.  Similarly,  no  such  correlation  can  be  made  with  any  of  the 
other  parameters  associated  with  the  PWG-fracture  model,  except  the  relative 
notch  sensitivities  Rng  and  a^. 

e.  Mi-Fracture  model  (Figures  131-134):  No  correlation  could  be  found 
between  the  exponential  parameter,  n,  and  notch  sensitivity,  Figures  131-132. 

For  laminates  containing  circular  holes,  most  values  are  between  0.1  to  0.4; 
while  for  laminates  containing  straight  cracks,  Che  values  range  between  0.2  to 
0.45,  all  with  very  large  scatter.  Neither  could  any  correlation  be  found 
between  the  composite  fracture  toughness,  H_/o  ,  and  notch  sensitivity,  Figures 
133-134.  The  values  of  Hc/oo  range  between  0.8  to  1.2  for  laminates  containing 
both  circular  holes  and  straight  cracks. 

f.  Effect  of  Model  Parameters  Selection  on  Predictions  (Figure  135):  Finally, 
many  studies  address  whether  the  constants  associated  with  the  different  fracture 
models  are  independent  of  either  material  system  or  laminate  configuration.  The 
results  shown  previously  clearly  indicate  that  all  the  constants  strongly  depend 
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figure  131.  Correlation  between  notched  •erengch  retlo  (for  2t/W  •  0.23)  end  exponential  paraaeter,  n, 
(Mt-f racture  aodel)  for  graphlte/epoxy  laalaatee  containing  circular  holes. 


132.  Correlation  between  notched  strength  ration  (for  2c /V  •  0.3)  and  exponential  paraaeter.  n, 
(ML-fr secure  aodel)  for  boron/elualnua  laainacee  containing  straight  cracks. 
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Figure  133. 


Correlation  between  notched  acrength  ratio  ((or  2K/U  •  0.23)  aod  eoapoalce  fracture 
HC/oo»  (ML-(raccure  aodel)  for  graphlce/epoxy  laalnatea  containing  circular  hole*. 


toujhne: 


Figure  1J4.  Correlation  between  notched  atrengch  ratio  (for  >c/U  •  0.3)  end  composite  fracture  toughneaa, 
*C/0o’  •©<•!)  for  boron/altaalnua  laninatea  containing  atralght  cracka. 
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on  laminate  configuration  as  well  as  other  factors  such  as  fabrication  procedures 

fibers'  properties,  etc.  Consequently,  the  applicability  of  a  given  fracture 

model  parameter  is  restricted  to  the  specific  set  of  data  for  which  it  has  been 

determined.  Figure  135  shows  a  representative  notched  strength  data  of  a  [02/i45] 

laminate  compared  with  the  WN-fracture  model  using  relatively  small  and  large 

values  of  d  and  a  within  the  range  of  values  listed  previously.  The  character- 
o  o 

istic  distances,  d  and  a  ,  of  the  [0,/±451  laminate  have  intermediate  values. 

O  O  Z  s 

The  comparison  clearly  indicates  the  degree  of  error  which  can  result  by  choosing 
arbitrary  values  for  do  and  aQ  and  applying  them  to  all  laminate  configurations 
of  a  given  material  system. 


Figure  133.  Effect  of  th#  character la tic  dimensions  (Wtt-frsctur*  model)  on  notch**  strength  of  boron/ 
•lunlnta  laalnstaa  containing  straight  crack*. 
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Symbol  Definition 


A 


ij 


a 

a 


c 


a 

cm 

aKra 

aKra 

a 

m 


Orthotuplc  in-plane  stiffnesses 

Length  of  intense  energy  region  (WEK-f racture  model). 

Crack  tip  damage  zone  size  at  failure  (WEK-f racture  model). 

Radius  shift  parameter  for  notch  sensitivity  factor  (PWG-fracture 
model) . 

Generalized  radius  shift  parameter  (PWG-fracture  model). 

Crack-length  shift  parameter  for  notch  sensitivity  factor  (PWG- 
fracture  model). 

Generalized  crack-length  shift  parameter  (PWG-fracture  model). 

Generalized  notch  sensitivity  factor  for  cracks  (PWG-fracture 
model) . 

Radius  shift  factor  for  exponential  parameter  (PWG-fracture  model, 


Characteristic  dimensions  (WN-f racture  model;  average  stress 
criterion) . 


C 

C* 


Radius  notch  sensitivity  factor  (PWG-fracture  model). 

Radius  notch  sensitivy  factor  on  master  curve  (PWG-fracture 
model) . 

Critical  hole  notch  sensitivity  factor  (PWG-fracture  model). 


c  Half  crack  length 

c'  Shifted  half  crack  length  by  a^  (PWG-fracture  model). 

c*  shifted  half  crack  length  by  a^  (PWG-fracture  model). 

c  Corrected  half  crack  length  (PWG-fracture  model) . 

c 

c*  Critical  corrected  half  crack  length  (PWG-fracture  model). 

+  Efforts  were  made  in  this  report  to  use  the  same  nomenclature  which 
appeared  in  the  papers  reviewed.  Therefore,  in  a  few  cases  the  same 
symbol  may  have  different  definitions  or  the  same  definition  may 
have  different  symbols. 
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Q- 1 


c^  Notch  Insensitive  half  crack  length  (PWG-fracture  model), 

c*  Critical  notch  insensitive  half  crack  length  (PWG-fracture  model). 


Reference  half  crack  length  (PWG-f racture  model). 
Effective  characteristic  dimension  (WN-f racture  model). 


d* 

o 


E 

EL 

et 

Euy 


Characteristic  dimension  (WN-fracture  model,  point  stress 
criterion) . 

Critical  characteristic  dimension  for  holes  and  cracks  (PWG- 
fracture  model). 

Characteristic  distance  for  "general  fracture-toughness  parameter'! 


* 


Young's  modulus 
Longitudinal  lamina  stiffness 
Transverse  lamina  stiffness 

Ultimate  (secant)  longitudinal  laminate  stiffness  (PS-fracture  mode)) 
Effective  transverse  laminate  stiffness 
Effective  longitudinal  laminate  stiffness 


tu 


Ultimate  tensile  strength  of  unnotched  laminate  (PS-fracture 
model) . 


G 

C, 


-LT 


yx 


Shear  modulus 

Energy  release  rate  (mode  I) 

Lamina  in-plane  shear  modulus 
Effective  in-plane  laminate  shear  modulus 

Composite  fracture  toughness  (ML-f racture  model). 


K* 


Crack  notch  sensitivity  factor  (PWG-fracture  model). 

Mode  I  stress-intensity  factor  (PS-fracture  model). 

Crack  notch  sensitivity  factor  on  master  curve  (PWG-fracture 
model) . 
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n 


Kt 


x IC 


kq 

KQe 


eQ 


eQe 


Kt 

CO 

*T 


Critical  crack  notch  sensitivity  factor  (PWG-fracture  model). 
Mode  I  stress-intensity  factor 
Mode  I  fracture  toughness 

Critical  stress-intensity  factor,  fracture  toughness 
Modified  critical  stress-intensity  factor  (WN-fracture  model). 
Elastic  critical  stress-intensity  factor  (PS-fracture  model). 
Strain-intensity  factor  (PS-fracture  model). 

Critical  strain-intensity  factor  (PS-fracture  model). 

Elastic  critical  strain-intensity  factor  (PS-fracture  model). 
Stress  concentration  factor  for  a  finite  width  plate 
Stress  concentration  factor  for  an  infinite  plate 
Proportionality  constant  (K-fracture  model). 

Length  of  specimen 


m  Exponential  parameter  (PWG-fracture  model), 

m*  Exponential  parameter  on  master  curve  (PWG-fracture  model). 


mi 


n 


n 


Critical  exponential  parameter  (P.WG-f racture  model). 

Order  of  singularity,  exponential  parameter  (ML-fracture  model). 


Q 


c 


General  fracture-toughness  parameter  (PS-fracture  model). 


R 

R* 


K 

Rt 


Hole  radius 

Shifted  hole  radius  by  a  (PWG-fracture  model). 

cn 

Critical  corrected  hole  radius  (PWG-fracture  model). 

Critical  notch  Insensitive  hole  radius  (PWG-fracture  model). 

Relative  notch  sensitivity  factor  for  holes  and  cracks  (PWC-f ractuie 
model) . 


A 
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Critical  relative  notch  sensitivity  factor  for  holes  (PWG- 
fracture  model). 

Reference  hole  radius  (PWG-f racture  model). 

Polar  coordinates 

Par-field  applied  stress  (PS-fracture  model). 

Non-dimensional  strength  at  a  given  notch  radius  (PWG-f racture 
model ) . 

Notched  laminate  strength  (PS-fracture  model). 

Transformation  matrix  (PS-fracture  model). 

Laminate  thickness 

Coordinate  measured  from  center  of  notch  perpendicular  to 
applied  longitudinal  stress 

Finite  width  correction  factor 

Coordinate  measured  from  center  of  notch  directed  in  the  loading 
direction 

Width  of  specimen 


Fiber  orientation  angle 
«  1/ (m*-l)  (PWG-f racture  model). 

“  (m*-m) /(m*-l)  (PWG-f racture  model). 

Matrix  of  constituents  properties  (PS-fracture  model). 
Shear  strain  (in  laminate  reference  system). 

Shear  strain  (in  laminae  reference  system). 

Normal  strains  (in  laminate  reference  system). 

Normal  strains  (in  laminae  reference  system). 

Far  field  strain  of  notched  laminate  (PS-fracture  model). 


Far  field  strain  of  notched  laminate  at  failure  (PS-fracture 
model ) . 

Ultimate  tensile  strain  of  unnotched  laminate  (PS-fracture  model). 

Ultimate  tensile  strain  of  fibers  (PS-fracture  model). 

Material  constants  (PS-fracture  model). 

Angle  (Figure  76) 

-  (1  +  Rm_1  R'1  C-1]-1  (PWC-fracture  model). 

o 

■  1  on  master  curve  defined  by  C*  and  m*  (PWG-fracture  model). 

-  (1  +  cm  ^cq  m  K  *)  ^  (PWG-fracture  model). 

-  X^on  master  curve  defined  by  K*  and  m*  (PWG-fracture  model  ). 

-  (1+(R*) (mi-1)  (R  )“mi  (C*)-1]-1  (modified  PWG-fracture  model). 

C  Ox 

Shear  modulus  of  matrix 

shear  modulus  of  filament 

Poisson's  ratio 

Poisson's  ratio  of  matrix 

Poisson's  ratio  of  filament 

Effective  laminate  major  Poisson's  ratio 

“  R/ (R+do)  (WN-fracture  model,  "point  stress"  criterion  for  holes). 

“  R/ (R+ao)  (WN-fracture  model,  "average  stress"criterion  for  holes). 

•  c/(c+dQ)  (WN-fracture  model,  "point  stress"  criterion  for  cracks). 

-  c/(c+aQ)  (WN-fracture  model,  "average  stress"  criterion  for  cracks). 

Crack  tip  damage  zone  size  in  terms  of  stress  (PS-fracture  model). 

Critical  crack  tip  damage  zone  size  in  terms  of  stress  (PS-fracture 
model) . 

Crack  tip  damage  zone  size  in  terms  of  strain  (PS-fracture  model). 


Critical  crack  tip  damage  zone  size  in  terms  of  strain  (PS- 
fracture  model). 

Axial  stress 

Far-field  longitudinal  stress 

Notched  tensile  strength  of  a  finite  width  plate 
Notched  tensile  strength  of  an  infinite  plate 
Unnotched  strength 

Normal  stress  component  in  the  y-direction  (loading  direction) 
Stress  distributions  in  the  laminate  reference  system  (Figure  53) 

Abbreviations 
CRIT.  Critical 

FVC  Finite  Width  Correction  Factor 

K  Karlak  (-fracture  model) 

ML  Mar-Lin  (-fracture  model) 

PS  Poe-Sova  (-fracture  model) 

PWG  Pipes-Wetherhold-Cillespie  (-fracture  model) 

SCF  Stress  Concentration  Factor 

SIF  Stress  Intensity  Factor 

WEK  Waddoups-Eisenmann-Kaminski  (-fracture  model) 

WN  Whltney-Nuisraer  (-fracture  model) 


Units 

GPa 

Glga  Pascal,  10^  N/m2 

in 

inch  («  25.4  mm) 

KN 

1000  Newton 

lb 

Pounds  (*  4.4  3  Newton) 

ksl 

1000  lb /in 2  (»  6.894  MP 

126 


millimeter 


N  Newton  (■  0.2248  lb) 

psi  lb/ln^ 
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SUMMARY 

In  order  to  determine  the  nature  of  failure  mechanisms  a  number  of  fatigue  tests  were  per¬ 
formed.  The  test  specimens  partly  have  artificial  delaminations  between  different  layers 
of  the  multidirectional  laminates  made  from  T300/91i*C  prepregs .  For  better  understanding 
of  the  strength  degradation  in  fatigue  a  damage  model,  based  on  the  delamination  propagation, 
starting  from  the  free  eages  between  the  plies  of  a  multidirectional  laminate, has  been 
developed.  These  defects  propagate  due  to  interlaminar  stresses  up  to  an  area,  which  is 
critical  in  the  case  of  tension-compression  fatigue  against  buckling  or  shearing  of  parts 
of  the  delaminated  test  specimen. 


LIST  OF  SYMBOLS 

{ a )  Extensional  stiffness  matrix,  KN/mm 
An  Stiffness  of  the  laminate  in 
x-direction,  KN/mm 
a  Strip  delamination  size,  crack 
length,  mm 

a0  Platie  zone  at  the  crack  tip,  mm 
do/dN  Delamination  growth  rate,  mm/cycle 
IB)  Coupling  stiffness  matrix,  KN 
b  Width  of  the  specimen,  mm 
c,  Displacement  at  the  edge  of  the 
specimen  in  z-direction,  mm 
Cm  Measured  value  of  the  crack  opening 
displacement  at  the  edge  of  the 
specimen,  mm 

(Di  Bending  stiffness  matrix,  KN/mm 
d  Laminate  thickness,  mm 
Ei  Modulus  of  a  lamina  in  fiber 
direction,  KN  mm2 

Et  Modulus  of  a  lamina  in  transverse 
direction,  KN/mm7 
GLt  Lamina  shear  modulus,  KN/mm2 
h  Lamina  thickness,  mm 
u  Plate  curvature,  1/mm 
M  Resultand  moment,  acting  on  the 
laminate  Nmm/mm 

n  Resultand  force,  acting  on  the 
laminate,  N/mm 
N  Number  of  load  cycles 
R  Stress  ratio,  R  -6 l  ICf 
w,  Displacement  In  z-direction,  mm 
x. y,z  Cartesian  coordinates 


&  Phase  lag 

Midplane  strain 
vtJ1vLT  Poission’s  ratio 

gnig(2i  Curvature  of  delaminated  parts  (1) 
and  (2 )  ,  1/mm 
6 a  Stress  amplitude,  N/mm? 

0,,  Calculated  lamina  normal  stress, 
i=  1,2,3,  N/mm2 

e,,  Calculated  lamina  shear  stress 
i  ,  j  =  1,2,5,  N/mm2 

e*0y6i  Laminate  normal  stress  in  x-,  y-, 
and  z-direction 
Yield  stress,  N/mm2 
Cl  Axial  strength,  N/mm2 

subscript 

c  Critical 
d  Damaged 

f  Failure  by  static  loading 
F  Fatigue  failure 
L  Longitudinal 
T  Transverse 
R  Residual 

superscript 

u.l  upper,  lower 
y  Yield 

(1)12)  Portion  1,  2 


1  INTRODUCTION 


In  the  failure  analysis  of  composite  laminates,  one  of  the  most  serious  problem  has  been 
the  propagation  of  interlaminar  cracking,  commonly  known  as  delamination.  Delaminations 
may  be  formed  during  manufacture  due  to  incomplete  curing  or  the  introduction  of  foreign 
particle;  they  may  result  from  impact  damage;  and  they  may  result  from  interlaminar 
stresses  existing  at  discontinuities  or  at  the  stress-free  edges  of  a  loaded  composite 
structure.  This  mode  of  failure  is  a  major  cause  for  the  deterioration  of  laminate 
structural  properties,  including  its  strength,  stiffness,  reliability  and  durability. 
Furthermore,  delaminations  may  grow  under  increasing  load  or  cyclic  loading.  Delamination 
growth  redistributes  the  stress  in  plies  of  a  laminate,  and  may  influence  residual  stiffness, 
residual  strength  and  fatigue  life  expecially  under  compressive  fatigue  loads.  Hence,  a 
fatigue  analysis  for  composite  laminates  should  take  in  to  account  the  presence  and  the 
growth  of  de laminat ions . 


The  exact  formation  and  growth  mechanisms  of  interlaminar  cracking  in  laminates  are  not 
well  unterstood.  The  general  belief  is  that  a  certain  distribution  of  small  interface 
flaws  with  a  size  in  the  order  of  the  fiber d iameter  exist  in  laminates  prior  to  loading 
Under  a  critical  loading  condition,  which  include  thermal  (curing)  stresses,  swelling 
stresses  of  moist  laminates  and  external  loadings,  some  flaws  would  grow  and  coalesce 
with  each  other,  forming  a  single  crack  of  macroscopic  proportion.  Such  an  event  would 
constitute  the  onset,  or  the  initiation  of  the  macro-crack. 


JL 


T ••  fiber  composite  laminate  made  of  {0/-^5/90]s  stacking  sequence  the  damage  develop- 
i:  under  tension  load  usually  starts  with  the  appearence  of  matrix  cracks  in  the  off- 

-plies  transverse  to  the  load  direction.  As  tension  load  or  load  cycles 
:>  ,the  matrix  cracks  in  the  off-axis  plies  grow  up  to  the  adjacent  layer  and  a 

, r*  f  ri  ?,t.  ic  pattern  of  matrix  cracks  forms  as  shown  in  Fig.  1.,  see  intern  1.  A  stable 
i  vn  of  regularly  spaced  matrix  cracks  develops  in  each  off-axis  ply  predictable  by 
a  simple  one  dimensional  model  [l]  .  Riefsnider  named  this  pattern  the  characteristic 

•  i-t-n-*  state  (CPS)  for  matrix  cracking.  Where  the  matrix  cracks  terminate,  local  interface 
•'ri’ks  frequently  form  at  the  ply  boundaries.  For  the  purpose  of  illustration  a  [  0 , 9  0 ]„ 
lirinate  is  presented  in  Fig.  1.  During  the  cyclic  loading  localized  delaminations 

">'■>  ]  ps  in  the  boundary  between  the  0°-and  90°- layer  (see  item  No.  2)  mainly  at  the  edges 
f  !  Mi-*  specimen,  but  in  the  interior  also  (see  for  example  item  No.  5)  (2jt(3j 

barl.v  investigators  focused  mainly  on  experimental  and  analytical  work  of  free  edge  induced 
1  nr. irat ion  in  laboratory  test  specimens.  A  bibliography  of  this  work  is  contained  in 
?  ‘i.  [;ij  .  After  one  has  an  understanding  of  the  stress  distribution  near  a  free  edge  by 
f Moment  calculations,  the  next  step  in  treating  the  delamination  problem  is  to  ex¬ 
plain  Me*  mechanism  and  the  criterion  for  the  edge  delamination.  One  of  the  most  promising 

♦  rhri  iqurp  for  characterizing  delamination  growth  is  based  on  the  strain  energy  release 
r-M*  r  during  delamination  growth.  Measured  critical  values  of  G  have  been  used  to  pre- 
d;  *t  the  onset  and  the  growth  of  edge  de laminations  in  composite  laminates  [5]  ,  (6). 

Tr  M--  present  study  a  technique  was  developed  to  characterize  the  onset  and  growth  of 
: ions  in  composite  laminates.  First,  the  damage  that  developed  in  unnotched 
i  1  :  -  /-fJ^/O/oo  j  graphite-epoxy  laminates  under  static  and  cyclic  tension  and  com- 
a.ling  was  determined  by  nondestructive  test  methods.  During  test  loading  crack 
i  r:  placement  was  monitored  to  relate  delamination-crack  opening  with  delamination 
.  •  .  ’.'"Xt  ,  stress  distribution  generated  from  a  finite  element  analysis  was  correlated 
with  ;•  -->rvc.'i  damage.  The  resulting  test  data  and  analysis  were  used  to  derive  a  closed- 
f-rm  equation  for  the  edge  crack  opening  assoziated  with  the  delamination  growth. Finally 
n  '■  "rrclarion  between  delamination  size  and  residual  compression  strength  of  fatigued 
n  r'inenf!  was  established. 

2 .  MATERIAL  AND  TEST  EQUIPMENT 

Laminates  of  [  0p/  +  *45/02/-M5/0/90  ]  stacking  sequence  were  fabricated  from  91J1C/T?00B 
rr'5r h i ! e/epoxy  ^repreg^tape  in  a  computer-controlled  autoclave  according  to  the  manufac¬ 
turers  recommended  cure  procedure.  Nominal  fiber  volume  of  the  laminate  was  60j +  2%. 

The  rix* eon-ply  laminate  had  an  average  ply  thickness  of  0,125  mm.  The  380  by  380  mm 
p- were  bonded  with  fiberglass  tabs  of  1,0  mm  thickness  for  the  clamping  reinforcement, 
'if:' •  .c  pannls  were  cut  by  a  diamond  saw  to  specimens  of  different  size  depicted  in  Fig. 2. 

. >f  the  +  ype  II  specimens  had  teflon  tabs  or  an  inflation  agent  embeded  between  the 
t'*:;:;. nr ios  of  different  layers.  The  coupons  were  stored  and  tested  under  laboratory 
'•nui  i  +  ions  (21+  1°C  and  50+  5X  relative  humidity)  3  to  6  months  arter  fabricat  ion  /  Some 
cr  i :r.»’nr;  were  stored  in  a  climatic  chamber  for  different  times  to  obtain  different 
n:  : M  >re  distributions  and  moisture  contents. 

Al!  »«‘sV.s  were  conducted  on  a  closed-loop  Schenck  hydraulic  testing  machine.  For  displace- 
measurement s  strain-gauge  transducers  (SOT)  were  mounted  on  the  specimens.  To  prevent 
si i ,  a  fast  drying  glue  was  applied  on  the  transducer  mounts  where  touched  the 

i  men . 


:•■•**  m*.'isurement  of  the  transverse  contraction  and  transverse  crack  opening  displacement 
L1!  on  some  specimens  we  used  the  DFVLR-MDR-Transducer  (magnetic  field  depending 
resistor)  (7J  .  The  mounting  of  the  transducer  on  the  test  specimen  is  shown  in  Fig.  3. 

After  o.ich  step-loading  or  after  certain  load  cycles,  the  specimen  was  removed  from  the 
t s t p r  and  delamination  size  measurements  (C-scan)  were  made  by  an  ultrasonic  testing 
m it  with  a  narrow  waveband  emitter  [8]  .  For  an  exact  ajustment  of  the  specimens  in 

+  ■'••■'Sting  machine,  the  end  of  the  specimens  and  the  clamping  device  contain  fitting 

Vice,  p-yo-penetrant-enhanced  radiography  was  not  used  because  of  the  unknown  influence 
» he  penetrant  fluid  on  the  crack  propagation.  A  Fhilips  SEM  505  scanning  eletron 
microscope  was  used  to  study  and  document  the  topographic  features  of  the  delaminated 
specimen  fracture  surfaces. [ 9 J  Prior  SEM  examination,  gold  was  sputtered  onto  the  fracture 
surface  to  obtain  optimum  resolution  of  the  topografic  features  and  to  minimize  static 
chirring  by  the  SEM  beam. 

OBSERVATION  OF  DELAMINATION  GROWTH 

Puri  nr,  both  quasi-static  tension  tests  and  constant  amplitude  tension-tension  fatigue 

on  [  Op/  +  1i5/02/-45/O/90]  laminates  the  same  type  of  d  image  developed.  First  a  few 
icMated  transverse  cracks  forffied  in  the  90°-  plies.  As  these  cracks  grew  in  transverse 
Erection  the  crack  spacing  decreased  by  forming  new  transverse  cracks.  These  were 
1  cl  lowed  almost  by  the  onset  of  small  delaminations  along  the  edge  at  the  transverse  crack 
tip  as  seen  in  Fig.  H.  The  length  of  the  small  delaminations  grew  up  to  the  interconnection 
with  an  other  delamination  during  increasing  static  load  or  load  cycles.  Finally  the 
^lamination  grew  much  more  rapidly  along  the  length  of  the  specimen.  In  all  cases 
Jo  laminar  10ns  on  both  sides  of  the  specimen  extended  along  the  entire  specimen  length 
c*»tween  the  grips.  0 


5-3 


t.  he 

\j  1 1 


i  l  lust  ra* c  t  tv  -  ;-r  h  of  t1-3  lo 1 nminat ion  a  few  ultrasonic  C-.-’oans  were-  made  'u  type  II 

•  specimen  -  1  r.  I  r.r  circles  iir.pl  ants  between  both  of  the  9-1  -d**gro»*  layers 

?  ho  cent  ral  i  ?f.  i i  :n .  Fir.  ^  shows  U  C-scans  of  a  stepwise  tension  loaded  specimen . 

rirth  har.  i  C-acan  was  made  after  final  failure.  This  specimen  frae*urod  on  the  lower 
>  with  a  V- oii-are  1  do  1  aminot.  i  on  inside  the  measurement  area.  Delamir.  ition  started  at  a 

;■!  at  '■u*  '•  r-.:  at  t  -Mi  ^  lros  In  the  boundary  between  the  Q-degre'*  and  y  -degree 

■  ;»::d  rro>;  r  \  ‘-'•■•aosr*  ■'■V  f  he  clamping.  effect  the  delaminat  i  on  sbap«  was  twi  rt  ed  .  There 

nc  de ’ aminat i on  growl h  during  tension  loading  starting  from  the  implants.  In  the  same 
a  test  specimen  with  an  artificial  blowing  agent  delamination  degraiar  .■  1  during  static 
sion  leading.  Fir.  h.  The  internal  delaminat_on  started  from  the  sharp  crack  tips  of 
artificial  blowing  agent  de lamination  and  grew  very  slowly.  The  specimen  failed  at  a 
imat*  strength  of  lOdO  N/mm?  after  30?  of  90-degree  layer  were  delaminated.  There  is 
r  i  ;7  • .  i  fie  art  d  i  f  f*:  ••«*»«.>>  I  e*  w**er  the  ultimate  strength  of  artificial  delaminated  cod  no 
aminatel  specimens.  The  load  at  edge  delamination  onset  as  well  as  the  load  at  dela- 
atMr  rb«  nr*1  ire  length  and  the  ultimate  load  depend  upon  the  moisture  content  and 
.  i  ure  profile  of  t.  hrt  specimen.  Fig.  7.  The  tests  were  conducted  in  all  cases  <-.t  a 


t  •  r “fi v  ire  >f  -'1 


1  durit'.r  compressive  loading  delamination  started  from  the  implants  as 
i- 7 ,  "  !••'!  am irat  5  n  was  observed  up  to  the  failure  lead  at  the  edges  r  the 

n  h i  Mr  ill  growing  of  the  central  delamination  was  observed. 


<  Jiffcrences  in  delamination  growth  between*  static  tension  t*=st  on  one 

)  i  r  r  ion  fat  i  rue  test  on  the  other  hand.  Compared  to  those  behavior  in 

rr  r  I  -  •:  fatigue  ♦,  o§t-  delaminat  ions  appeared  between  both  the  Op  degree  and 
•  m-  •  ■« t-r  >  >•••  1  i  at-  between  the  4 5 -degree  and  0  degree  layer’s.  Fig.  9  shews  the 

c-  •  nV  i  •  u  m  of  a  spec  imen ,  which  was  loaded  in  tens  i  on-compress  i  on  by  a  stress 

•  c.' •' ’■  :iJ  with  a  stress  ratio  R  =  -1.  After  30.000  cycles  at  the  urp^r 
•  ''  o;  -  •  n/sim*  the  delaminated  crack  between  the  0-  degree  and  9°  degree 

-r  w.-.  !,-■  f'rick  opening  size  is  about  0,1  mm.  Fig.  9a .  The  sai?«*  de  lamina  tod 

is  r b-’wri  in  Fig.  To  aft^r  ur.laoding  to  a  stress  of  0  N/mnv' .  At  this  leading  the 
■  if  crack  opening  is  0,02  mm  because  of  the  curing  stress.  Under  compressive  stress 
,,  -  -op'.iN/  mm’  delaminat  ions  were  observed  between  the  *1 5  -  degree  and  0  -degree  layers. 


r -'•i-r-'s  s  in  the  cycling  test  with  a  stress  ratio  R  r>[  /o,1  -  -1  is  shown 

!»  Mv  ’  i ,  with  d?  •  -  0  i  -  *100  N/mnr  )  by  the  aid  cf  C-scar.  r  whi'h  were 

a  f r  c  >•  vv‘  numbers  of  -yelps  until  the  test  piece  finally  fract  m-d .  IV* 

■- i'1'  :  n  •  f  tfc  •dr*-  del  am  in  a*,  im-s  in  the  90*  layers  shows  up  clearly  after  .'"'.ono  cycles 
i  1  t  br  •::>•!•••  Me  in  de  1  cm  i  rat  ion  area  at  both  the  test  piece  ed/e  and  at  My-  art  i  fi- 
1 -»r j  r  -i  :•  l  l  as  at  the  natural  delaminat  ions  which  happened  to  l-e  ir«*srnt. 

■  :>'•  i  r.iM  .»  ..  w-T*  i'  .?c  examined  with  a  microscope.  It  war,  seen  that  de  I  ami.  nations 

-  ".?  s  ir,  t  inf  erf  vn  between  the  0-degree  and  the  ^5-degree  layers  also. 

i '  *  !.e  ticx*-  il  lu;»  rat  ions  is  shown  the  planimetered  area  of  the  de  lamination  propagation 
during  fatigue  tests.  Fig.  11  find  Fig.  12.  Doth  specimen  with  embeded  teflon  tabs  or 
I:  f.'atjon  agent  were  tested  with  a  stress  ratio  R  =  -1  and  an  upper  stress  of  o”  =  h 0 0  N/mm*  . 

growth  of  the  different  delaminat  ion  modes  were  planimetered  as  the  number  of  load 
cycles  inc rased. 

>i .  DEL  AMI  NATION  GROWTH  MODET. 

The  monitoring  and  the  accurate  description  of  the  state  of  damage  as  a  function  of  static 
1/  ading  and  time  or  fatigue  loading  is  required  for  successful  development  of  procedures 
for  predicting  the  residual  strength  and  fatigue  life  of  composite  structures.  Most  of 
th*»  available  nondestructive  test  (NOT)  methods,  as  for  example  penetrant,  enhanced  X-ray 
radiography  or  ultrasonic  C-scanning,  require  interruptions  of  testing.  One  have  to 
remove  the  specimens  out  of  the  testing  machine  or  at  least  to  stop  the  test.  The  state 
of  damage  in  composites  can  be  measured  indirect  without  interruption  during  testing  by 
thA  measurement  of  the  change  of  compliance.  (5)  One  of  the  most  promising  techniques 
for  characterizing  delamination  growth  is  based  on  the  rate  of  strain  energy  released, 
il,  wtfh  delaminat  :rr.  growth.  Measured  critical  G  values  have  been  used  in  sophisticated 
analyses  fbj  to  predict  the  onset  of  edge  delaminations  in  unnotched  composite  laminates. 

Jn  our  present  study,  a  technique  was  developed,  employing  the  simple  line-plasticity 
model  (Dugdale-B.arenbl.att  model). 

Fire*,  during  qu as  i -slat  i  o  tension  loadings  and  during  fatigue  tests  the  fransverse 
rie  forma  t  ion  in  My  t  hi  r,ki  !*■*::;:  direction  at  the  edge  of  the  specimen  was  measured  by  a 
special  home-made  displacement  -  transducer.  Next,  stress  distribution  generated  from  a 
fir. ’Me  element  analysis  war  correlated  with  the  observed  damage.  Then,  test  data  and 
analyses  were  used  to  derive  a  closed-form  equation  for  the  characterization  of  the 
de  lamina  t.  i-^n  cnsM.  ,  m.1  ‘he  del  ami  nation  growth  in  relation  to  the  crack  opening  displace- 
rr-  i  f  v  f  tqn  d*>  1  em i nut  i on  ('.vpp). 

*i.l  MEAfUREMKHT  OF  CRACK  OPENING  DISPLACEMENT 

For  the  measurement  of  transverse  deformation  in  thickness  direction  of  the  specimen  we 
ur-’d  *.he  rEVER-MPR-Transducer  (Magnetic  field  depending  resistor).  The  mounting  of  the 
t.  r  .-lr.r.dueer  on  the  t^st.  specimen  is  shown  in  Fig.  3. Three  halve-c;  rc  lo  small  springs  of 
m  |  >*-r?d  j  nr  st  i  ffness  w«r^  glt>erl  in  the  presented  form  at  the  specimen.  Tlie  transducer 
w.r  mounted  by  a  vessel  with  two  incorporated  helical  springs  which  supported 


V  -jir-sitf*  edge  of  the  test  specimen.  The  ter»r  r  '  -  *  *  '  i- 

•.  r.-i  constant  amplitude  tension-tension  fatigue.  Fig.  1'  f  •  *•  -t- 

I!  •  quasi-static  loading  was  interrupted  bv  unloading-  1  «■  "'ll 

i-i-.c k  opening  displacement  due  to  curing  stre.-res  after  ui:  1  .  1  like 

first  test  results  of  a  very  extensive  test  stria. 

ic  tension  test  was  run  at.  a  constant  load  rate  -  f  i  f  ‘  N/mro*/ 

stress  of  a,  =  758  N/mm* .  The  stress-deformation  re  cor  •  •  r  :  .c.l1* 

hand  side  the  nonlinear  lateral  contraction  of  t  h*-  rj  f  •  r  et  *} ■.  kness- 

1  .  r,  it.  The  maximum  contraction  was  measured  as  t  •*  -  '  ■  •  •'  r  ’be 

was  a  creep  deformation  in  ^-direction.  After  rib'  ut  1°  r- ’he 

•  ~  ipn  was  pos  it  i  v.Then  the  specimen  war-  unloaded.  !  ■  '•  the 

r*4l  ax  from  c*,  =  +*,7um  to  +3»5um  .  After  1  minute  the  rj . .  :  loaded 

■  • :  i  •>''><'!’  \  minute  mere,  at  a  *  otai  test  time  op  i<\  r  i*n’.  •  **  ,  y 

:  uv.  starts  as  seen  in  Fig.  t:,a.  The  crack  opening  dir  p  »  “  to 

grew  during  constant  loading  with,  do.' re  as  i  nr  rate  »;  i  ■  •  :’i-m 

.  luring  ‘his  period  the  st  ress-de  J  •urr  i'n  :  r*:i  n  o,  i  unloa- 

-cimen  war  recorded  at  vari'us  i  n*  «-rv  *.  1  :*  a*  re*-u:  in  . 

•  •  h'r*  was  interrupted  for  16  hours.  Thii:  ‘he  ut  ?»?■’  •  ••r-  >  i  ‘  h 

l  o.  -755  N/mm*  .  After  20  minuter;  nv-re  :  hr-  v*"  r  was  in'  •  •  ••  :  •: 

o.  =8*6  fl/mm*  for  3  hours.  Af‘er  ‘his  *  ir>-  i.lin  g:  w  1  •  •  •  -.  <■  about 

•  :  o.i  ing  displacement  Cm  -  1 6 l r. .  TV  :  ■  f  f 1  >•  d<  I  ••  r  .  • 

ui.  Th*»n  the  specimen  was  .^ar^f-iil  u*  •  •  i  in’  two  r.:‘  ‘  "•  ° 

•  >•  fr-m  microscope  (f5EM)  the  distance  i  et  f  *.•■»  cra'T:  f  > ;  “dge. 

r  nlitule  tension-tension  tests  the  urn^-i'Tirvf  i r>rt  :  h 

•  :>:?d  at  various  load  cycle  interval;-.  Fig.  *’  civ  wr  a  ‘  v 

:  rr-o .  The  first  cycle  starts  with  t  h-'  i.v  inn  f  ‘in*  ’  ..  '  nr. 

;  -65O  N/mm?  at  a  stress  rate  of  1  Kfl.miri.  rvn  this  v*  1  •  ■  ‘  ho 

in.:  starts  with  an  upper  stress  of  c?  ‘  "  N  r>?-:  and  the  of 

tress  ratio  R  =  0,1).  The  stress-def  rmnt  i- r.  r*-  \  ".ir*  dur  in-  . 

•  1  at  various  load  cycle  intervals  s°-rs  to  he  ■  uual  to  th.at  >.  .  :  -a 

quasi -static  test.  From  the  beginning  of  eycK  loid.r.r  t!  <•  i  -  -  .  n'4- 

i  :i  of  the  edge  turn  to  positive  values  as  sh  .‘wr:  in  Fig.  1  ■ 

I  acenent  increased  during  t  h"  t-'t-al  1  <"»  1  i  vg  *  u»  'h 

"0  load  cycles  and  increased  f*r  higher  \  \i  .if.- 

•  >..  thee"  was  a  small  mistake  in  th*»  1  -ad  c  n’  r  • !  which.  •  <  ■!  '  u  .*>  :■ 

ITIOO  load  cycles  the  stress  amplitude  was  raised  t  -  :n  :•*  ‘f-  r -■  *: 

>  cycler,  the  measured  crack  opening  displr-.v>rv*n‘  remain  ■*  ••!;•■  t  .  ilv 
.  *  >i*»  freo  edge  de  ]  ami  nat.  ion  was  formed,  lay  it.  bet  wren  /  0  ' 1  •  a  :  load 

;  the  load  history  and  the  d irp  1  acemer.*  history  of  ’  ’  ‘ ! - r*i  n - 

•  states  of  the  damage.  Thr  sinusoidal  cy-  ’it  g  1 t. •  *-  .  1  :  ;  *  -  •  .at 

: :  ■  *  t-y  the  magnitude  of  the  amplitude  i/f  c«  ,  the  moan  va  >  ■  •  •-  :  ihe 

t  h«  stress  cycling  and  the  displacrne.ut  resr-rns** .  At  ‘  iv  •  : 'fate 
-f  tiaonri*  f.  '  ”  S\  "lRo°.  At  t!ie  second  a*  e  II  tim  phas-' 

•  ’  n»  lead  cycles  to  13  =0°  at  2203  cycles .  Retw*^*n  the;-.*'  •  1  •  ■  i'U-s 

•  •  i  »  lar  ions l  the  lateral  displacrmen*  ?r»ci  1 1  a1-  1  ei  and  ‘  h-»  !  1  . .’  i'-'n 

•  i  1  Tf-ruvrit  osc  illation,  both  of  the  ram,  fr'^yjoticy )  supoi’  i*i;  .  'Ives 

•>:  -i?  :  i^C'  .  Therefore  the  amplitude  of  the  measured  <J  i  r-j  l  ■  ■  •••  ?  is 

...  .  T  t  is  the  birth  of  the  first  small  i  nt  orl  uni  na; . w;:i:n 

•  o’-  the  r.prrinert  from  the  tip  of  the  t.rancwerse  oracle  •'  -ply 

’  •  •  .  As  mrjorved  with  th^  stereo  microscope  Id  ose  small  in’  vi  •  i 

anchor  during  increasing  load  cycles.  The  union  of  \  ■  ^-.inar 

•  :  ••  d  "free  edge  dei aminat i on”  ,  effected  the  i rsc rear-  i nr.  of  •  t  .  ;:.ude 

!’  ‘he  displacement  as  seen  L11  Pig.  18  at  the  state  Ill.  TV  M:/ 5 1 ;  on 

•m  **‘'00  and  3010  load  cycles.  The  size  of  the  d*» ] .uminat  ion  •.••  •.  c  c  i  r  ored 
-  «  >n.  The  distance  between  the  crack  tip  of  cWaminnt  ion  and  ;  u-  -Ire 
w  i:--  measured  by  oEM. 


•  )'..ive  predictions  of  the  onset  and  growth  of  del'iminati  u.-  '  *  ‘uite 

■,'K.A  was  used.  Because  delaminations  f«'rm  in  uunot  cli»'d  '  a  ?•. :  u  •  *  •  •  1  u 
c-iinar  stresses  that  develop  at  the  edg«v- ,  a  ■p..,a:;l-Mirr,e- I  ir/ut.  ::ua! 

,  "‘.vs  is  war,  performed.  Some  detailed  of  » lv  .ra’.ysi  s  ar**  t  -•  • ;  1  i  !•.  ( 1 0 )  . 
.  •  •  is  shown  in  Figure  19-  Only  1  / h  of  the  sect  ion  nee d*  c  ^  <  il  vse  j 

,  •  m’ct  r'.v .  This  quarter  includes  8  layers  and  half  Uk>  wid‘h  of  •  v*  t  »st 
-•-■t  at  5  mm  under  the  assumption  that  the  edge  stresses  would  ’uli  off 


illustrated  in  the  following  pictures.  Figure  20  shows  the  de format  ion 
*  h<°  test  piece  section  under  a  tensile  loading  of  a,  -  1  KN  nr'.  On  the 
-n  1  side  picture  edge  ,  the  displacements  of  the  section  in  t!v 
idi«:ated.  All  originally  level  Surfaces  exhibit  distortions,  in  t.h« 

»  dge  under  tensile  stressing;  these  dist  ortions  indicate  ‘  !v*  >  ro':on''e 
os  at  the  edges.  The  distribution  of  tensile  stress  On  in  *  r.e 
•a  f  the  laminate  resulting  from  an  applied  tensile  stres«*  o,  ag»-.»e 
•  which  would  also  be  obtained  from  the  layer  theory.  Fipurf»  21.  Tl-je 
T"  layers,  which  carry  the  main  portion  of  the  load,  are  Uu.J  higher 
i  moan  tensile  load,  and  do  not  change  markedly  towards  the  edge  uf  the 
stresses  <3 2?  arising  as  a  result  of  the  lateral  deformation  in  the 
sre  in  equilibrium  with  each  other  and  fall  nff  to  zero  at  th**  edge,  r»s 
1  rig.  22.  Contrary  to  this,  the  peeling  stresses  On  in  tb*»  layer, 
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r. .  ^ *: -mr*t  from  zero  in  the  middle  of  the  test  piece  in  the  compressive  stress  area 

: ;  i  hen  pass  through  zero  at  about  0.2  mm  from  the  edge  and  display  a  stress  peak 
-I.-!*,  on  the  edge  between  the  90°-and  the  0°-layer,  the  height  of  which  cannot  be 

by  use  of  finite  element  methods.  An  investigation  conducted  in  the 
•  i  ;  •  ..r  on  the  effect  of  the  mesh  size  of  the  elements  t  ( 10  ]  assuming  a  linear  law  of 

•  •  ;  ■  1  -  y ,  indicates  tint  the  peeling  stress  at  this  point  and  at  the  edge  between  the 

j  ;r>:  »h?  0°  layers, albeit  here  as  compressive  stress  o?1  under*  applied  stress 

>  .  •  .  :r,<T  infinitely  high  as  the  edge  is  approached.  The  same  conclusion  is  also  valid 

i  - ;  -».■  rhear  .-.Presses  e,3  ,  which  are  illustrated  in  Fig.  21 J  whereas  the  shear  stresses 
o,.  ?nd  within  and  between  the  individual  layers  on  Fig.  2b  and  Fig. 26  fall  off  to 
zero  at  the  edge.  Of  course  the  manufacture-induced  intrinsic  stresses  between  the 
’■’’dividual  layer's  are  included  in  the  investigation.  With  the  aid  of  the  distribution 
w  of  the  secondary  stresses  033  and  S13  at  the  edge  of  a  multilayer  laminate, 

h.o  r-r-f’k  initiation  and  the  crack  propagation  already  described  which  are  observed  in 

-i t  acting  or  in  fatigue  testing  on  an  undamaged  test  piece  can  be  very  well  explained. 
r,r  ]r\  fracture  mechanics  it  may  be  assumed  that  no  infinitely  high  stresses  are  generated 
v  a  r '■*.*!  1  material  (in  this  case  the  matrix  resin  between  the  fiber  layers)  because 
1  i-iftic  deformation.  1  will  deal  with  the  possiblitiy  of  using  fracture  mechanics 
.•  d  'ui  -i^ lamination  propagation  at  a  later  point. 

a*  ions  w.-rc  also  carried  out  on  test  pieces  with  a  delamination  crack  on  the  edge 
•  ■  r.r-  ‘ - 1  avers .  The  peeling  stress  c3J  and  the_shear  stress  013 

*  h,-»  .** j •  0 k  tip  are  of  the  same  order  compared  with  the  corresponding  edge  effect  stresses. 

•  *  ♦  r**r r  -1  i ;•  r'ibut ions  and  the  deformation  of  the  cross-section  showed  reasonable 

-  ’  r  ;  .  wjth  »  he  observed  damage  that  developed.  Indeed,  examining  O33  and  interlaminar 

■  -r  distributions  are  helpful!  in  identifying  likely  del  ami  nation  sites.  However 

•  . :  ■  :  *  i  n  .-'f  interlaminar  stress  distributions  by  finite  element  analysis  are  very 

•'  ;■  v.  m  -r  only  useful  for  special  cases  of  modelling  damage  growth  qualitatively 

r  re  ,-n.i  t  ude  of  calculated  peak  stresses  at  the  edge  varied  with  mesh  size. 

•  •  •  r  •  ,  an  alternative  approach  based  upon  the  yield  stress  criterion  was  used  to 

:  :!  ‘  !!;»'  d«l  ami  nation  size. 

ATI  f-M  OF  THE  DELAM I  NATION  OPENING  DISPLACEMENT 

■  r  t  h"  results  of  finite  element  calculations  there  exists  a  stress  concentration 

‘■■dge  the  considered  multilayered  laminate  in  the  interface  between  the  90-degree 

•  e_i.^  -r-.ri  pi  ir><5 .  Tne  behavior  of  the  matrix  material  in  the  presence  of  this  stress 

•  :  r  ■,«  :  >r;  interface  can  be  estimated  under  the  assumption,  that  the  extend  of  the 

•  !;  rt  rerr- ..train  behavior  islimited  on  a  small  plastic  zone  around  the  stress 

•.  1 . •  last  ic  materials  cracks  can  form  and  propagate  with  very  low  velocities  and 

••]•  --rate  slowly  [ll]  .  Therefore,  estimation  of  the  crack  tip  velocity  and  the  crack 
■  -»v.  any  r\mo  become  necessary  to  estimate  the  lifetime  of  a  structure. 

*r  £>laatic  stress  analysis  shows  that  stresses  at  the  90-0-degree-interface  and  at 
.  -r  --k  t  ip  become  unbounded.  [10)  .  However,  materials  exhibit  a  yield  stress  Oj  above 
v.k  :  *h  they  deform  plastically,  and  thus  there  must  be  a  plastic  zone  around  the  stress 

rv”en‘.  rat  ion  which  limits  the  size  of  any  stresses.  One  can  model  this  plastic  behavior 
in  a  .*.  imple  line-plasticity  model  (Dugdale-Barenblatt  model). 

T  •  m  > lyse  the  crack  opening  displacement  c  of  the  delamination  (CODD)  as  a  function 
.  ;•  ;t  rip  delaminat  ion  size  a  and  of  the  axial  stress  o*  of  the  laminate 


(  1 )  c  ^  f In. 6,» 


simple  model  was  used.  Considering  a  free  edge  strip  delamination  of  the  size  a  the 
m'-Trured  value  of  the  crack  opening  displacement  cm  consists  of  (see  Fig.  27) 


(  P  )  r  „  s  c  ,  »  r.  j  •  c  j  ♦  r.  t  =  c  *  c  3  . 


>.  ,  i  r-  Mte  displacement  due  t.o  different  transverse  contractions  v12  of  the  single  layers 
i *  u h  unbalanced  ( unsymmet ric )  delaminated  strip  induced  curvature,  r  means  the 
;!  i  "p  1  .» cement.  due  to  curing  stresses  induced  curvature;c,  the  displacement  due  to 
transverse  contraction  v,j  and  the  displacement  due  to  yield  stress  at  the  vicinity  of 

u:i-'  ’•rack  tip. 

T  ■ii'ulato  » h«  displacement  c^,  parts  of  the  general  constitutive  equations  for  the 
(see  for  example  (12)  ) 


:-"’1  , 
irf 
f  :"ne 
Th  a  s 


where  N  and  M  are  the  resultant  forces  and  moments  acting 
’ailed  the  extensional  st if fness  matrix ,  coupling  stiffness 
ss  matrix,  respectively,  e?  is  called  the  midplane  strain 
curvature  can  be  expressed  for  the  y-direction  as 


on  the  laminate, 
matrix  and 
and  k  the  plate 


d '  a 


.•a in  coupling  terms  of  the  general  constitutive  equations  (3)  for  thi3  case  are 

M  y  s  B| j  E J  ♦  D2j  ky 

‘■h*4  momer*.  My  =  ?.  Combining  Eq8 .  (*0  and  (5)  and  integration  leads  to 

SI  312  ( a » Oo ) 

,,,  n"'e’  io.q0]! 


the  subscripts  denote  the  coordinate  axis  and  the  superscripts  (1)  and  (2)  denote 
^'laminated  parts  of  the  specimen.  Therefore  the  displacement  due  to  different 
v -me  contraction  is; 


2  l  o'"  oi1;  I 


“m^nf  c  j  due  to  curing  stresses  induced  curvature  can  be  expressed 


•‘  ?  -j  '<i  ♦<?  Ha*a0J 

§’  arid  9L'are  the  individual  curvatures  of  the  delaminated  parts  of  the  specimen, 
ransverse  contraction  of  the  whole  laminate  is 

■  3  *  tl  v,3  =  tl  v,3  d 

d  -  d^  *  d  ^  is  the  thickness  of  the  specimen.  At  least  the  displacement 
constant  yield  stress  can  be  estimated  like  a  clamped  beam  with  partly  uniform 
isM'ihution  of  the  length  a0  as 


,  1.) 

o’  si  / 

t  ♦ 

1  t-'Vi 

»»  I 

’  a  \ 

to,  / 

»tir>g  Fqs.  8  to  11  into  Eq.  2  yields 


1  i  B'.’l  e‘.f  \  la-a,r 

!  -  Vl)dl  .  _L 

IS1’1. 

If. 

/  '  .  '\i 

fu±  a  ’ 

L 

\  °o  ~r 

2 

J  ? 

1  y  « 

3  1  ID  *0,1 

1 1>«  o',1,1 ) 1 

l  2  o0  , 

I  1 ' 6' 

r: 

•  0,/a,, 

1  " 

all  projerties 

are  known  exept 

the 

length  of 

the  plastic 

7.one  o„  . 

Ar  sum  i  ng 

••  i  •;  *  material  property,  we  can  estimate  the  value  of  o0  by  the  mear.ur  cm-  nt  -f  f(T, 
inr-f  ~.pf>r  i m»n  at  delamination  onset  stress  o,  .  In  this  case  th'*  ''rack  length 
and  therefore  from  Eq .  12 

•  fir*  0,1. 

'  tested  laminates  of  the  stacking  sequence  [o  / -^5/0/90  Js  fabricated  from 

r'^'P  graphite  epoxy  prepreg  tape  th  amount  of  calculated  plastic  zone  is  o,  -  0,5mm. 
f«cMve  modulus  properties  of  each  unidirect ional  ply  used  for  calculation  are 

EL  =  150  kN/mm? ,  ET =  10.000  kN/mm’ 

GLT  5  5  kN/mm*  and  vlT  =  0,30  . 

Properties  the  extensional  stiffness  A„  ,  the  coupling  stiffnesses  bV?  and  b\\: 
■,incr  stiffnesses  Djj  and  D1,2,  and  the  curvatures  due  to  curing  stresses  s'”  and  9I:)  * 

ilc'ilnf*d  by  using  classical  lamination  theory  for  both  delaminated  parts. 
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~>ur  example  these  properties 

A i  •  -  l  7  v  1 2 6 I  KM/ mm 
-c  -  *,h7t>5  KM 
r-  ,5  <«0  KM,  mm 
b"'  -  0  ,  >i7  5!J  <h  1/mm 
3  j  r  ' )  ♦ 1  - r>  r.N/nurr 


-0,1 7 39  KM 

0,9136  KN/mm 
0,000956  1/mm 
0,758  KM/mm? 


curvatures  q  were  estimated  for  a  difference  of  temperature  AT  =  100  K  between  test  tempera- 
*  and  glass  transition  point.  The  strip  deamination  size  a  was  calculated  by  the  Eo.12  as  a 
fur.  r-.ion  nf  the  crack  opening  displacement  c  m  .  The  results  were  compared  with  the  size 
■>f  *■»  !ge  delamination  for  different  delamination  states  on  several  test  pieces.  Therefore 
;>  --'i nation  sir.es  were  observed  by  C-scans  and  by  a  scanning  elektron  microscope  (SEM) 
if'iM’  ending  the  test  and  cutting  the  specimen  transverse  to  the  load  direction.  A 

ari  r on  made  between  the  measured  and  the  calculated  stress-displacement  curves  at 
i : :  1  ••  r^nt  delamination  states  i s  shewn  in  Fig.?8  for  a  quasi-stat  ic-t.est  and  in  Fig.  29  for 
r  I  t.ppt-.  The  estimation  was  carried  out  under  the  assumption,  that  ai  (  in  the 

;  0  >>f  ci p  a*  Fq .  12}  is  growing  linear  with  the  increasing  stress  q,  up  to  the  critical 

*  v-.'S  o.r  of  ‘he  uelamination  onset.  (  a0  =  q0ox/g,c  )  -  A  good  correlation  was  achieved 

'  *  be  measure!  and  calculated  stress-displacement  curves. 

«'•  ire  r till  working  on  the  formulation  of  delamination  growth  as  a  funktion  of  time  t, 

:  c> ,  .  t  }  and  a  =  f  (o0.R.N.f  )  ,  and  on  the  problem  of  strain  energie  release  rate 

ia‘o.1  with  the  delamination  growth  also  in  relation  to  the  change  of  specimen  axial 
lance  .  (S]  . 

:'i"*  objective  is  to  combine  the  degradation  due  to  growing  delamination  size  in  tension- 
ms  version  cycling  with  the  residual  compressive  strength  of  fatigued  and  partly  buckled 

;  mens . 

FATIGUE  AND  RESIDUAL  STRENGTH  MODEL 

•  :rir.g  quasi-static  compressive  testing  or  fatigue  testing  with  considerable  amplitudes 
f  '••'mrroccive  stress,  failure  occurs  by  the  buckling  or  kinking  of  individual  fibers, 

•”  f undies  or  fiber  layers  which  are  more  or  less  supported  by  the  matrix  resin.  The 
"■'V*  i. vr-  effect  is  reduced  as  the  delamination  progresses,  under  both  tensile  and 
1  f;  —  >  v**  loading.  Fracture  occurs  finally  when  a  portion  of  the  cross-section  buckles 
:  > b*  pr-m-ininr  cross-section  is  no  longer  sufficient  for  the  transmission  of  the  applied 
The  effect  of  the  ant  ibuckl  ing  support  which  holds  the  test  piece  either  over  the 

■  ‘  ; : area,  or  -d  a  center  line,  or  on  the  2  edges  of  the  gauge  area,  must  be  taken  into 

•  • :  >  rn-  ior; .  Whilst  the  delamination  state  of  a  test  piece  in  multi-layered  CPRP  has  no 

influence  on  the  residual  tensile  strength  or  on  the  tensile  fatigue  strength, 

■  ?  pieces  with  de 1 ami nat tons  are  especially  sensitive  to  compressive  loading. 

•  *  he  denomination  state  of  a  multi-layered  test  piece  during  cyclic  loading  is  to  be 
•■'■*ribe  ,  a  formulation  from  fracture  mechanics  may  be  used  for  the  purpose  [  1 3  »  1  ^  )  .The 
r^wth  of  crack  length  (  or  in  this  case  the  delamination  size  a)wit’.i  increasing  number 

f  cycles  M  is  shown  by 


wh^re  0  is  a  reference  stress  in  the  vicinity  of  the  delamination  crack  tip;  g,  a?  and 
m  are  constants  to  be  determined  empirically.  By  integration  and  insertion  of  limits 
for  M  =  0  with  a  =  a,  and  for  N  -  Na  with  a  =  a„  the  number  of  loading  cycles  N  is 
obtained  which  corresponds  to  the  delamination  size  a«j  : 

'  *■ r  }  Nd  s  (o'rt  n’1  -  /  |1-ml  g0* 

:>'W  we  assume  that  the  test  piece  fractures  due  to  compressive  loading  when  the  delamination 
"  y-'-r  a  reaches  a  critical  value  { 1 5  ]  .  As  shown  in  Appendix  A  the  critical  delamination 
rir.e  a  for  a  cyclic  loaded  specimen  can  be  expressed  as 


re  0e  is  the  critical  compressive  stress,  of  the  compressive  strength  of  the  virginal 
r\ -vr (men ,  A  the  cross  section,  h  the  thickness  of  delamination  strips,  E  the  modulus  of 
the  non  delaminated  and  E,  the  modulus  of  the  delaminated  cross  section.  In  the  case  of 
fit ; rue  failure  the  specimen  will  fracture  after  N  load  cycles  at  the  compressive  cycling 
~treps  amplitude  oF  -  ol  when  the  critical  delamination  size  art  -  nF  is  reached. 


care  of  residual  strength  test  of  a  fatigued  specimen  .after  Mq  load  cycles  the 
il  damage  size  -  ap  can  be  expressed  by 


*  •  7 


A_  _E_ 
h  E, 


where  is  the  residual  strength.  The  damage  state  D  can  be  formulated  as  a  fun  ti  or.  of 
dimensionless  number  of  load  cycles. 


Substituting  N*  and  NF  by  Eqs.  15,  17  and  l8  yields 


For  our  further  considerations  we  will  use  the  simple  expression  for  the  damage  state  due 
to  de lamination 


(21)  0 


( 0* -gw  y 

\0t-OF  / 


The  validity  of  this  formula  was  controlled  by  fatigue  tests. 
f .  FATIGUE  TESTS  AND  RESIDUAL  STRENGTH. 

Fif igue  tests  were  carried  out  on  type  I  specimens  with  a  gauge  length  of  15  nr  f  - 

rtM  ermine  the  number  or  load  cycles  to  fracture  under  constant  stress  amplitude  w:*h  the 

::r  rpss  ratios  of  R  =  -1.0  and  R  =  0.1.  The  results  are  illustrated  in  Fig.  V' . 

The  specimens  were  supported  over  their  entire  area  against  structural  bucklinr.  In  '-rder 
r:  determine  the  damage  state  of  the  test  pieces,  additional  specimens  were  r:;1  .;'*'***  i  *  o 
fatigue  loads  at  R  =  -.1.0.  After  certain  number  of  cycles  N  the  residual  s*.t  • 

checked.  The  results  of  the  residual  strength  tests  are  plotted  versus  the  number  T 

'-•■‘•iir.g  cycles  N*  in  Fig.  31.  The  residual  compressive  strength  falls  off  :  ra  u  •  1  v 

v.- i Mi  the  approach  to  the  number  of  cycles  to  fracture  Nf  ,  down  to  the  fatigue  r*  r*v  «;  <*P . 

'I  ho  residual  tensile  strength  shows  no  significant  change  in  comparison  with  the  static 
strength. 

Trv  illustration  of  residual  strength  investigation  in  this  form  is  very  conn!**.  !n  Fir.1;? 
the  damage  state  D  of  the  test  specimens  is  shown  with  reference  to  the  pr'»rr‘»?:*. 

delamination.  Since  the  test  results  are  naturally  scattered,  the  equation  ,  P 
expanded  by  a  number  of  terms  which  take  account  of  this  scattering. 


)  D 


Ot  *  (  1  -0C) 


(fl 


Q>  -qp  \  1-m 

Of -Of  / 


Wh:rr*»  of  -  of,  ♦  oC 

•>nl  eft  1  statistical  scatter  for  a  given  failure  probability  P  based  on  static  r  ♦.  rei.pt  h. 

<*:  -  average  damage  of  a  virginal  test  piece  with  a  failure  probability  of  V  *  , 
ft  -  statistical  scatter  for  a  given  survival  probability  Pq  based  on  the  toil  cycles 
to  failure. 

7h-*  corresponding  scatter  distributions  and  curves  of  equal  survival  probability  P0  -  90*, 
r  7-  and  10%  are  shewn  in  Fig.  32.  The  damage  D  to  be  determined  in  the  residual  strength 
investigations  is  calculated  from  equation  (20)  and  entered  into  Fig.  32.  The  exponents 
m  determined  from  the  test  results  as  m  s  0,^5.  As  can  be  deduced  from  Fig.  3?,  the 
•  U  1  am j nat ion  progress  can  be  described  satisfactorily  with  the  aid  of  equations  (17), 

'  !■”*;  an  (21),  taking  account  of  the  scattering  of  the  test  results. 

7.  NUMMARY 


in  r  he  present  study  a  technique  was  developed  to  characterize  the  onset  and  t.  ho  growth  of 
laminat ions  in  multidirectional  composite  laminates.  The  delamination  that-  farmed 
un  static  and  cyclic  loading  in  the  interface  between  the  layers  of  the  ]  ami  nat-- • 
v * ■ ! .  h-t^rminod  by  ultrasonic  C-scans.  During  test  loading  the  transverse  contract  i  n  and 

ty,-  ‘ rnnrverse  crack  opening  displacement  due  to  peeling  stresses  were  monitor*-  '.  The 
or  •‘■rvM  onset  of  de  laminat  ions  was  correlated  with  the  stress  distribution  generated 
fr'-m  a  finite  element  analysis.  The  test  data  and  the  analysis  were  used  to  derive  a 
r l ed- form  equation  for  the  delamination  size  assoziated  with  the  onset  of  cracks  and 
t h“  'rick  opening  displacement.  Those  defects  progated  slowly  owing  to  interlaminar  stresses 
up  to  an  area,  which  is  critical  against  buckling  of  parts  of  the  delaminated  test 

in  the  case  of  tension-compression  fatigue.  A  correlation  between  de laminat ion 
sir*  and  residual  compression  strength  of  fatigued  specimens  was  derived. 
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APPENDIX  A 

A  delamination  propagation  model  under  peeling  stresses  Q33  (mode  I  condition)  may  be 
expressed  in  a  form  similar  to  a  crack  propagation  model  in  metals  and  is  given  by  the 
Eq.  14*  Let  the  initial  defect  characterized  by  a  delamination  size  a, =  0  and  the  size 
a^  after  a  number  of  loading  cycles  Nd  .  The  integration  of  Eq .  yields 


(Al)  =  a’d^/ll-m)  g  o* 


Assuming  the  edge  delaminations  ai  in  the  multilayer  composite  develope  symmetrical 
from  the  edges  to  the  center  planeof  the  specimen  as  shown  in  Fig^  Al.  Than  we  may  define 
the  delaminated  area  A*  =  a, h i  and  the  non-delaminated  area  A  -  A  -  A*  ,  where  A  -*b  d 
means  the  total  cross  section.  The  load  distribution  on  both  parts  of  the  cross-section 
is 

(A2)  Pa  *  !>»•,  *  Pi 

i:1 

It  was  observed  in  tension-compression  fatigue  tests  and  residual  strength  tests  that 
parts  of  the  cross  section  failed  by  buckling.  The  buckling  load  of  the  delaminated  section 
i  which  buckles  at  first  can  be  expressed  by  [15]  ,  (l6] 


(A3)  Pj 


k,  A* 


TC2Ei 
12(1-  vLV 


where  k;  is  the  compressive  buckling  coefficient  for  flat  plates  with  one  single  supported 
edge  and  three  clamped  edges  as  shown  in  Fig.  A  2  [l6j  .  The  value  of  the  buckling  coeffi¬ 
cient  is  2,0  >  k,  >1,2!!  During  fatigue  loading  the  delamination  will  grow  up  to  the 
delamination  size  ae  .  Than  the  strip  becomes  instable .  It  will  be  assumed  that  the 
postbuckling  load  of  the  delaminated  strip  remain  constant  and  equal  the  buckling  load 
during  increasing  strain  (see  Fig.  A3  and  [17]  ).  With  increasing  number  of  load  cycles 
the  delamination  growth  will  continue  and  the  buckling  load  will  decraese  after  Eq .  (A3). 

See  Fig.  AU,  When  a  certain  part  of  all  delaminated  strips  become  instable  the  remaining 
non-delaminated  cross  section  will  faile  because  of  exceeding  its  static  strength  (or  the 
failure  strain) .  Those  considerations  are  transferable  on  the  behavior  of  a  delaminated 
laminate  on  residual  strength  test. 

The  residual  strength  at  a  certain  number  of  load  cycle  N  can  be  expressed  from  Eq.  (A2) 


(A'O  Oh  --  £  O.;  -  aS  (l-  -£2^1 ) 

where  fl**  =  P,  /  a,  h;  and  Pj  are  the  buckling  loads  after  Eq.  (A3).  To  form  a  very  simple 
relationship  between  residual  strength  and  delamination  size  we  assume  that  h,  =  h  =  const, 
E",  o,  h,  =  oN  h  (where  oN  the  delamination  state  after  N  number  of  load  cycles). 

<3 *«  =  E,  e  and  Oj  =  E  e  .  With  regard  to  buckling  stress  and  the  above  simplifications 

the  residual  strength  yields 


r  E,  1  I  2hj  , 
IT  T  \  TTSi'  "  h) 


where  k,  =  2 
be  expressed 


was  used.  Assuming  I<3n /o»  - 1  )7  A7 /4h*>2h?  / the  critical  delamination  size 
for  the  residual  strength  and  for  the  failure  load  cycle  as 

—  and  oF  *  |l-  )  —  A.  respectively .  (see  Eqs .  ( 16 )  to  ( 18 )  ), 

h  E,  »  Of  I  h  E, 


can 


T^siic  stress  0.  =700N/mm‘ 


C,  =  800N/mm2 


O,  =  9C0  N I  mm 


Fig.  6  Growth  of  the  edge  delaminations  and  of  the  artificial 

inflation  agent  delamination  under  increasing  tensile  loading 


T1 


T 
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Fig.  7  The  influence  of  moisture  profile  and  moisture  content  on  the  formation 
of  free  edge  delamination  during  tensile  tests 
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Fig.  8  Growth  of  artificial  delamination  under  step-by-step  compressive 
loading.  Ultrasonic  C-scans  of  a  test  specimen  containing  teflon 
circles  implants 
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Fig.  9  Crack  types  resulting  from  edge  effects  after  3.10  load  cycles. 
Stress  amplitude  ©a=  +  500  N/mm2 ,  Stress  ratio  R  =-l. 
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Fig.  10  Delamination  states  on  increasing  numbers  of  tension-compression 
loading  cycles.  Ultrasonic  C-scans  of  a  test  specimen  containing 
teflon  circles  implants 


Fig.  18  Load  cycling  and  displacement  response  for  three  various  states 
during  fatigue  test 

State  I  :  Creep  phase.  State  II:  Onset  of  delamination 
State  III  :  Growth  of  delamination 


ng.  19  Coordinate  System 


re  formations  of  the  cross-section  result inp  froir.  loadir.p  or 
3x  =  1000  N/mm* 


Pig.  22  Transverse  stresses  o2?  in  the  individual  layers  resulting  from 
longitudinal  stress  O,  =  1000  N/mm* 


Pig.  23  Transverse  stresses  033  in  the  individual  layers  resulting  from 
longitudinal  stress  o,  =  1000  N/mm* 


Fig.  26  Shear  stresses  <323  in  the  individual  layers  resulting 
longitudinal  stress  o*  2  1000  N/mm* 
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Fig.  27  Composition  of  the  displacement 
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Fig.  30  Fatigue  strength  diagram  for  graphite/epoxy  laminates. 
Stress  ratio  R  =  -1  and  R  =  0.1 
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Fig.  31  3-N-curve  and  residual  strength  after  cyclic  loading. 


Distribution  of  Load 

Cycles  to  Failure 


APPENDIX  K 

Zur  Auswahl  eines  CFK-Mehrschi chten- 
laminats  fur  Versuche  im  Zugbereich 

von 

Or.-Ing.K.  Rohwer 

Institut  fUr  StrukturmechanFk  .Braunschwei g 


1 


DEUTSCHE  F0RSCHUN6S-  UNO  VERSUCHSANSTALT  FOR  LUFT-  UNO 
RAUMFAHRT  E.V.,  INSTITUT  FOR  STRUKTURMECHANI K 


Interner  Bericht 
IB  131-82/01 


Zur  Auswahl  eines  CFK-Mehr sch i ch ten  1  am i na ts  flir  Versuche 
im  Zugbereich. 


Braunschweig,  Januar  1982 

Der  Bericht  umfaBt: 
19  Seiten  mi t 
3  Tabellen  und 
8  B  i  1  d  . rr 


Der  Bearbeiter: 

{ Dr . -  I ng . K .  Rohwer ) 


2 


G1  i  ederung  Sei te 

1  .  Aufgabenstel 1 ung  4 

2.  Oberlegungen  zur  Auswahl  4 

2.1  E i nschrankende  Voraussetzungen  4 

2.2  Anzahl  der  moglichen  Komplexionen  5 

2.3  Bewertungskri teri um  6 

3.  Fal 1 untersuchung  7 

3.1  Berechnete  Laminate  7 

3.2  Bewertung  7 

4.  Entscheidung  8 

5.  Schrifttum  8 


3 


A 


1  .  Aufgabenstel 1 ung 

Strukturen  aus  kohl ef aserverstarktem  Kunststoff  (CFK)  werden  in 
der  ube rw i egenden  Mehrzahl  aus  mehreren  Schichten  m  i  t  unter- 
schiedlicher  Faserrichtung  lamelliert,  um  das  Tragvermbgen  des 
Werkstoffes  moglichst  gut  zu  nutzen.  Ein  Hindernis  fur  eine  um- 
fassende  Anwendung  dieses  vi e 1 versprec henden  Materials  auch  in 
Primarstrukturen  ist  die  immer  noch  unzureichende  Kenntnis  des 
Versagensvorgangs  unter  mechani scher Belastung  und  bei  der  Wir- 
kung  un  tersch  i  edl  i  cher  Temperatur-  und  Feuc  hteei  nf  1  iisse .  Fur 
die  Erweiterung  des  Wi ssensstandes  sind  unter  anderem  Laborver- 
suche  an  Probestaben  aus  Mehrsch i ch ten  1  ami na ten  erforderl ich. 

Bei  Versuchen  an  solchen  Probestaben  unter  axialer  Zugbelastung 
haben  Pagano  und  Pipes  (11  ,( 2  1  ein  Oel ami n i eren  entlang  des 
freien  Randes  beobachtet.  Sie  fiihrten  dieses  Phan  omen  auf  das 
Auftreten  i n ter  1  am i na rer  Scha  1  s pannungen  c^gzurtick.  Andere  Unter- 
suchungen  von  Puppo  und  Evensen  l 3 1  ,  Pipes  und  Pagano  (4  1  und 
Rybicki  t 5l  haben  hohe  i n ter 1  ami na re  Schubspannungen  o^an  den 
freien  Randern  von  Proben  mit  Winkellagen  ergeben.  Art  und  GroBe 
der  Randspannungen  sind  von  dem  gewahlten  Lagenaufbau  abhangig. 

Es  ist  daher  moglich,  durch  eine  giinstige  Wahl  der  Schi  chtenf  ol  ge 
die  Randspannungen  zu  reduzieren  und  damit  das  Auftreten  von 
Delaminationen  einzuschranken.  Delaminationen  kbnnen  die  Ober- 
tragbarkeit  der  MeBergebnisse  auf  das  Verhalten  groBflachiger 
Strukturen,  bei  denen  freie  Rander  nur  eine  untergeordnete  Rolle 
spielen,  erschweren.  Fur  eine  Reihe  von  anstehenden  Zugversuchen 
an  Mehrsc h i c h ten verbunden  aus  CFK  soil  deshalb  im  folgenden  ein 
Lagenaufbau  mit  moglichst  geringer  Delaminationsgefahr  ausge- 
wahl t  werden . 

2  .  Oj) e_r  1  e gungen  zur  Auswahl 

2. 1  EjnschrankendeVoraussetzungen 

Die  he rzus te 1 1 enden  Proben  sollen  380  mm  lang  und  mindestens 
10  mm  breit  sein  und  aus  16  Schichten  zu  je  0,125  mm  Dicke  be- 
stehen.  Eine  Einschrankung  fur  den  Lagenaufbau  besteht  in  der 
Forderung  nach  Symmetric  und  Aus geg 1 i chenhe i t ,  d.h.:  jede  der 
beiden  symmetri schen  Lam i na tha 1 f ten  soil  die  gleiche  Anzabl  von 
Winkellagen  mit  Faserrichtungen  von  +6  und  von  -0  enthalten. 
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Die  Einhaltung  dieser  Forderung  bewirkt  die  Entkopplung  der  Ver- 
formungen.  Unter  axialer  Zugbelastung  wird  daher  die  Achse  des 
Zugstabes  gerade  bleiben. 

Die  Belastung  der  Probe  soil  an  beiden  Enden  durch  Klemmbacken 
iibertragen  werden.  Damit  die  f U r  CFK  recht  hohen  Krafte  auch  ein- 
geleitet  werden  konnen,  mlissen  die  Fasern  der  beiden  auBen  lie¬ 
ge  n  d  e  n  Schichten  in  Langsrichtung  {  ©  =  0°)  orientiert  sein. 

Unter  den  genannten  Voraussetzungen  kann  man  noch  liber  die  Faser- 
richtungen  von  sieben  Schichten  entscheiden;  dabei  muB  allerdings 
die  Forderung  nach  Ausgegl  ichenhei  t  beachtet  werden. 

2.2  Anzahl_dermoglichen_Komglexionen 

Die  durch  ihre  Reihenfolge  un tersch i ed 1 i c hen  Kompl ex  i  onen  von 
n  Schichten  heiBt  Permutation  Pn  .  Die  Anzahl  aller  mbgllchen 
Permutationen  ergibt  sich  zu 


Fiir  n  =  7  erhal  t  man 

P?  =  7  !  =  5040  . 

Nun  wird  in  der  praktischen  Anwendung  aber  kaum  ein  Laminat  mit 
sieben  verschiedenen  Faserrichtungen  hergestellt  werden.  Einige 
Lagen  werden  gleiche  Orientierung  aufweisen.  Tritt  unter  den  n 
Elementen  a  me  1  das  Element  a,  P  mal  das  Element  b  u.s.w.  auf, 
dann  ergibt  sich  die  Gesaamtzahl  der  Permutationen  zu 

n  ! 

P  =  - 

n  a!  P !  • 


Um  die  Zahl  der  moglichen  Permutationen  zu  verringern,  wird  in 
einer  ersten  Wahl  vorgegeben,  daB  das  Laminat  je  eine  +45°-, 
eine  -45°-  und  eire  90°-Schicht  sowie  4  (+1)  0°-Schi chten  auf¬ 
weisen  soli.  Damit  ergibt  sich 

P  =  7  !  /  4  !  =  5-6-7  =  2'0  . 
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In  einer  zwe i ten  Wahl  wird  die  90°-Schicht  durch  eine  zusatzliche 
O^-Schicht  ersetzt,  so  daB 


Auch  diese  Zahl  ist  fiir  eine  sy s tema t i sc  he  Analyse  noch  zu  groB. 

Die  groBe  Anzahl  moglicher  Komplexionen  verdeutlicht  die  An- 
passungsfahi gkei t  von  geschichteten  Faserverbundstrukturen  ,  stellt 
aber  auch  eine  Schwierigkeit  bei  der  Wahl  eines  giinstigen  Schicht- 
aufbaues  dar.  Man  kann  nur  intuitiv  wenige  auswahlen  und  diese 
miteinander  vergleichen. 

2.3  Bewertungskri teri um 

Angestrebt  wird  ein  Laminat,  fLir  das  die  Gefanr  einer  Randdela- 
mination  unter  Zugbeanspruchung  moglichst  gering  sein  soli.  Am 
freien  Rand  treten  Spannungen  o,,  ,  o33  und  o13  auf.  Es  wird 
angenommen,  daB  nur  die  beiden  1 etztgenannten  Komponenten  die 
Delamination  beeinflussen. 

Finite-Element-Rechnungen  haben  ergeben,  daB  die  groBten  inter- 
laminaren  Schub-  und  Schal spannungen  an  den  Randpunkten  auftre- 
ten,  an  denen  die  Faserrichtung  wechselt.  Andererseits  weisen 
Konvergenztests  I  6  1  darauf  hin,  daB  die  Spannungen  hier  unend- 
lich  groB  werden.  Das  resultiert  aus  der  Fiktion  einer  unend- 
lichen  Elastizitat  der  Werkstoffe.  Re ales  Material  baut  in  der 
Regel  die  Spannungsspi tzen  durch  n i chtl i neares  Material verhal ten 
ab.  Allerdings  muB  damit  gerechnet  werden,  daB  an  den  Stellen, 
an  denen  man  Spannungsspi tzen  berechnet,  d.h.  also  wo  die  faser¬ 
richtung  wechselt,  eine  Delamination  initiiert  werden  kann. 

Die  berechneten  Spannungen  an  den  Randpunkten  sind  nicht  ver- 
wendbar.  Stattdessen  wird  vorgeschlagen,  entlang  der  Schicht- 
grenze  nach  innen  zu  gehen  und  die  dem  Rand  am  nachsten  ge- 
legenen  verlaBlichen  Spannungswerte  fur  die  Beurteilung  der  RiB- 
gefahrdung  heranzuz i ehen .  Um  einen  Vergleich  anstellen  zu  konnen, 
mu'ssen  die  so  bestimmten  interlaminaren  Schubspannungen  0^3 
und  Scha 1 spannungen  o33  zu  einem  Einzelwert  verkniipft  werden. 
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In  Anlehnung  an  v.  Mieses'  Kriterium  wird  als  Vergl eichswert 


Ov  =  sign  (  O33I  V°332  +  3(o13) 


gewahlt.  Je  groBer  ov  ist,  umso  groBer  wird  die  Delaminations- 
gefahr  e  i  ngescha tzt . 

3 .  Fall untersuchung 
3.1  BerechneteLami nate 

Aus  der  Vielzahl  moglicher  Schichtaufbauten  sind  die  folgenden 


v  i  er 

ausgewahl t 

worden  : 

1  ) 

1  0  2  ’ 

+  45,  0 

2 ,  -45  ,  0  , 

90  ls 

2) 

io5, 

+  45,  - 

45  ,  90  1  s 

3) 

I  02> 

90,  02 

,  +45,  0, 

"45  1  s 

4) 

1 02, 

+  45,  0 

2  ,  -45 ,  02 

>s 

Zur  Berechnung  der  Spannungen  unter  Axiallast  wurde  die  in  l  6 1 
beschri ebene  E 1 emen tte i 1 ung  benutzt.  Ergebnisse  fur  die  ersten 
drei  Laminate  sind  in  den  Bildern  1  bis  6  dargestellt.  Man  er- 
kennt,  daB  das  laminat  2  wegen  gleichzeitig  hoher  Schub-  und 
Scha I spannungen  am  Obergang  zwischen  der  90°-  und  der  -45°- 
Schicht  ungeeignet  ist.  Es  wird  daber  von  vornherein  ausge- 
schal tet. 

Die  Bilder  7  und  8  zeigen  jeweils  einen  Vergleich  der  Rand- 
spannungen  zwischen  den  Laminaten  3  und  1,  sowie  4  und  1.  Im 
Laminat  4  fallt  der  sehr  geringe  Scha 1 s pa nnungsante i 1  auf. 

3.2  Bewertung_ 

Line  Analyse  der  berechneten  Spannungswerte  entlang  der  Schicht- 
grenzen  zeigt,  daB  bereits  der  nachstinnere  Wert,  eine  Element- 

_  "3 

breite  =  (1  /  2  561-Lami  na  td  i  c  ke  »  7.8  *10  mm  void  Rand  entfernt, 

als  verlaBlich  angesehen  werden  kann.  Die  Genauigkeit  kann  da- 
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durch  kontrolliert  werden,  daB  die  Spannungen  nicht  iiber  die 
Schi chtgrenzen  hinweg  gemittelt  wurden.  Sie  muB  aus  Gleichge- 
wichtsgrlinden  ubereinstimmen.  Fur  das  benutzt  WeggroBenverfahren 
cha rak ter i s i eren  daher  die  relativen  Spa nnun gsd i f f erenz en  den 
Fehler  in  den  Ergebnissen. 

Die  Tabellen  1,  2  und  3  enthalten  die  berechneten  Spannungswerte 
an  den  Schi chtgrenzen .  Als  Fehler  wurde  jeweils  die  absolute 
Spannungsdifferenz  auf  den  absoluten  Mittelwer*  bezogen.  Die 
fast  liberall  gute  Genauigkeit  schon  in  sehr  geringem  Abstand 
vom  Rand  wird  auf  die  Verwendung  von  Elementen  mit  quadrati- 
schen  Ansatz  zurlickgefuhrt.  Die  Tabellen  enthalten  auBerdem  noch 
den  jeweiligen  Vergleichswert  ov  . 

4 .  Entscheidung 

Das  Laminat  1  weist  in  der  90°-Schicht  einen  Wert  von  ov  =  0, 12353 
auf.  Vorausgegangene  Versuche  haben  auch  gezeigt,  daB  dort  die 
Delamination  einfallt.  Das  Laminat  4  hat  zwischen  zwei  0°-Schichten 
einen  ov-Wert  von  0,061  34  ,  liegt  daher  erheblich  glinstiger. 

Das  Laminat  3  hat  einen  maximalen  °v-Wert  von  nur  0.04442. 

Dieser  Wert  liegt  im  Gegensatz  zum  Laminat  4  an  einer  Schicht- 
grenze,  wo  die  Fa  ser  r  i  c  h  tung  wechselt.  Das  konnte  einen  ungiin- 
stigen  EinfluB  haben.  Da  andererseits  das  Laminat  4  wegen  der 
fehlenden  90°-Schicht  und  der  damit  sehr  geringen  Quersteifig- 
keit  weniger  praxisnah  ist,  wird  fur  die  Versuche  das  Laminat  3 
a  usgewa  hi t . 
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Faser- 

rtchtung 

033  [  kN/  mm  2  1 

Irel  Fehler) 

0 13  (  kN  /mm 2 1 

|  rel  Fehler ) 

Ov*sign(o3j)  Vo33J»  3(o|3)  2' 

0° 

0.00078  (0,2%) 

0.00592  (0,0%) 

0.01029 

0r 

0.03613  (0,4%) 

-0.06809 

+45° 

-0.07231  (0,2%) 

-0.12806 

0° 

0.00512  (0,0%) 

-0.04698  (0,0%) 

0.08154 

0° 

-0.07202  (0,2%) 

-0.12555 

-45° 

HBflBIlillM 

0.03351  (0,4%) 

-0.05804 

0° 

0.00625  (0,0%) 

0.10294 

90° 

Bn 

0  (-) 

0.12353 

90° 

Tabelle  1  :  Verg 1 e i chswer te  ov  fur  das  Laminat  1  infolge 

Axiallast  o ^  =  1  KN/mirr 


Faser- 

richtung 

0  33  [kN  /mm  2] 

(rel.  Fehler) 

0,, !  kN  /mm  2 ) 

(rel  Fehler) 

Ov=  sign  (033)-\/o332*  3  (Oi3)2 

0° 

0.00127  (1,2%) 

0.00104  (0,0%) 

0.00221 

0° 

0.04036  (2,3%) 

0.00233  (0,0%) 

0.04056 

_j 

90° 

0.04373  (2,2%) 

0.00448  (0,0%) 

0.04442 

.  J 

0° 

-0.00962  (0,1%) 

0.00937  (0,0%) 

-0.01887 

0" 

-0.05945  (1,7%) 

0.04143  (0,6%) 

-0.09318 

*45° 

-0.06939  (1,5%) 

-0.08303  (0,2%) 

-0.15968 

0° 

-0.07541  (1,4%) 

-0.08400  (0,3%) 

-0.16387 

-45° 

mm 

0  (  -  ) 

-0.08455 

-45° 

Tabelle  2  :  Vergl eichswerte  ov_  fur  das  Laminat  3 
infolge  Axiallast  =  1  KN/mm*- 


Faser- 
ri rhFung 

0,,  (kN/mm  2 1 

( re l  Fehler ) 

0.,(  kN/mm2  ) 

( rel .  Fehler] 

ov=  sign(o331^/o332*  3(0, 3)2’ 

0° 

0.00176  (0.4%) 

0.00565  (0.0%) 

0.00994 

0° 

-0.01819  (4.0%) 

0.03244  (0.4%) 

-0.05905 

+  45° 

-0.01701  (4.2%) 

-0.05770  (0.2%) 

-0.10138 

n° 

0.00658  (0.0%) 

-0.03521  (0.0%) 

0.06134 

0° 

-0.01396  (4.9%) 

-0.05745  (0.2%) 

-0.10048 

-45° 

-0.01179  (5.8%) 

0.03183  (0.4%) 

0.05637 

0° 

0.01274  (0.0%) 

0.00561  (0.0%) 

0.01602 

0° 

mm 

BB 

0.01528 

0° 

o 

Tabelle  3  :  Vergl ei chswerte  Ov_  fur  das  Laminat  4 
infolge  Axiallast  =  1  KN/mm? 


0FVLR131/R-0184 


-14- 


1 


-15- 


31  /R-  0187 


-16- 


i 


k 


T 


I  DFVLR  131/R-0189 


Bild  7:  Randspannungen  cf13  urid  a33  unter 
Axiallast  fur  das  Laminat : 

4  0°2.90°.  0°2.  -45°.  0°,  -45°  ]  s 
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